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I.  INTRODUCTION 

Electro-optical  and  infrared  sensors  detect  and  identify  targets  by  discrimination 
of  their  thermal  contrast  against  the  cluttered  background.  In  the  case  of  targets  at  sea, 
this  means  distinguishing  the  target  radiance  from  the  sea  or  sky  background  radiance, 
attenuated  by  transmission  through  the  atmosphere,  together  with  the  radiance  emitted  by 
the  atmosphere  along  the  path.  The  performance  of  electro-optical  sensors  is  greatly 
influenced  by  atmospheric  factors  related  to  the  marine  boundary  layer  and  the  scenario 
geometry.  Polarized  images  of  the  research  vessel  POINT  SUR  (R/V  POINT  SUR)  and 
of  the  Netherlands  Oceanographic  ship  HM  Tydeman  have  shown  that  use  of  a  polarizer 
can  improve  contrast  temperature  difference  AT  against  background.  [Ref.l]  What  action 
is  still  needed  is  a  comparison  of  predicted  detection  range  with  polarization  and  without 
polarization  to  see  if  the  increased  contrast  increases  the  detection/  recognition  range 
more  than  the  loss  of  energy  (photons)  decreases  it. 

This  thesis  addresses  the  prediction  of  the  detection/recognition  range  without 
polarization.  The  concept  of  the  approach  is  that  at  the  maximum  detection/recognition 
range  the  projected  apparent  target  background  temperature  difference  (ATapp)  must  equal 
the  apparent  contrast  temperature  difference  required  by  the  sensor.  The  sensor 
performance  is  described  by  the  parameter  Minimum  Detectable  Temperature  Difference 
(MDTD)  or  Minimum  Resolvable  Temperature  Difference  (MRTD)  as  appropriate.  The 
apparent  temperature  contrast  AT  available  at  given  range  is  computed  from  the  zero 
range  target  temperature  converted  to  apparent  blackbody  radiance,  and  background 
radiance  computed  using  the  SEARAD[Ref.  2]  atmospheric  propagation  code,  a  modified 
form  of  the  standard  MODTRAN  code.  SEARAD  was  also  used  to  adjust  the  radiance  for 
atmospheric  attenuation  and  path  radiance.  SEARAD  is  a  DOS -compatible  program 
developed  at  NCCOSC-NRaD  (Naval  Command,  Control  and  Ocean  Surveillance 
Center-Naval  Research  and  Development)  for  application  to  naval  problems. 


The  source  of  meteorological  parameters  for  the  predictions  was  based  on  the 
EOMET95  experiment  that  was  conducted  on  May  15  to  24,  1995  in  which  polarized 
image  data  was  recorded  from  Monterey  Bay  Aquarium  Research  Institute  at  Moss 
Landing,  CA.  The  Research  Vessel  POINT  SUR  (R/V  POINT  SUR)  was  used  as  a  target 
and  as  a  measurement  platform.  The  ship  skin  temperature  distribution  was  recorded 
every  20  seconds  throughout  the  experiment  using  a  set  of  15  thermistors  on  the  ship 
skin.  Meteorological  and  sea  surface  data  measured  on  board  the  POINT  SUR  were  used 
as  input  to  the  atmospheric  propagation  computations.  The  meteorological  parameters 
and  observation  range  selected  in  this  experiment  were  used  as  inputs  to  the  "SEARAD" 
code  to  calculate  the  atmospheric  infrared  (8  to  12  urn)  transmittance  and  path  radiance. 

After  accessing  the  transmittance  and  path  radiance  calculated  from  SEARAD, 
Planck's  radiation  law  was  applied  to  find  the  apparent  radiance  difference  and 
corresponding  temperature  difference  (AT)  at  the  detection  range.  Generic  common 
module  FLIR  parameters  tabulated  by  Shumaker,  Wood  and  Thacker  [Ref.  3:p.  8-59] 
were  used  to  estimate  MDTD/MRTD  appropriate  to  a  generic  "Marine  Patrol  Aircraft" 
FLIR  sensor,  and  the  conventional  Johnsion  criterion  used  to  deduce  the  equivalent 
operational  target  spatial  frequency. 

In  selecting  the  scenario  for  this  analysis  the  ship  location,  heading  and  aspect 
with  respect  to  the  polarizing  images  in  the  EOMET95  image  data  base  were  collated  and 
processed;  the  collated  data  are  listed  in  Appendix  H  and  Appendix  I  to  this  thesis,  for 
later  analysis.  After  a  brief  summary  of  the  relevant  radiation  propagation  theory,  the 
EOMET95  experimental  procedure  and  measurement  are  discussed,  including  the 
procedure  for  target  aspect  and  apparent  dimensions,  and  the  computation  and  modeling 
of  the  FLIR,  MRTD,  and  MDTD  functions  are  described.  Computation  of  target  contrast 
and  MRTD  and  MDTD  range  estimates  are  then  described,  for  selected  target-sensor 
scenarios  from  the  EOMET95  data  base  and  followed  with  discussion  of  the  results. 


II.  INFRARED  RADIATION  FUNDAMENTALS 

Optical  radiation  covers  the  ultraviolet  (UV),  visible,  and  infrared  (IR)  portions  of 
the  electromagnetic  spectrum.  The  ultraviolet  portion  ranges  from  about  0.1  to  0.38  urn. 
The  visible  portion  is  from  approximately  0.38  to  0.76  urn  in  wavelength.  The  infrared 
portion  is  divided  into  the  near  infrared  or  short-wavelength  infrared  (SWIR)  region 
(from  0.77  to  3  urn),  the  middle-wavelength  infrared  (MWIR)  region  (from  3  to  8  urn), 
the  long-wavelength  infrared  (LWIR)  region  (from  8-14  urn),  and  the  far  and  extreme 
infrared  region  (from  14  to  1000  urn).  The  radiant  energy  received  by  an  optical  system 
from  a  scene  object  can  be  generally  broken  down  into  two  components:  (1)  radiant 
energy  due  to  the  self-emission  of  the  object  and  (2)  radiant  energy  reflected  by  the  object 
due  to  external  radiation  sources.  The  objects  being  observed  by  the  optical  system  can  be 
categorized  as  either  the  target  or  the  background.  [Ref.  4:p.  37]  The  radiation  emitted  and 
reflected  from  the  target  and  the  background  traverse  through  the  atmospheric  medium, 
where  absorption  and  scattering  take  place  by  molecular  constituents  of  the  atmosphere 
and  the  aerosol  particles  suspended  in  it.  Re-emission  and  scattering  from  these 
atmospheric  components  also  contribute  path  radiance  to  the  received  power.  The 
absorption  ,  scattering,  and  emission  are  conventionally  computed  using  standard  DOD 
propagation  computer  codes.  For  this  study  the  developmental  SEARAD  code  was  used. 
SEARAD  is  described  further  in  Appendix  D. 


A.         THEORY  OF  BLACKBODY  RADIATION 

This  section  summarizes  commonly  known  material  which  can  be  found  in  many 
tutorial  sources,  for  example  "  The  Infrared  Handbook"  [Ref.  5].  All  objects  with  a 
temperature  above  absolute  zero  emit  radiation.  Relative  spectral  intensity  of  the  radiation 
is  dependent  primarily  on  the  temperature  of  the  object  and  on  the  radiant  properties  of 
the  material   the   object  is  made  of  (in  particular,   the   spectral   emissivity   of  the 


material). [Ref.  6:p.  24]  A  blackbody  is  defined  as  an  ideal  body  or  surface  that  absorbs  all 
radiant  energy  incident  upon  it  at  any  wavelength  and  at  any  angle  of  incidence,  so  that 
none  of  the  radiant  energy  is  reflected  or  transmitted.  Therefore,  the  term  blackbody 
refers  to  a  perfect  absorber  of  radiation.  Because  all  radiant  energy  is  absorbed  by  a 
blackbody  and  none  is  reflected  or  transmitted,  it  will  appear  to  be  "black". [Ref.  4:p.  44] 

1.  Planck's  Radiation  Law 

Planck's  radiation  law  can  be  used  to  calculate  the  spectral  radiance  of  a 
blackbody  source.  Planck's  equation  satisfies  thermodynamic  requirements  and  gives  the 
spectral  radiant  exitance  Mbb(v,T)  of  a  blackbody  in  terms  of  temperature,  T,  and 
radiation  frequency,  v: 


2nh  v5 
^(7^-1) 


—  Jit l   V 


where 

h  is  Planck's  constant, 

c  is  the  speed  of  light, 

k  is  Boltzmann's  constant  (1.38054xl0"23  J-K/1) 

This  equation  can  also  be  written  in  terms  of  wavelength,^  : 

2;t/zc2 
A  (e  - 1) 

Equation  2.1  and  2.2  express  the  spectral  radiant  exitance  of  a  source  radiating  in  a 

vacuum.  For  any  other  media,  c  should  be  replaced  by  the  light  velocity  in  the  media,  c/n, 

where  n  is  the  refractive  index  of  that  medium. 

2.  Emissivity 

The  radiation  from  all  bodies  in  a  wavelength  range  under  consideration  is 
dependent  firstly  on  the  temperature  and  secondly,  in  the  case  of  non-black  bodies,  on  the 
material  composition  and  surface  condition.  The  ratio  of  the  emitted  radiation  of  any 


given  temperature  radiator  to  the  emitted  radiation  of  a  black  body  of  the  same 
temperature  is  known  as  the  emissivity.[Ref.  6:p.  50] 

M=eMBB        ,  2.3 

3.  Stefan-Boltzmann  Law 

The  total  radiated  energy  from  a  unit  area  of  a  surface  in  a  unit  of  time  over  a  2n 
steradian  solid  angle  is  called  the  radiant  exitance  of  that  surface.  For  blackbody  source 
the  radiance  is  proportional  to  the  fourth  power  of  the  absolute  temperature  of  the  source. 
This  gives  [Ref.  4:p.45] 

M(T)=ct  T4  ,  2.4 

where  a  is  the  Stefan-Boltzmann  constant  =5.6697xl0'12  W-cm"2K"4 
For  use  with  non-black-body  sources,  this  law  is  modified  by  the  inclusion  of  an 
"effective  emissivity",  8,  giving  the  form 

M(T)=eo  T4 ,  2.5 

where  e<l.  For  in-band  radiance,  the  radiance  N  can  be  calculated  using  the  universal 
blackbody  curve  or  the  proportional  radiation  table  (Appendix  B)  with  Equation  2.6 

soT4 

N  = Aq  ,  2.6 

where 

£  =  emissivity  of  the  source. 
Aq  =  difference  of  two  q  values  from  the  proportional  radiation  table  (Appendix  B) 

The  proportional  radiation  q  =  f  (X,T)  represents  the  fraction  of  the  radiant  exitance 

emitted  by  a  blackbody  at  temperature  T  at  all  wavelengths  up  to  the  selected  value  of  X. 

It  can  be  described  in  Equation  2.7 

?M(X,f)dA 

q  =  — ,  2.7 

^  M(X,T)dA 


If  q  values  are  not  shown  on  the  table,  they  can  be  calculated  by  interpolation  between 
two  q  value  entries. 

4.  Kirchhoff  s  Law 

When  radiation  is  incident  upon  a  body,  some  of  it  is  transmitted,  some  absorbed, 
and  some  is  reflected.  Thus,  the  ratios  of  each  of  these  to  the  incident  power  must  add  up 
to  unity  [Ref.  1 1  :p.  1  -29]  given 

a(X)+p(X)+x(X)=l  ,  2.8 

where 

a(A,)=  spectral  absorptivity 
p(A,)=  spectral  reflectivity 
x(X)=  spectral  transmissivity 

5.  Lambert-Beer  Law 

Atmospheric  transmittance  is  a  function  of  wavelength,  which  can  be  described  by 
the  Lambert-Beer  law 

x(X)=e^  ,  2.9 

where  x(X)  =  spectral  atmospheric  transmittance, 
X  =  wavelength, 
u.  =  extinction  coefficient, 
R  =  path  length. 
The  extinction  coefficient  depends  on  the  atmospheric  molecular  composition  and 
aerosol  concentration. 

B.    INFRARED  RADIATION  SOURCES 

Optical  sensors  function  by  sensing  the  radiation  emitted  by  an  object  or  the 
radiation  reflected  from  the  object  being  illuminated  by  an  external  radiation  source. 
Detection  of  a  jet  engine  plume  by  an  infrared  detector  is  an  example  of  the  first  category 


whereas  the  visible  detection  of  an  area  illuminated  by  the  sun  is  an  example  of  the 
second  category.  In  both  cases  these  sensors  are  passive,  as  they  do  not  generate  their 
own  power.  Active  sensors  on  the  other  hand  generate  their  own  power  and  detect  targets 
from  the  reflected  energy  of  this  power.  Electro-optical  systems  may  operate  using 
natural  and  artificial  sources  of  power. 

C.         ATMOSPHERIC  PROPAGATION  OF  INFRARED  RADIATION 

The  attenuation  of  optical  radiation  due  to  the  atmosphere  arises  from  the 
individual  or  collective  effects  of  the  following  phenomena 

1 .  Molecular  absorption 

2.  Molecular  scattering 

3.  Aerosol  absorption 

4.  Aerosol  scattering 

The  two  main  causes  of  attenuation,  however,  are  molecular  absorption  by  several  minor 
constituents  of  the  atmosphere  and  scattering  due  to  atmospheric  molecules  and  the 
particles  in  the  present  atmosphere  (aerosols).  Molecular  absorption  occurs  mainly  in 
several  more  or  less  narrow  absorption  bands  and  is  due  to  transmittence  in  certain 
molecules,  matched  to  that  transition  from  one  state  of  vibration  and  rotation  to  another, 
thereby  absorbing  or  emitting  a  photon. [Ref.  4:p.  75] 

1.  Absorption  Effects 

Radiation  absorption  due  to  molecules  in  a  gaseous  medium  can  be  divided  into 
two  processes: 

( 1 )  Atomic  absorption  is  due  to  the  transition  of  electrons  in  an  atom,  and  usually 
requires  a  photon  with  an  energy  of  a  few  electron  volts.  It  generally  takes  place  in  the 
visible  or  ultraviolet  region  of  the  spectrum. 


(2)  Molecular  absorption  is  due  to  a  transition  between  electronic,  vibrational,  or 
rotational  energy  states  of  a  gas  molecule  or  due  to  a  combination  of  all  of  these.  This 
mechanism  is  dominant  in  the  infrared  range. [Ref.  4:p.  76] 

2.  Scattering  Effects 

Scattering  is  the  process  by  which  the  energy  in  an  electromagnetic  wave  is 
intercepted  and  reradiated  into  An  steradians  solid  angle.  It  results  from  the  interaction  of 
the  wave  field  with  the  electron  oscillators  in  the  scattering  medium.  These  are  excited  by 
the  incoming  wave  field,  behaving  as  forced  harmonic  oscillators  which  reradiate  at  the 
frequency  of  the  incident  wave.  The  secondary  wave,  however,  shows  a  phase  lag  from 
the  primary  dependent  on  the  difference  between  the  wave  frequency  and  the  oscillator 
resonance  frequency.  Any  inhomogeneity  in  refractive  index  can  cause  scattering;  in  the 
atmosphere  the  component  gas  molecules,  aerosol  particles  of  various  sizes,  and  fog,  rain 
and  hail  drops  are  all  effective  scatters.  The  spatial  distribution  of  the  scattered  radiation 
is  strongly  dependent  on  the  relative  magnitudes  of  the  particle  size  and  the  wavelength. 
For  particles  very  small  compared  with  the  wavelength,  the  scattering  is  approximately 
isotropic;  as  the  ratio  of  size  to  wavelength  increases,  the  scattering  is  concentrated  more 
into  the  forward  hemisphere.  For  very  large  scattering  objects  forward  scattering 
dominates,  and  secondary  lobes  develop  at  other  angles  in  the  radiation  pattern.  In  the 
small  particle  regime,  the  scattering  can  be  described  by  the  relatively  simple  theory 
developed  by  Rayleigh,  which  predicts  scattering  proportional  to  the  second  power  of  the 
volume  of  the  wavelength.  This  describes  scattering  by  atmospheric  molecules.  For 
particle  sizes  greater  than  one  tenth  of  the  wavelength,  the  Rayleigh  approximation  is  no 
longer  adequate;  the  more  comprehensive  Mie  theory  (1908)  is  required.  [Ref.  7:p.  12-22] 

a.  Rayleigh  Scattering 

When  the  atmosphere  does  not  contain  large  size  particles,  but  only  the 
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molecules  of  its  gaseous  constituents,  the  scattering  phenomenon  is  dominated  primarily 
by  these  molecules.  The  intensity  of  the  scattering  depends  on  the  number  and  index  of 
refraction  of  these  gaseous  molecules  and  on  the  scattering  angle,  wavelength,  and 
polarization  of  the  incident  radiation.  The  attenuation  of  a  light  beam  due  to  Rayleigh 
scattering  is  expressed  by 

I=I0exp(-CTmR)  ,  2.10 

where 

I0  is  the  magnitude  of  the  original  radiation 
I  is  the  magnitude  of  the  scattered  radiation 
am  is  the  Rayleigh  scattering  coefficient 
R  is  the  path  length 
The  Rayleigh  scattering  coefficient  am  is  the  product  of  two  parameters 

<V=<*rN   ,  2.11 

where 

aT  is  the  Rayleigh  scattering  cross  section, 
N  is  the  molecular  density. 

The  Rayleigh  scattering  cross  section,  ar  ,  is  defined  as  the  cross  section  of  an  incident 
wave,  so  that  the  total  power  scattered  by  a  gaseous  molecule  in  all  directions  is  equal  to 
the  power  flowing  across  the  cross  section.  It  is  given  by  the  expression 

8;r3(»2-l)2 
ff'=      3N>Z*       '  2-12 

where 

n  is  the  index  of  refraction  of  the  gaseous  medium 

N  is  the  molecular  number  density. 

Note  that  the  scattering  cross  section  ar  is  a  function  of  wavelength  and  proportional  to 

A."4  .  It  therefore  dominates  more  in  the  ultraviolet  and  visible  regions  than  in  the  infrared 

region.  The  formula  also  indicates  that  ar  is  a  function  of  the  index  of  refraction  of  air  and 

molecular  number  density.  [Ref.  4: p.  78] 


b.         Mie  Scattering 

Aerosol  scattering,  or  Mie  scattering,  is  concerned  with  the  radiation 
interaction  with  spherical  particles  of  various  sizes.  The  transmitted  beam  irradiance  here, 
similar  to  Rayleigh  scattering,  is  expressed  by 

I=I0exp(-aaR) ,  2.13 

where  I,  I0,  and  R  are  as  defined  previously,  and  aa  is  the  aerosol  attenuation  coefficient, 
which  is  a  function  of  aerosol  density.  Aerosol  size  distribution  does  not  change  much 
with  altitude,  whereas  the  aerosol  density  does.  The  relationship  for  aerosol  attenuation 
coefficient  at  altitude  h  is  given  by 

M{h) 

Where  M(h)  denotes  the  aerosol  number  density  at  height  h  and  aa(0)  is  the  aerosol 
attenuation  coefficient  at  sea  level. [Ref.  4:p.  81]  These  absorption  and  scattering  losses 
must  be  compounded  when  treating  propagation  through  the  atmosphere. 
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III.       EOMET95  EXPERIMENT  AND  DATA  COLLECTION 

A.         INTRODUCTION 

The  data  for  this  thesis  were  obtained  from  an  ocean  field  experiment  conducted 
from  15  to  24  May  1995  off  the  coast  of  Monterey.  The  EOMET  95-Spring  Cruise 
involved  the  R/V  POINT  SUR,  which  supported  shipboard  measurements  of  atmospheric 
and  oceanic  planetary  boundary  layer  data  made  by  the  NPS  Boundary  Layer 
Meteorology  group.  In  addition,  measurements  from  stationary  platforms  located  in 
Monterey  Bay  and  its  shoreline  were  collected.  The  R/V  POINT  SUR  is  a  135  foot  ship 
(dimension:  41.5Lx9.75Wx8.8H  (m))  owned  by  the  National  Science  Foundation  and 
operated  by  Moss  Landing  Marine  Laboratory.  Of  particular  interest  was  the  availability 
of  the  POINT  SUR  during  the  mornings  of  the  cruise  for  imaging  operations,  which 
placed  the  ship  at  various  ranges,  and  aspect  angles,  which  enabled  the  collection  of 
polarized  IR  imager  data  coincident  with  the  environmental  measurements.  During  these 
operations,  measurements  were  made  of  the  degree  of  polarization  and  contrast 
improvement  factor  for  ship  images  as  a  function  of  range  and  aspect  angle.  The 
environment  impacts  these  measurements  in  the  following  ways: 

(l)The  observed  wave  spectrum  and  wind  speed/direction  affects  the  sea  surface 
degree  of  polarization. 

(2)The  SST  (Sea  Surface  Temperature)  and  air  temperature  impacts  the  target  to 
background  contrast. 

(3)The  air  temperature,  humidity,  and  aerosol  content  affect  the  IR  path  radiance 
and  extinction  coefficient. 


11 


B.         DATA  MEASUREMENTS 

1.  Imaging  System 

The  NPS  AGA  780  Thermovision  Thermal  Imaging  System  was  set  up  and 
operated  from  Moss  Landing.  The  site  was  the  MBARI  Building  located  at  approximately 
6  meters  above  sea  level.  The  estimated  horizon  range  from  this  location  was  4.5  nautical 
miles.  The  AGA  780  imaging  system  is  a  dual  band  serial  scanning  infrared  imager 
capable  of  operation  in  the  3-5  and  8-14  micron  bands.  The  system  utilizes  a  7x7  degree 
field  of  view  lens  or  an  optional  3.5  degree  lens  for  each  channel  to  focus  incoming 
infrared  radiation  through  a  polarizing  filter  wheel,  then  through  a  dual  rotating  prism 
system  onto  the  corresponding  detector.  The  detectors  are  cooled  to  77  K  by  dewars  filled 
with  liquid  nitrogen.  The  filter  wheel  was  rotated  to  provide  vertical,  horizontal  and  zero 
polarization  settings  during  data  collection.  The  detectors  produce  a  current,  which  is 
called  "thermal  value"  and  is  proportional  to  the  intensity  of  the  radiation  received.  The 
signal  is  then  amplified  and  converted  into  a  video  signal  for  monitor  display. 

2.  Calibration  Technique 

The  AGA  780  imaging  system  was  calibrated  with  a  laboratory  blackbody  source. 
The  system  output  or  "thermal  value  "  is  recorded  against  the  blackbody  temperature 
through  a  range  of  temperatures.  A  calibration  curve  is  generated  which  takes  into 
account  the  scanner,  lenses  and  filters  (Moretz,  1994).  All  the  data  collected  during 
EOMET  95  Spring  Cruise  was  either  calibrated  or  can  be  calibrated. 

3.  Target  Temperature 

To  provide  a  means  of  determining  target  temperature,  the  R/V  POINT  SUR  was 
instrumented  with  thermistor  temperature  sensors  at  numerous  points  on  the  ship's  skin. 
Each  thermistor  was  fastened  to  the  skin  by  use  of  beryllium  adhesive,  which  is 
extremely  conductive.  The  thermistor  was  then  covered  with  white  epoxy  to  assure 
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maximum  albedo  from  sunlight.  The  sensors  were  linked  by  shielded  cable  to  a  data 
collection  computer  in  the  instrumentation  spaces  aboard  ship. 

4.     Visual  Imaging 

Coincident  images  in  the  visual  range  were  made  from  the  IR  imager  site  at  the 
MBARI  building.  These  were  made  with  a  high  quality  video  recorder  system.  The 
polarization  analysis  of  these  data  will  be  presented  in  a  further  report  or  thesis.  [Ref.8] 
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IV.     MRTD  AND  MDTD 

Minimum  Detectable  Temperature  Difference  (MDTD)  and  Minimum  Resolvable 
Temperature  Difference  (MRTD)  are  two  standard  descriptions  of  FLIR  performance. 
They  are  defined  as  those  temperature  contrast  differences  which  make  a  standard  bar 
chart  just  resolvable  or  a  square  target  just  detectable,  respectively,  to  a  trained  observer. 
Each  of  these  provides  a  sensitivity  measurement  of  the  FLIR  system.  Both  of  them 
employ  the  concept  of  the  Equivalent  Blackbody  Temperature  of  a  non-blackbody  object. 
The  equivalent  blackbody  temperature  is  that  temperature  a  blackbody  must  have  to  emit 
the  same  radiance  as  the  target.  This  allows  the  use  of  a  single  value  to  represent  the 
target's  temperature  and  emissivity.  To  complete  the  comparison  of  apparent  AT  with  the 
system  performance  AT  requires  numerical  values  of  these  quantities  as  functions  of 
spatial  frequency  and  range. 

A.         MINIMUM  RESOLVABLE  TEMPERATURE  DIFFERENCE 

MRTD  (often  also  abbreviated  MRT)  is  a  function  of  spatial  frequency.  It  gives 
the  AT  between  the  hot  (cold)  and  ambient  temperature  bars  of  a  standard  4  bar  (7:1 
aspect  ratio)  chart  required  to  make  the  bars  just  resolvable  as  a  function  of  the  spatial 
frequency  of  the  bars.  MRTD  is  a  measure  of  the  performance  of  the  entire  FLIR  system. 
Since  the  MRTD  includes  the  observer,  the  measurement  is  subjective.  The  computed 
system  MRTD  shows  an  asymptotic  value  beyond  which  the  detector  angular  subtense  is 
greater  than  the  reciprocal  of  the  spatial  frequency,  or  the  "spatial  period"  and  finer  detail 
cannot  be  distinguished.  For  a  typical  FLIR  (Common  Module),  the  instantaneous  field  of 
view  (IFOV)  is  approximately  0.25  mRad  so  that  the  in-scan  detector  subtense  Ax  is 
approximately  0.25  mRad  and  the  cut  off  spatial  frequency  (asymptote)  is  approximately 
l/Ax=l/0.25=4  cyc/mRad.  If  we  use  the  "rule  of  thumb"  asymptote  between  0.7/Ax  and 
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0.9 /Ax,  we  would  find  that  the  asymptote  is  approximately  3  cyc/mRad.  The  effective 
IFOV  depends  on  the  magnification  of  the  optics  that  may  be  variable  in  some  FLIRs. 
Shumaker[Ref.  3:p.  8-52]  gives  a  typical  form  for  the  MRTD. 


MRT(v)  =  2-™RT 


v  AxAv 


-,!/ 


[teFrXo,Nss]-1A  ,  4.1 


MTFs(v) 

where  SNRT  =  signal  to  noise  ratio  threshold, 

NET  =  noise  equivalent  temperature, 

px  =  noise  filter  factor, 

Ax  =  in-scan  detector  subtense  (mRad), 

Ay  =  cross-scan  detector  subtense  (mRad), 

MTF  =  modulation  transfer  function, 

v  =  spatial  frequency  for  which  MRT  is  desired  in  cyc/mRad, 

L  =  length-to-  width  ratio  of  the  bar  (7), 

te  =  eye  integration  time  in  seconds, 

Fr  =  frame  rate  s"1, 

Nos  =  overscan  ratio, 

Nss  =  serial  scan  ratio. 
Shumaker  [Ref.3:p.  8-7]  also  provides  an  expression  for  the  Noise  Equivalent 
Temperature  Difference  (NETD),  the  target  apparent  contrast  AT  which  yields  an  output 
SNRofl,as 


10(FOVxFOV,FrNosNss) 


NET    = : 4  2 

l^ND77sc)  2DAxA^D**7cs      /^T^o 

where  FOVx  =  in-scan  FOV  in  mRad, 

FOVy  =  cross-scan  FOV  in  mRad, 

ND  =  number  of  detectors, 

r|  sc  =  scan  efficiency, 

3N/3T  =  derivative  of  Planck's  equation  in  W/cm2  -K-Sr, 
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D**       =  the  band  average  detectivity  in  cm  Hz1/2  /W, 
D       =  aperture  diameter  in  meters, 
x0        =  transmission  of  the  optics. 


B.         MINIMUM  DETECTABLE  TEMPERATURE  DIFFERENCE 

The  MDTD  of  a  FLIR  system  gives  the  temperature  difference  between  an 
isolated  square  target  and  a  uniform  background  that  renders  the  square  just  detectable,  as 
a  function  of  the  size  of  the  square.  In  order  to  maintain  similarity  between  graphs  of 
MDTD  and  MRTD,  MDTD  is  plotted  as  a  function  of  the  spatial  frequency  and  range. 
Shumaker  [Ref.3:p.  8-67]  provides  a  typical  form  for  MDTD 


(SNRT)  (NET)(  Qr  +  r2  VAxAy)^ 
MDT{v)=V- ii A  A J—       ,  4.3 


QT 


-(r*+r2B+QT)teFrNosNs 


4 


where  QT  =  target  angular  subtense  (mRad2), 

rs  =  system  resolution  (mRad), 

rB  =  back-end  resolution  (mRad). 
The  other  definitions  of  the  parameters  in  Equation  4.1  such  as  SNRT,  NET,  Ax, 
Ay,  te ,  Fr  Nos,  Nss  are  the  same. 

C.         MRTD  AND  MDTD  COMPUTATION 

Shumaker[Ref.3:p.8-59]  provides  the  following  set  of  parameters  that  are 
postulated  to  be  appropriate  to  the  Common  Module  FLIR  in  NFOV  (Narrow  Field  of 
View)  mode.  px=(l+(2v»  r^)'^'^  is  computed  from  the  equivalent  back-end  resolution 

rB  which  is  given  for  this  example  as  0.335  mRad.[Ref.  3:p.  8-60]  The  resulting 
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computation  from  MATLAB  is  shown  in  Figure  5  which  shows  MRTD  and  MDTD  vs 
spatial  frequency.  MATLAB  code  is  shown  in  Appendix  J  section  B. 


FOVx 

6.86° 

FOVY 

5.16° 

L 

7 

Fr 

30  sec"1 

Nos 

1 

rE 

0.3  mRad 

Nss 

1 

ND 

180 

rd 

0.08 

*lsc 

0.75 

T0 

0.70 

Ns 

2 

K 

0.1  sec 

SNRT 

2.5 

MTF 

Tabulated 

f 

0.20  m 

Ay 

0.25mRad 

D* 

4xl010 

D 

0.10m 

Ax 

0.25  mRad 

NET 

0.1  °K 

Table  1.  Parameters  of  NFOV  FLIR 


Spatial  Frequency  (  v  ) 
(cyc/mRad) 

Ax=  0.25  mRad 
MTF 

0.00 

1.00 

0.25 

0.89 

0.50 

0.79 

0.75 

0.64 

1.00 

0.53 

1.25 

0.43 

1.50 

0.35 

1.75 

0.27 

2.00 

0.20 

2.50 

0.10 

3.00 

0.04 

Table  2.  MTF  for  NFOV 
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V.  DATA  COMPUTATIONS  AND  ANALYSIS 

A.         APPARENT  AT  COMPUTATION 

1.  Data  Summary  of  Infrared  Measurements  in  EOMET95 

The  intent  of  the  infrared  imaging  measurements  made  during  the  EOMET95 
cruise  was  to  observe  the  variation  with  range  of  the  sea  surface  degree  of  polarization 
and  of  the  apparent  target/background  temperature  contrast  improvement  with 
polarization  filtering  for  the  R/V  POINT  SUR,  and  to  attempt  to  observe  emission 
polarization  effects  on  facets  of  the  ship  image.  To  attain  these  ends  the  available  ship 
operation  periods  were  allotted  alternately  to  two  series  of  measurements; 

(1)  Rotation  of  the  POINT  SUR  by  small  increments  of  heading  angle  through 
360°  at  fixed  range  and  bearing  from  the  sensor,  for  a  selection  of  ranges  and  bearings. 
For  each  target  aspect  and  range  images  were  recorded  with  vertical  and  horizontal 
polarizing  filter  and  no  polarizer. 

(2)  Recording  of  vertical  and  horizontally  polarized  and  unpolarized  images  of  the 
POINT  SUR  at  selected  aspect  at  a  set  of  selected  ranges  along  a  constant  bearing  line  to 
a  maximum  observation  range. 

Reference  ship  skin  and  sea  and  air  temperature  measurements  were  recorded 
during  the  entire  observation  period,  as  described  below.  Global  Positional  System 
(GPS)  data  for  ship  heading  and  position  and  radiosonde  vertical  atmospheric  profiles 
were  recorded  on  the  R/V.  POINT  SUR.  These  values  were  accessed  with  the  assistance 
of  the  NPS  Boundary  Layer  Meteorology  Group.  The  appropriate  data  selections  are 
included  in  the  appendices  to  this  thesis. 

For  the  purposes  of  this  thesis  certain  data  points  were  selected  from  this 
extensive  data  base,  for  calculation  of  predicted  detection  or  recognition  range. 
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2.  Ship  Temperature  Computation 

Ship  skin  temperature  samples  were  collected  from  15  thermistors  attached  with 
epoxy  cement  to  the  surfaces  of  the  larger  facets  of  the  ship,  located  and  distributed  as 
indicated  on  Figure  12  and  13.  The  thermistor  voltages  were  recorded  continuously  at  20 
second  intervals.  The  records  were  listed  by  time  entries  in  seconds  from  midnight 
GMT  for  each  day.  These  temperature  files  are  listed  in  Appendix  A.  The  format  of  the 
output  files  is  illustrated  by  the  sample  file  of  Table  3.  In  each  file  the  first  header  line 
gives  the  file  name.  The  second  line  gives  the  Month,  Day  and  Year  from  the  computer 
clock,  followed  by  the  Hour,  Minute  and  Second  of  the  file  on  the  next  line.  The  next  line 
gives  the  DAC  voltage  (digital  value),  start  channel,  stop  channel,  number  of  samples 
averaged,  and  reference  channel  number  on  each  board.  Note  that  thermistor  channels  1 
and  9  always  read  25.00  °  C  since  they  are  reference  channels.  The  position  data  used  in 
this  analysis  (extracted  from  the  ship  GPS  data  file  listed  in  Appendix  H  for  the  period  of 
IR  measurement)  are  for  the  period  from  approximately  1830  to  1840  May  16  1995 
GMT. 

May  16.1 

Mon,Day,Year=  5  16  1995 
Hour,Min,Sec=  8  27  36 
volts,istart,istop,nsam,iref  1.49939  1  16  100  1 

1326456  25.001  12.856  13.403  11.963  14.615  11.677  12.359  11.634  25.001  14.242 
22.956  12.010  11.617  11.431  15.410  25.874 

1326476  25.001  12.864  13.397  11.959  14.646  11.673  12.371  11.633  25.001  14.200 
22.937  12.001  11.607  11.426  15.229  25.535 

Table  3.  File  Identification  for  Ship  Temperature  Thermistor  Reading 

These  thermistor  readings  provide  ship  temperature  at  selected  points;  from  these  an 
average  ship  temperature  must  be  derived  as  input  to  a  range  estimation.  This  has  been 
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done  by  comparing  previous  detailed  radiometric  measurements  made  of  the  temperature 
distribution  over  images  of  the  R/V  POINT  SUR  [Ref.9]  with  the  concurrent  thermistor 
spot  temperature.  From  these  data  each  facet  area  on  the  ship  profile  has  been  correlated 
with  a  representative  thermistor  location,  as  shown  in  Figure  13.  Based  on  this 
correlation,  an  area  weighting  factor  is  assigned  to  each  recorded  thermistor  temperature 
and  a  weighted  average  temperature  computed  for  the  appropriate  ship  aspect.  While  it  is 
realized  that  the  temperature  distribution  will  change  with  change  of  ship  operating 
conditions  and  meteorological  conditions,  this  gives  a  relatively  simple  technique  for 
modeling  the  ship  signature,  roughly  comparable  to  methods  previously  used  in  Tactical 
Decision  Aid  models,  and  has  been  used  here  for  the  input  ship  temperature  for  the 
sample  case.  Table  7  shows  ship  average  temperature  for  computation  by  equation  5.1. 


Average  temperature  = '■ — '- ,  5.1 


3.  Ship  GPS  Position  and  Heading  Data  Correlation 

The  times  of  recording  of  the  ship  GPS  data  and  the  skin  temperature  data  are  not 
synchronized.  The  temperature  is  recorded  every  20  seconds  while  the  position  is 
recorded  every  10  minutes.  For  accuracy,  in  target  signature  modeling,  recorded 
simultaneous  ship  heading  and  position  values  are  selected  from  the  GPS  files  (Appendix 
H),  and  the  appropriate  temperature  values  interpolated  from  the  temperature  files  ( 
Appendix  A).  As  an  example,  to  compute  a  detection/recognition  range  for  the  target  at 
heading  120  degrees,  the  ship  location  and  atmospheric  parameters  are  interpolated 
between  the  measurement  points  in  Appendix  H  closest  to  that  heading,  those  entries 
highlighted  in  Appendix  H.  Data  were  recorded  at  1830  GMT  at  heading  172  degrees, 
and  1840  GMT  at  106  degrees.  Interpolating  to  120  degrees  assuming  a  constant  rotation 
rate  of  66  degrees  in  10  minutes,  or  0.1 1  degrees  per  second,  gives  us  a  time  of  18:37:53 
GMT  or  11:37:53  local  time.  To  read  the  corresponding  skin  temperatures  we  convert  this 
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time  into  seconds  after  the  start  of  file  at  08:27:36  at  file  number  1326456.  This  lapse 
time  of  1 1416  seconds  corresponds  to  file  number  1337872.  The  closest  skin  temperature 
set  to  this  is  file  number  1337876,  shown  in  the  sample  listing  below,  highlighted  in 
Appendix  A,  Page  65.  It  is  assumed  that  no  significant  changes  in  meteorological 
parameters  will  occur  in  the  4  seconds  discrepancy,  so  that  the  data  selected  for  this  time 
will  form  a  consistent  set. 

1337856  25.001  16.392  17.862  17.334  15.415  15.885  14.338  15.873  25.001  18.294 
24.326  15.281  16.169  15.352  22.299  28.721 

1337876  25.001  16.423  17.864  17.291  15.361  15.830  14.368  15.909  25.001  18.333 
24.164  15.296  16.151  15.346  22.393  28.829 

1337896  25.001   16.458  17.844  17.241   15.295  15.771   14.383  15.946  25.001   18.366 
24.093  15.304  16.117  15.305  22.319  28.933 

Table  4.  Thermistor  Sample  Output  File 


Side 

S  Pixel 

I  Temperature 

ETemperature  /ZPixel 
=  Average  Temperature 

Port 

59 

976.196 

16.54569 

Starboard 

59 

1028.287 

17.42859 

Total  Average 
Temperature 

118 

2004.483 

16.98714 

Table  5.  Ship  Average  Temperature  vs  Pixel  for  the  Sample  Case  Selected 

4.  Temperature  and  Radiance  Conversion 

The  average  of  ship  temperature  computed  for  this  sample  period  is  273  + 
16.987144  =  289.987144  K.  Ship  temperature  and  sea  background  temperature  were 
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converted  to  in-band  radiance  by  application  of  the  Planck  function.  The  ship 
temperature  is  then  converted  to  in-band  radiance  Ns  at  zero  range 

Ns  =  eaT4Aq/7t 

=  0.95x5.6697  xl0"8x(289.98714)4x0.24861  17/tt  =  30.14347844  (W/m2Sr) 
where 

8=0.95,  emissivity  of  R/V  POINT  SUR  and  Aq=q2-ql  is  the  fraction  of  the 
thermally  emitted  radiance  lying  in  the  band  from  A,,  to  X2  at  source  temperature  T. 

XJ  =  8^x289.98714  =  0.23198971  cmxK  =>ql=  0.1314394 

X2T  =  12^x289.98714  =  0.3479845  cmxK  =>ql=  0.380051 1 

and  Aq=q2-ql  =0.2486117 

The  q  values  have  5  digits  as  reported  in  the  radiation  table  (Appendix  B)  to  keep 
the  fraction  of  the  radiant  exitance  as  close  as  possible  to  the  integral  of  Planck's  Law. 

5.  Transmittance  and  Path  Radiance 

The  SEARAD  program  was  used  to  obtain  the  transmittance  and  path  radiance 
between  the  target  and  the  FLIR  sensor.  Appendix  D  describes  how  to  install  SEARAD 
and  to  run  the  example  case.  There  are  some  points  which  need  to  be  noted  carefully 
when  inserting  radiosonde  profile  data  in  the  input  data  file.  The  maximum  number  of 
atmospheric  levels  that  can  be  inserted  is  34.  The  position  of  every  data  point  is  important 
for  SEARAD,  so  it  is  essential  that  each  must  be  inserted  in  the  correct  column  and  row 
of  the  input  table.  The  Windows  text  editor  may  cause  a  problem  after  the  data  have  been 
inserted  and  saved.  It  is  recommended  that  we  use  the  simplest  editor  available  (i.e.  the 
MS-DOS  editor).  To  run  SEARAD  it  is  not  necessary  to  follow  all  the  steps  listed  in 
Appendix  D.  "Tape  5"  is  the  data  file  which  will  be  read  into  SEARAD  and  executed; 
thus  we  need  only  modify  the  "Tape  5"  data.  The  necessary  steps  to  run  SEARAD  are 
shown  below: 

•  Go  to  DOS  program;  it  should  show  C:WINDOWS>  . 

•  Under  C:WINDOWS  type  "cd  c:\". 

23 


When  it  shows  C:>  then  type  "cd  searad". 
It  shows  C:\searad>  then  type  "edit  tape5". 
It  shows  MS-DOS  editor. 
Modify  the  data  and  save  it;  then  exit. 
It  should  return  to  C:\searad>. 

Under  C:\searad>  type  "searad"  which  means  to  execute  SEARAD. 
We  can  use  WORD  or  WORDPERFECT  to  read  the  output  file  "out"  under 
the  SEARAD  directory. 
The  inputs  of  the  SEARAD  code  are  tabulated  in  Appendix  F.  The  transmittance 
for  this  typical  case  is  x  =0.735  and  path  radiance  Np=  7.58028  W/m2Sr.  Ship  radiance 

Ns+   at  full  range  =NS  r  +Np 

Ns+p=NsT+Np=30.14347844x0.735+7.58028=29.7357366(W/m2Sr) 
It  is  necessary  to  convert  target  radiance  at  maximum  range  to  equivalent  blackbody 
temperature.  In-band  radiance  at  full  range  can  be  converted  into  equivalent  blackbody 
temperature  by  use  of  Planck's  radiation  law  using  an  interactive  computer  program  and 
look-up  table  as  shown  in  Appendix  C.  As  an  example,  when  Ns+p=  29.7357366  W/m2Sr, 
the  effective  Ts+p  =  284.68  °K.  From  the  SEARAD  output  (Appendix  G  shows  a  typical 
case),  the  background  radiance  Nb  at  full  range  is  21.93238  W/m2Sr.  By  conversion  with 
the  computer  program  in  Appendix  C,  the  corresponding  apparent  temperature  is  found  to 
be  268.7  K.  The  apparent  temperature  difference  at  full  range  is  then  Ts+p-Tb  =15.98  °C. 
The  whole  data  set  for  the  typical  case  are  shown  in  Table  8. 

B.         MRTD  AND  MDTD  VS  RANGE  COMPUTATION 

1.  Range,  Elevation  and  Azimuth  Angle 

From  GPS  data  (Appendix  H)  we  must  compute  range  and  elevation  and  azimuth 
angles.  There  is  a  7  hours  time  difference  between  GMT  and  local  time  (PACIFIC 
Summer  Time).  For  example  the  time  of  0800  for  Monterey  will  be  1500  for  GMT. 
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Appendix  H  contains  the  GPS  data,  for  the  time  period  in  which  the  experiment  took 
place.  Figure  3  and  4  show  the  orientation  of  the  ship  and  relation  between  bearing  angle 
and  ship's  aspect.  Since  we  treat  the  ship  bow  and  stern  aspects  as  the  same,  we  only 
consider  heading  from  zero  to  1 80  degree  in  this  thesis. 

dx  =  (ship  longitude)-(sensor  longitude)  (nmi) 

dy  =  (ship  latitude)-(sensor  latitude)  (nmi) 

a  =  arctan  (dx/dy)  (degree) 


range  =^dx2  +  dy2  (nmi) 

0  =  arctan  ((sensor  altitude)  /  range)  (elevation  angle) 
To  compute  bearing  we  need  to  consider  two  different  kinds  of  situation;  those  in  which 
the  ship  is  to  the  north  or  south  of  the  sensor  position.  Range  is  computed  by  conversion 
of  GPS  location  difference  into  distance. 

a.  Condition  1 

dy>0.  This  is  illustrated  by  Figure  3. 

If  heading  >  1 80  -  a  ,  then  p  (bearing)  =  heading  -1 80  +  a 
If  heading  <  1 80  -  a  ,  then  p  =  1 80  -  a  -  heading 

b.  Condition  2 

dy<0.  See  Figure  4. 

If  heading  <  a  ,  then  p  =  a  -  heading,  if  heading  >  a,  then  p  =  heading  -  a 
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If  p  <  90  then  <j>  =  90  -  p,  if  p  >  90  then  $  =  180  -  p 

1  degree  latitude  at  37  degree  N  latitude  is  59.97  nmi 

1  degree  longitude  at  37  degree  N  latitude  is  47.995  nmi 

1  nmi  (nautical  mile  )  =  1.852  km  =1852  m 

The  computations  are  described  in  Appendix  J  section  A,"MATLAB"  code.  After 
computation  with  previously  selected  data  (see  Table  8)  we  have: 

a  =  38.6737  (degree) 

P  (bearing)  =  21.3263  (degree) 

range  =  0.7681  (nmi)  =  1.4225  (km) 

<j)  =  1.1985  (rad)  (azimuth  angle) 

G  =  0.0042  (rad)  (elevation  angle) 


2.         Projected  Area  of  the  Ship  A, 

From  previous  calculation  we  compute  projected  area. 

AT   =  Iw  sin  6  +  hw cos 6  sin  <f)+ hi  cos  6 cos  ^ 

=  41.5  •9.75»sin(0.0042)+8.8»9.75«cos(0.0042)sin(l.l985)+41.5«8.8«cos( 0.0042)  cos(l.l 985) 
=  214.4513w2 


where 


1    =  41.5  m  (ship  length) 

w  =  9.75  m  (  ship  width) 

h  =  8.8  m  (ship  height) 

4=  1.1985°  (azimuth  angle) 

6>=  0.0042°   (elevation  angle) 
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Figure  1  shows  the  orientation  of  a  ship  target  model  for  calculating  the  projected  area. 

3.  Critical  Dimension  D 

The  comparison  of  the  temperature  contrast  of  a  ship  target  with  the  standard  "bar 
target  "  for  which  the  performance  parameter  MRTD  is  derived  requires  the  definition  of 
an  equivalent  resolved  spatial  frequency  required  on  the  target  image  for  target 
recognition  (or  detection).  In  applying  the  long  standing  Johnson  criterion  [Ref.3:p.  2-7] 
of  resolved  cycles  on  the  target,  we  adopt  Moser's  [Ref.  3:p.  2-6]  convention  for  the 

"critical  dimension"  of  the  ship  =  ^Aj =  14.64  m.  Dividing  the  critical  dimension  Dcby 
the  number  of  equivalent  bars  N  on  the  target  gives  the  physical  dimension  x  to  De 
resolved  (ie.  pixel  length)  to  accomplish  the  task,  i.e.  x  =D/N.  From  the  Johnson 
Criterion  (see  Table  6),  the  number  of  bars  N  necessary  for  pure  detection,  N  =  1.  The 
number  of  bars  N  necessary  to  classify  the  ship  is  4  eye  x  2  bars/cyc  =  8  bars.  The 
number  of  bars  N  necessary  to  identify  the  ship  is  6.4cyc  x  2  bar/cyc  =  12.8  bars.  The  one 
dimensional  resolved  "footprint "  required  on  the  target  is  then: 

a.  For  detection  the  physical  dimension  x  =D/N=  1 9. 3/1 =14.64  m 

b.  For  classification  the  physical  dimension  x  =D/N=19. 3/8=1. 83  m 

c.  For  identification  the  physical  dimension  x  =D/N=  19. 3/1 2.8=1. 144  m 


Discrimination  Level 

Cycles  Across  Minimum 
Dimension 

Number  of  Bars 

Detection 

1.0 

1.0 

Classification 

4.0 

8.0 

Identification 

6.4 

12.8 

Table  6.  Johnson  Criterion  [Ref.  5:p.  54] 


27 


4.  Conversion  of  Spatial  Frequency  to  Range 

The  fundamental  spatial  frequency  fc    characterizing  the  ship  can  be  written  in 
lines  per  milliradian  as 

/<  =  20^(CyC/mRad)-  5'2 

where  R  is  the  range  in  meters  from  sensor  to  target.  The  relation  of  the  spatial  frequency 
to  range  R  and  physical  dimension  x  is  shown  in  Figure  2.  The  MRTD  and  MDTD  curves 
as  a  function  of  spatial  frequency  in  NFOV  for  a  generic  Common  Module  FLIR  have 
been  given  in  Chapter  IV  and  are  shown  plotted  in  Figure  5.  The  spatial  frequency  vs 
range  relationship  of  Equation  5.1  has  been  used  with  Equation  4.1  and  4.2  to  compute 
MRTD  and  MDTD  vs  range,  for  the  R/V  POINT  SUR  as  target,  using  the  required 
"cycles  per  critical  dimension"  Nc,  in  the  form  x  =  Dc/Nc,  where  Nc  expresses  the  Moser 
definition  of  "critical  dimension".  The  detailed  calculations  are  in  Appendix  J, 
"MATLAB  code".  Figure  6  shows  MRTD  and  MDTD  for  detection  as  function  of  range, 
for  the  R/V  POINT  SUR  as  target.  The  target  available  contrast  temperature  difference 
against  background  ATapp  is  also  shown  on  the  same  range  scale.  The  intersection  points 
of  AT  with  the  MRTD  and  MDTD  curves  yield  the  maximum  ranges  for  detection 
according  to  the  two  criteria.  Figure  7  shows  the  corresponding  MRTD  for  classification 
of  the  R/V  POINT  SUR,  according  to  the  Johnson  criterion  of  required  resolution  line 
pairs  on  the  critical  dimension.  MDTD  cannot  be  used  for  classification.  Figure  8  shows 
the  equivalent  curve  for  MRTD  for  identification  of  the  POINT  SUR.  The  value  of 
ATapp(R)  in  these  three  figures  is  the  simple  estimate  of  AT0  multiplied  by  the 
transmittance  (see  Section  VI.  A). 
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VI.  DATA  ANALYSIS  AND  CONCLUSIONS 

The  computations  of  target  temperature  difference  AT,  MDTD  and  MRTD  for  the 
experiment  are  described  in  the  previous  chapter.  Meteorological  data  (including 
radiosonde  measurement)  for  SEARAD  code  inputs  were  taken  from  the  meteorological 
data  that  were  recorded  during  the  experiment,  and  are  shown  in  Appendix  H  and 
Appendix  I.  The  estimated  MRTD  curve  and  MDTD  curve  were  derived  from 
Shumaker's  example  8-13[Ref.  3:p.  8-59]  for  a  generic  FLIR  system  in  NFOV.  The  NET 
was  chosen  as  0.125  K,  between  0.1  K  to  0.25  K,  due  to  degradation  of  the  FLIR  sensor 
after  being  in  service  for  several  years.  This  sensor  is  a  generic  common  module  FLIR 
model,  with  the  parameters  of  Table  2.  The  condition  for  detection  at  maximum  range  is 
that  the  apparent  target/background  temperature  difference  as  modified  by  atmospheric 
path  attenuation  and  path  radiance  will  match  the  required  MDTD  or  MRTD  at  the  same 
range.  We  have  already  expressed  the  MRTD  and  MDTD  parameters  as  functions  of 
range.  We  are  now  required  to  calculate  the  apparent  temperature  difference  between 
target  and  background  as  a  function  of  range  from  the  target.  Two  methods  of  differing 
complexity  are  available  for  this  computation;  the  greatly  simplified  approach  of 
assuming  that  the  temperature  difference  AT  attenuates  through  the  atmosphere  in  the 
same  way  as  the  irradiance,  and  the  more  complicated  but  more  rigorous  approach  of 
converting  the  AT  into  a  radiance  difference,  computing  the  transmitted  radiance 
difference  due  to  attenuation  and  path  radiance,  and  reconverting  AL  back  to  ATapp 
through  inversion  of  Planck's  Law.  Sample  computations  of  both  types  have  been  made. 
These  are  compared  in  the  following  paragraphs,  using  the  sample  atmospheric  data  set  of 
Table  8. 
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A.  APPARENT  AT  BY  ATTENUATION  OF  SOURCE  AT 

SEARAD  with  the  "typical  case"  (see  Table  8)  scenario  is  found  to  give  a 
transmittance  of  0.735.  Transmittance,  averaged  over  a  bandwidth  which  does  not  contain 
absorption  edges,  typically  is  written  as  an  exponential  function  of  range  following 
Beer's  Law.  From  transmittance  we  calculate  attenuation  coefficient  averaged  over  the 
assumed  bandpass  of  8  to  12  urn,  of 

,,r  ,  ^  -lnr 

^        R 

0.735  =  e"MR  =>  ln(0.735)  =  -n(l. 444km)  =>  \i  =  0.2132  km"1 
For  transmittance  x  =  0.735  over  a  range  of  1.4444  km  for  horizontal  range  1 .4225 
km,  the  attenuation  coefficient  n  =  0.2132  km"1.  The  apparent  temperature  difference  as  a 
function  of  range  is  then  computed  with  the  form  ATx  t(R)  =  AT*  e'**  using 
MATLAB.  For  a  given  target,  the  functions  MRTD,  MDTD  and  ATapp  can  be  plotted  on 
the  same  scale  as  functions  of  range.  The  intersection  of  the  ATapp  plot  with  the 

appropriate  MRTD  or  MDTD  curve  occurs  at  the  maximum  range  for  the  detection, 
classification  or  identification,  according  to  the  defined  criterion.  Figure  6  shows  The 
MRTD-derived  detection  range  as  21  km  and  the  MDTD-derived  value  as  19  km.  The 
classification  range  (from  MRTD)  is  seen  in  Figure  7  to  be  8.8  km.  Figure  8  shows 
identification  range  as  6.2  km. 

B.  RADIANCE  COMPUTATION  METHOD 

The  second  method  of  computation  of  effective  AT  is  to  convert  the  target  and 
background  temperatures  to  radiances,  and  compute  the  attenuated  AT  and  appropriate 
path  radiance  for  selected  range,  and  reconvert  the  result  into  ATapp.  This  method  was 

applied  for  a  set  of  ranges  appropriate  to  the  measurement  scenario  and  data  set.  Tables 
9.1  to  9.5  show  a  tabulation  of  the  computed  parameters  for  the  sensor  elevation  of  6 
meters,  the  AGEMA  imager  location  at  Moss  Landing.  The  locations  were  selected  from 
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those  for  which  observed  images  were  recorded  in  EOMET95,  so  that  the  results  may 
later  be  applied  to  comparisons  of  polarization.  All  the  computed  components  of  AT,  the 
difference  between  the  apparent  values  of  ship  average  temperature  corrected  for 
attenuation  and  background  temperature,  are  shown  in  this  table.  The  independent 
variable  is  the  zenith  angle,  used  to  define  the  path  geometry  for  this  sensor  location.  The 
zero  range  ship  radiance  is  also  input.  The  SEARAD  program  computes  the 
corresponding  range,  transmittance  and  path  radiance,  and  apparent  ship  radiance,  and  the 
background  radiance.  The  radiance  contrast  is  then  converted  back  into  temperature 
difference.  Table  9.5  shows  an  anomaly  in  range  and  a  path  that  does  not  intersect  the  sea 
surface  beyond  13.3  km;  this  range  defines  the  effective  horizon  for  the  selected  sensor 
elevation  of  6  meters.  It  follows  that  we  cannot  compute  predicted  maximum  range  for 
detection  or  recognition  for  ranges  greater  than  this.  To  allow  range  predictions  to  be 
extended  to  greater  ranges  a  further  set  of  data  was  computed,  for  sensor  elevation  raised 
to  250  meters.  This  allows  us  to  see  range  limitation  by  atmospheric  effects  rather  than  by 
geometry.  The  results  of  these  computations  are  shown  as  Tables  10.1  to  10.6.  From 
these  tables  we  can  see  that  AT  attenuates  almost  to  zero  at  long  range.  AT  is  shown 
plotted  against  range  in  Figure  9  and  Figure  1 0.  Figure  9  represents  the  original  situation 
with  sensor  at  6  meters  elevation;  Figure  10  represents  the  same  target  and  atmospheric 
scenario,  but  with  the  sensor  raised  to  an  elevation  of  250  m.  This  revised  curve  allows 
the  intersection  of  the  AT(R)  curve  with  MRTD(R)  to  be  observed,  and  the  Maximum 
Detection  Range  extracted.  This  curve  is  shown  as  Figure  11.  Figure  14  shows  AT  vs 
MRTD  for  classification  and  identification.  The  apparent  MRD  is  44  km. 

C.         SEARAD 

The  SEARAD  code  input  file  requires  the  meteorological  data  inputs  wind  speed 
(Ws),  wind  direction  (W[)),  visibility  (VIS),  air  mass  character  (AM)  and  the  sea 
temperature.  A  sample  input  file  for  SEARAD  is  shown  in  Appendix  F.  The  contents  of 
the  three  header  lines  in  Appendix  F  are  defined  by  the  "cards"  listed  in  Appendix  E. 
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Since  SEARAD  reads  in  the  data  it  must  to  be  in  the  correct  column  and  row.  The  sample 
input  file  shows  "card  2C1"  which  is  for  radiosonde  data  to  be  inserted.  The  last  digit  of 
the  layer  altitude  must  to  be  in  column  1 0  and  the  last  digit  of  pressure  must  be  in  column 
20.  The  last  digit  for  the  temperature  of  the  layer  boundary  must  be  in  column  30  and  the 
last  digit  of  the  relative  humidity  must  be  in  column  40.  The  designated  specific  units 
must  be  in  columns  61,  62  and  63.  From  the  sample  input  file  we  can  see  they  should  be 
lined  up  column  by  column  and  row  by  row.  The  sample  output  from  SEARAD 
(atmospheric  transmittance  rand  path  radiance  Np)  are  shown  in  Appendix  G.  To  read  in 
a  new  model  atmosphere  from  radiosonde  data,  the  sensor's  altitude  can  not  be  greater 
than  the  upper  layer  boundary.  SEARAD  will  use  the  altitude  of  the  highest  layer 
boundary  rather  than  the  sensor  elevation.  To  use  a  high  altitude  sensor,  it  is  necessary  to 
use  a  higher  layer  boundary  for  the  calculation  of  atmospheric  transmittance  r  and  path 
radiance  Np  . 

D.         DISCUSSION  OF  THE  DATA  AND  CONCLUSIONS 

Based  on  a  sample  scenario  selected  from  the  data  base  of  the  EOMET95 
measurement  series  in  Monterey  Bay  off  Moss  Landing,  ship  skin  temperature  samples 
measured  with  a  set  of  transistors  were  used  with  historical  thermal  distributions  of  the 
R/V  POINT  SUR,  as  measured  radiometrically,  to  estimate  an  average  ship  temperature 
for  the  selected  measurement  time.  From  this  together  with  the  sea  background 
temperature  computed  by  SEARAD,  a  "zero  range"  contrast  temperature  difference  AT 
was  computed  as  289.987  -  286.600  =  3.3871.  This  AT  has  been  applied  in  two  ways  to 
estimate  the  Maximum  Detection  Range  for  the  R/V  POINT  SUR.  The  first  and  simpler 
method  is  to  assume  that  the  contrast  temperature  difference,  AT,  is  attenuated  by  the 
atmosphere  in  the  same  fashion  as  the  radiance,  ie  ATapp  =  AT0  e "•*  .  The  second  ,  more 
justifiable,  method  is  to  convert  the  temperature  difference  into  in-band  background 
radiance  and  target  radiance,  compute  the  corrections  for  attenuation  and  path  radiance 
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using  SEARAD,  and  reconvert  to  an  apparent  AT  as  a  function  of  range.  For  each  method 
ATapp  is  then  compared  with  the  MRTD  or  MDTD  of  the  assumed  sensor,  as  appropriate, 
also  expressed  as  a  function  of  range.  The  point  of  intersection  then  gives  the  maximum 
detection  or  recognition  range.  A  weakness  of  this  method  is  that  the  ship  radiance  is 
based  solely  on  the  thermal  emission  from  the  ship  while  the  background  calculation 
takes  into  account  the  reflected  radiance  also.  The  following  observations  are  drawn  from 
the  results. 

We  see  in  Table  9  (page  53)  and  Figure  9  that  the  apparent  temperature  difference 
at  short  range  initially  increases.  At  0.345  km  ATapp  is  21.18,  much  higher  than  the  zero 
range  value  of  3.387.  After  reaching  a  peak  value,  the  apparent  temperature  difference 
then  decreases  in  an  approximately  exponential  fashion  with  range. 

For  sensor  altitude  of  6  meters,  beyond  elevation  angle  0.044,  corresponding  to  a 
range  of  13.262  km  the  calculated  refracted  path  no  longer  intersects  the  earth  surface. 
The  path  radiance  cannot  then  be  calculated,  and  the  MDR  cannot  be  calculated.  This 
geometry  is  not  appropriate  for  estimation  of  maximum  detection  range.  The 
recalculation  with  assumed  higher  elevation  for  the  sensor,  required  a  redescription  of  the 
atmospheric  layer  structure,  and  allows  calculation  of  the  path  radiance  and  surface 
radiance.  We  see  in  Table  10,  with  a  sensor  altitude  of  250  m,  ATapp  has  a  value  of  0.007 
degrees  at  range  33.415  km,  with  zenith  angle  0.550  degrees.  The  variation  of  AT  with 
range  (Figure  10)  shows  small  differences  from  Figure  9.  This  must  be  ascribed  to  the 
reassignment  of  layer  structure  in  SEARAD.  The  initial  increase  of  ATapp  at  short  range 
persists,  and  the  curve  is  basically  similar  to  Figure  9,  except  that  it  is  extended  to  much 
longer  ranges,  so  that  estimation  of  MDR  is  possible.  A  comparison  of  the  ATapp  (R)  as 
computed  by  the  two  methods  is  shown  as  Table  7. 

A  direct  comparison  of  these  range  predictions  indicates  that  the  simple  Beer's 
Law  variation  of  ATapp  (R)  leads  to  an  underestimation  of  the  maximum  range,  relative  to 
the  SEARAD  method  including  computed  transmittance  and  path  radiance.  The  values 
shown  in  Table  7  are  comparable  to  values  obtained  in  other  range  models.  However  a 
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factor  of  approximately  2  is  apparent  between  the  sets  of  results  for  detection  with 
smaller  differences  for  classification  or  identification.  This  point  requires  further 
discussion.  The  selection  of  the  MDTD  or  MRTD  criterion  for  detection  appears  to  be  of 
much  lesser  importance  than  the  other  differences  between  the  model  approaches. 


Detection 

Classification 

Identification 

Range 

Range 

Range 

I.  Extinction 

MRTD:  21  km 

MRTD:  8.8  km 

MRTD:6.2  km 

only 

MDTD:  19  km 

II.  SEARAD 

MRTD:  47  km 

MRTD:  13  km 

MRTD:8  km 

Correction 

MDTD:  44  km 

Table.  7  Comparison  of  Two  Methods 

Some  insight  into  the  discrepancy  may  be  gained  by  comparison  of  the 
transmittances  calculated  with  the  Beer's  Law  approach  and  the  SEARAD  code.  Figure 
1 5  shows  this  direct  comparison.  It  can  be  seen  that  the  Beer's  Law  method  considerably 
underestimates  the  transmittance  compared  with  direct  computation.  This  is  probably  due 
to  the  band  averaging  over  the  8  to  12  band,  for  which  the  apparent  extinction  coefficient 
decreases  initially  approaching  a  constant  value  only  at  longer  ranges.  The  coefficient 
used  in  this  analysis  was  deduced  from  extinction  at  1 .4  km,  a  measurement  point  in  the 
EOMET95  data  set.  The  differences  observed  in  Figure  1 5  will  decrease  the  discrepancy 
in  predicted  range  but  only  in  a  minor  way,  as  can  be  seen  from  the  plot  of  apparent 
temperature  difference  shown  as  Figure  1 6.  The  SEARAD  method  shows  negative  values 
at  very  short  range  rising  to  values  of  up  to  20°C,  giving  an  increased  effective  AT0, 
before  decreasing  in  an  approximately  exponential  fashion.  The  greatest  factor  in  the 
range  difference  is  this  high  AT0.  The  SEARAD  and  Beer's  Law  methods  use  the  same 
input  radiometric  sea  surface  temperature  and  thermometric  target  temperature.  The 
SEARAD  code  then  computes  the  sea  surface  radiance  taking  into  account  the  sea  surface 
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roughness  and  emissivity,  and  the  reflection  of  sky  radiance.  The  corresponding  radiance 
temperature  is  then  computed  and  used  in  the  initial  temperature  AT0.  This  procedure 
leads  to  the  initial  negative  temperature  difference  seen  in  Tables  9  and  1 0,  and  in  Figure 
16.  No  equivalent  radiance  temperature  model  is  available  for  the  ship  target. 

The  increase  in  ATapp  at  short  ranges  is  an  unexpected  feature.  Examination  of 
Tables  9  and  10  shows  that  the  ship  radiance  corrected  for  path  radiance  and 
transmittance  decreases  slowly  and  monotonically  with  range.  The  background  radiance, 
however,  decreases  from  a  large  initial  value  to  a  minimum  occurring  at  92°  zenith  angle 
for  the  6  meter  sensor  and  95°  for  the  250  meter  sensor,  and  then  rises  slowly  at  longer 
range.  This  behavior  is  not  completely  understood.  It  is  non-intuitive  that  the  sea 
radiance  and  the  target  radiance  should  behave  differently  at  long  range.  The  initial 
behavior  at  short  range  may  be  related  to  the  variation  of  the  emissivity  and  sky  radiance. 
This  behavior  deserves  further  study. 

It  may  be  concluded  that  the  use  of  SEARAD  for  range  prediction  constitutes  a 
considerable  improvement  over  the  Beer's  Law  method.  However  a  valid  target  signature 
radiance  model  is  needed  to  give  confidence  in  the  predicted  range.  The  SHIPSIG  and 
THERMAL  CONTRAST  MODEL  have  been  used  in  some  previous  prediction  models, 
but  are  very  specific  to  target  and  conditions,  and  so  not  convenient  for  general 
application.  Addition  of  a  simple  facet  ship  signature  representation  to  the  ship 
orientation  model  is  recommended. 
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LOS. 


Figure  1 .  Orientation  of  Ship  Target  Model 
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2000^ 
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N     =  Number  of  Bars  in  the  Bar  Chart 
R     =  Slant  Range 
TAS  =  Target  Angular  Subtense 
&C    =  Critical  Dimension 
X     =  Physical  Dimension  to  be  Resolved 
*c    =  Spatial  Frequency  of  Bar 


Figure  2.  Relation  Between  Spatial  Frequency,  Range  R  and  Physical  Dimension  % 


37 
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Condition  1 


dx  =  (ship's  longitude)  -  (sensor's  longitude) 
dy  =  (ship's  latitude)  -  (sensor's  latitude) 
a  =  arctan  (dx/dy) 
4>  =  (90-P) 

0  =  arctan  (6  /  Range) 

1  degree  latitude  at  37  degree  N  latitude  is  59.97  nmi 

1  degree  longitude  at  37  degree  N  latitude  is  47.995  nmi 
1  nmi=  1.852  km  =1852  m 


Figure  3.  Relation  Between  GPS  and  Ship  Aspect  Condition  1 
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Sensor's 
position 
121  80  W 
36  80N 


Condition  2 

dx  =  (ship's  longitude)  -  (sensor's  longitude) 
dy  =  (ship's  latitude)  -  (sensor's  latitude) 
a  =  arctan  (dx/dy) 
4>  =  (90-p) 

0  =  arctan  (6  /  Range) 

1  degree  latitude  at  37  degree  N  latitude  is  59.97  nmi 

1  degree  longitude  at  37  degree  N  latitude  is  47.995  nmi 
1  nmi -1.852  km  =1852  m 


Figure  4.  Relation  Between  GPS  and  Ship  Aspect  Condition  2 
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Figure  5.  MRTD  and  MDTD  as  a  Function  of  Spatial  Frequency 


40 


3.5 

3 

2.5 

9>      2 

13 

CO 

&_ 

0) 
Q. 

E 

h-   1.5 

1 

0.5 

MRT,  MDT  &  dT  vs  Range  for  Detection 

i 

i 

A HT 

U  r 

MRT 

, 

I 

\ 

. 

MDT 

) 

I             : 

i            : 

i 

i            i 

5            ( 

•\            - 

7                 1 

3                 9 

Range (m) 


x1(T 


Figure  6.  MRTD,  MDTD  and  AT  as  a  Function  of  Range  for  Detection 


of  R/V  POINT  SUR,  AT  by  Beer's  Law 
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MRT  &  dT  vs  Range  for  Classification 
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Figure  7.  MRTD  and  AT  as  a  Function  of  Range  for  Classification; 


AT  by  Attenuation  of  ATC 
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MRT  &  dT  vs  Range  for  Identification 
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Figure  8.  MRTD  and  AT  as  a  Function  of  Range  for  Identification: 


AT  by  Attenuation  of  Source  Method 
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Temperature  Difference  vs  Range 
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Figure  9.  AT  vs  Range  for  Typical  case,  by  Radiance  Compensation  Method, 
(SEARAD)  Sensor  at  6  meters  Elevation 
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Figure  10.  AT  vs  Range  for  MDR  Estimate; 


Sensor  at  250  meters  Elevation 
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Figure  1 1 .  AT  vs  MRTD  &  MDTD  for  MDR  Estimate 
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Figure  12.  Locations  of  Thermistors 
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TEMPERATURE  DISTRIBUTION  TABLE,  R/V  POINT  SUR    PORT  BEAM 
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Figure  13.  Distribution  of  Thermistors,  and  Correlation  with  Pixels 
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<JT  vs  MRT  for  Classification  and  Identification 
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Figure  14.  AT  vs  MRTD  for  Classification  and  Identification 
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Extinction;  SEARAD  vs  Beer's  Law 
Beers  Law  from  Tau  at  1 .444  km 
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Figure  15.  Comparison  of  Extinction;  SEARAD  vs  Beer's  Law 
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Ship  Background  Temperature  Difference  vs  Range 

Beers  Law  vs  SEARAD 
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Figure  16.  AT  for  Beer's  Law  vs  SEARAD 
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Date 

95/05/16 

Time  (GMT) 

18:37:53 

Time  (local) 

11:37:53 

Ship  Heading 

120 

Wind  Speed  (knot) 

3.69 

Wind  Speed  (m/s) 

1.8966 

Wind  Speed  (knot)  24-H  Average 

3.1201 

Wind  Speed  (m/s)  24-H  Average 

1.6037 

Wind  Direction 

252 

Air  Temperature  (°C) 

12.30 

Relative  Humility  (%) 

93.34 

Pressure  (mb) 

1012 

Latitude 

36.81  N 

Longitude 

121.81  W 

Ship  Temperature  (°C) 

16.987144 

Ship  Temperature  (°K) 

289.987144 

Sea  Temperature  (°C) 

13.60 

Sea  Temperature  (°K) 

286.60 

Ship  Radiance  Ns  (zero  range) 

30.14347844 

Transmittance  t 

0.735 

Path  Radiance  Np 

7.58028 

Nsxx+Np=Ns+p 

29.7357366 

N     =>T 

284.68 

Background  Radiance  Nb 

21.93238 

Nb  =>  Tb 

268.7 

*■  s+p  —  *  b  —  ^  *■  app 

15.98 

Observe  Range(km)  horizontal 

1.4225 

Observe  Range  (km) 

1.444 

AT 

3.387144 

Bearing  (degree) 

21.326 

Azimuth  angle  (degree) 

68.674 

Elevation  (degree) 

0.24066 

Ship  Project  Area  (m2) 

214.4513 

Table  8.  Data  File  for  Typical  Case 
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Observed 
Range  (km) 

Zero 
Range 

0.006 

0.021 

0.039 

0.069 

Zenith 
Angle 

0.0 

158.500 

107.000 

99.000 

95.000 

Emissivity 

8 

0.9876 

0.8518 

0.6937 

0.5887 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

0.0 

0.9924 

0.9846 

0.9773 

0.9670 

Path  Radiance 

Np 

^c^c^^^t 

0.19415 

0.39946 

0.59549 

0.87565 

Nsxx+Np 
=  NS+P 

30.10853 

30.07872 

30.05470 

30.02438 

Ns+P  convert  to 
Temperature 

285.37 

285.32 

285.27 

285.22 

Background 
Radiance 

30.63620 

27.41545 

23.96328 

22.04542 

Nb  convert  to 
Temperature  Tb 

286.35 

280.25 

273.25 

269.05 

Ts+P-Tb  =  AT 

3.387144 

-0.98 

5.07 

12.02 

16.17 

Table  9.1  Data  File  for  Typical  Case  Radiance  Computation 
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Observed 
Range  (km) 

0.086 

0.115 

0.138 

0.172 

0.230 

Zenith 
Angle 

94.000 

93.000 

92.500 

92.000 

91.500 

Emissivity 
s 

0.5613 

0.5351 

0.5227 

0.5099 

0.4954 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

X 

0.9618 

0.9538 

0.9478 

0.9393 

0.9262 

Path  Radiance 

Np 

1.01743 

1.23811 

1.40434 

1.64105 

2.01083 

Nsxx+Np 
=  NS+P 

30.00941 

29.9889 

29.97432 

29.95481 

29.92926 

Ns+P  convert  to 
Temperature 

285.19 

285.15 

285.13 

285.09 

285.04 

Background 
Radiance 

21.46527 

21.19206 

21.14012 

20.68977 

20.70253 

Nb  convert  to 
Temperature  Tb 

267.75 

267.05 

266.95 

265.85 

265.95 

Ts+P-Tb  =  AT 

17.44 

18.1 

18.18 

19.24 

19.09 

Table  9.2  Data  File  for  Typical  Case  Radiance  Computation 
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Observed 
Range  (km) 

0.345 

1.740 

3.616 

4.071 

4.669 

Zenith 

Angle 

90.999 

90.200 

90.100 

90.090 

90.080 

Emissivity 

8 

0.4820 

0.4630 

0.4617 

0.4616 

0.4615 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

X 

0.9023 

0.6995 

0.5210 

0.4868 

0.4457 

Path  Radiance 

Np 

2.69034 

8.63780 

13.99925 

15.03495 

16.27719 

Ship  Radiance 

Ns+P 

29.88879 

29.72315 

29.70399 

29.70879 

29.71213 

Ns+P  convert  to 
Temperature  Ts+P 

284.96 

284.66 

284.62 

284.63 

284.64 

Background  Radiance 

Nb 

20.78745 

22.26335 

24.05040 

24.40729 

24.83836 

Nb  convert  to 
Temperature  Tb 

266.09 

269.40 

273.4 

274.1 

275.03 

Ts+p-Tb  =  AT 

18.87 

15.57 

11.22 

10.53 

9.61 

Table  9.3  Data  File  for  Typical  Case  Radiance  Computation 
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Observed 
Range  (km) 

5.501 

6.772 

9.167 

9.557 

10.008 

Zenith 
Angle 

90.070 

90.060 

90.050 

90.049 

90.048 

Emissivity 
s 

0.4614 

0.4613 

0.4611 

0.4611 

0.4611 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

T 

0.3951 

0.3298 

0.2366 

0.2243 

0.2110 

Path  Radiance 
Np 

17.80972 

19.78584 

22.60054 

22.97070 

23.37262 

Ship  Radiance 

Ns+P 

29.71940 

29.72715 

29.73248 

29.73188 

29.73289 

Ns+P  convert  to 
Temperature  Ts+P 

284.65 

284.67 

284.68 

284.67 

284.68 

Background  Radiance 

Nb 

25.37396 

26.06976 

27.06865 

27.20058 

27.34404 

Nb  convert  to 
Temperature  Tb 

276.15 

277.57 

279.56 

279.83 

280.11 

Ts+P-Tb  =  AT 

8.5 

7.1 

5.12 

4.84 

4.57 

Table  9.4  Data  File  for  Typical  Case  Radiance  Computation 
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Observed 
Range  (km) 

10.534 

11.174 

12.009 

13.262 

97.261 

Zenith 
Angle 

90.047 

90.046 

90.045 

90.044 

90.043 

Emissivity 
s 

0.4611 

0.4610 

0.4610 

0.4610 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

X 

0.1965 

0.1802 

0.1611 

0.1363 

0.0001 

Path  Radiance 

Np 

23.80918 

24.29645 

24.86820 

25.60797 

Ship  Radiance 

Ns+P 

29.73237 

29.72830 

29.72431 

29.71652 

Ns+P  convert  to 
Temperature  Ts+P 

284.66 

284.67 

284.66 

284.65 

Background  Radiance 

Nb 

27.50002 

27.67432 

27.87898 

28.14412 

Nb  convert  to 
Temperature  Tb 

280.43 

280.75 

281.15 

281.67 

Ts+P-  Tb  =  AT 

4.23 

3.92 

3.51 

2.98 

Table  9.5  Data  File  for  Typical  Case  Radiance  Computation 
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Observe 
Range  (km) 

0.259 

0.289 

0.500 

1.441 

2.875 

Zenith 
Angle 

195.000 

150.000 

120.000 

100.000 

95.000 

Emissivity 
s 

0.9877 

0.9869 

0.9539 

0.7180 

0.5882 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

0.9311 

0.9258 

0.8920 

0.7764 

0.6457 

Path  Radiance 
Np 

1.86920 

2.018832 

2.97817 

6.34750 

10.23864 

NsxT+Np 
=  N 

29.93578 

29.92565 

29.86614 

29.75089 

29.70227 

Ns+P  convert  to 
Temperature  Ts+P 

285.06 

285.04 

284.93 

284.71 

284.62 

Background  Radiance 

Nb 

30.59705 

30.57243 

29.77694 

25.48528 

24.50860 

Nb  convert  to 
Temperature  Tb 

286.27 

286.23 

284.76 

276.35 

274.35 

Ts+P-Tb  =  AT 

-1.21 

-1.19 

0.17 

8.36 

10.27 

Table  10.1  Data  File  for  MDR  Computation 
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Observed  Range 
(km) 

3.596 

4.805 

7.260 

9.783 

11.382 

Zenith 
Angle 

94.000 

93.000 

92.000 

91.500 

91.300 

Emissivity 
s 

0.5607 

0.5344 

0.5084 

0.4935 

0.4882 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

T 

0.5919 

0.5140 

0.3909 

0.2985 

0.2526 

Path  Radiance 

Np 

11.85932 

14.21330 

17.95142 

20.76567 

22.16049 

Nsxx+Np 

:  ^S+P 

29.70123 

29.70704 

29.73450 

29.76349 

29.77473 

Ns+P  convert  to 
Temperature  Ts+P 

284.62 

284.63 

284.68 

284.73 

284.76 

Background  Radiance 

Nb 

24.55046 

25.01152 

25.88875 

26.76083 

27.20808 

Nb  convert  to 
Temperature  Tb 

274.9 

275.49 

277.7 

278.95 

279.85 

Ts+P-Tb  =  AT 

9.72 

9.14 

6.98 

5.78 

4.91 

Table  10.2  Data  File  for  MDR  Computation 
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Observed  Range 
(km) 

13.776 

15.334 

17.098 

19.663 

21.299 

Zenith 

Angle 

91.090 

90.990 

90.900 

90.800 

90.750 

Emissivity 
c 

0.4818 

0.4785 

0.4756 

0.4725 

0.4709 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

0.1977 

0.1690 

0.1418 

0.1101 

0.0939 

Path  Radiance 

Np 

23.82643 

24.69582 

25.51815 

26.46917 

26.95469 

Nsxx+Np 
=  NS+P 

29.78579 

29.79006 

29.79249 

29.78796 

29.78516 

Ns+P  convert  to 
Temperature  Ts+P 

284.77 

284.78 

284.79 

284.78 

284.78 

Background  Radiance 

Nb 

27.73904 

28.01768 

28.28500 

28.59857 

28.76048 

Nb  convert  to 
Temperature  Tb 

280.95 

281.4 

281.9 

282.5 

282.8 

Ts+P-Tb  =  AT 

3.82 

3.38 

2.89 

2.28 

1.98 

Table  10.3  Data  File  for  MDR  Computation 
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Observed  Range 
(km) 

23.278 

25.743 

28.948 

31.421 

33.415 

Zenith 
Angle 

90.700 

90.650 

90.600 

90.570 

90.550 

Emissivity 

8 

0.4694 

0.4679 

0.4664 

0.4655 

0.4649 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

T 

0.0776 

0.0613 

0.0452 

0.0358 

0.0297 

Path  Radiance 
Np 

27:44248 

27.92743 

28.40245 

28.67854 

28.85722 

Nsxx+Np 
:  Ns+p 

29.78161 

29.77522 

29.76493 

29.75767 

29.75248 

Ns+P  convert  to 
Temperature  Ts+P 

284.77 

284.76 

284.74 

284.72 

284.71 

Background  Radiance 

Nb 

28.92441 

29.08870 

29.25099 

29.34601 

29.40779 

Nb  convert  to 
Temperature  Tb 

283.16 

283.47 

283.77 

283.95 

284.07 

Ts+P-Tb  =  AT 

1.61 

1.29 

0.97 

0.77 

0.64 

Table  10.4  Data  File  for  MDR  Computation 
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Observed  Range 
(km) 

35.809 

37.210 

38.962 

39.310 

42.769 

Zenith 
Angle 

90.530 

90.520 

90.509 

90.507 

90.490 

Emissivity 

8 

0.4643 

0.4640 

0.4637 

0.4636 

0.4630 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

0.0237 

0.0208 

0.0177 

0.0171 

0.0124 

Path  Radiance 

Np 

29.03059 

29.11495 

29.20565 

29.22193 

29.35741 

Nsxx+Np 
:  Ns+p 

29.74499 

29.74193 

29.73918 

29.73738 

29.73118 

Ns+P  convert  to 
Temperature  Ts+P 

284.70 

284.69 

284.69 

284.68 

284.67 

Background  Radiance 

Nb 

29.46800 

29.49740 

29.52908 

29.53478 

29.58232 

Nb  convert  to 
Temperature  Tb 

284.18 

284.24 

284.30 

284.31 

284.40 

Ts+P-  Tb  =  AT 

0.52 

0.45 

0.39 

0.37 

0.27 

Table  10.5  Data  File  for  MDR  Computation 
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Observed  Range 
(km) 

45.409 

48.877 

54.146 

58.641 

71.559 

Zenith 
Angle 

90.480 

90.470 

90.460 

90.455 

90.450 

Emissivity 
s 

0.4626 

0.4623 

0.4618 

0.4615 

0.4610 

Ship  Radiance 

Ns 

30.14347 

30.14347 

30.14347 

30.14347 

30.14347 

Transmittance 

X 

0.0097 

0.0070 

0.0043 

0.0028 

0.0008 

Path  Radiance 

Np 

29.43472 

29.51111 

29.58875 

29.63071 

29.68708 

Nsxx+Np 
=  NS+P 

29.72711 

29.72211 

29.71836 

29.71511 

29.71111 

Ns+P  convert  to 
Temperature  Ts+P 

284.67 

284.66 

284.65 

284.64 

284.63 

Background  Radiance 

Nb 

29.60958 

29.63663 

29.66433 

29.67943 

29.69996 

Nb  convert  to 
Temperature  Tb 

284.44 

284.50 

284.55 

284.58 

284.62 

Ts+P-  Tb  =  AT 

0.23 

0.16 

0.10 

0.06 

0.01 

Table  10.6  Data  File  for  MDR  Computation 
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Channel 

Location 

Temperature 

Pixel 

Tempx  Pixel 

1 

Reference 

25.001 

**** 

***** 

2 

Starboard  bow 

16.423 

10 

164.23 

3 

Aft  pilot  house 

17.864 

**** 

**** 

4 

Port  aft  of  stack 

17.291 

2 

34.582 

5 

Aft  on  stack 

15.361 

**** 

**** 

6 

Aft  port 

15.830 

26 

411.58 

7 

Aft  port  water 

14.368 

**** 

**** 

8 

Port  pilot  house 

15.909 

11 

174.999 

9 

Reference 

25.001 

**** 

**** 

10 

Bow  pilot  house 

18.333 

**** 

**** 

11 

Port  stack 

24.164 

5 

120.82 

12 

Aft  starboard 

15.296 

22 

336.512 

13 

Port  below  stack 

16.151 

5 

80.755 

14 

Port  bow 

15.346 

10 

153.46 

15 

Bow  air 

22.393 

**** 

**** 

16 

Starboard  stack 

28.829 

4 

115.316 

** 

Starboard  above 
water-line 

17.923 

23 

412.229 

Total 

********* 

***** 

118 

2004.483 

Table  1 1 .  Distribution  and  Location  of  Thermistors 
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APPENDIX  A.  THERMISTOR  SAMPLE  OUTPUT  FILE 

May  16.1 

Mon,Day,Year  =  5  16  1995 

Hour,Min,Sec  -  8  27  36 

volts,istart,istop,nsam,iref  1.49939  1  16  100  1 

1326456  25.001  12.856  13.403  11.963  14.615  11.677  12.359  11.634  25.001  14.242 
22.956  12.010  11.617  11.431  15.410  25.874 

1326476  25.001  12.864  13.397  11.959  14.646  11.673  12.371  11.633  25.001  14.200 
22.937  12.001  11.607  11.426  15.229  25.535 

1326496  25.001  12.858  13.394  11.954  14.699  11.662  12.415  11.647  25.001  14.162 
22.924  11.993  11.598  11.420  15.138  25.186 

1326516  25.001  12.861  13.404  11.950  14.752  11.662  12.437  11.625  25.001  14.145 
22.923  11.989  11.593  11.416  15.062  24.958 

1337856  25.001  16.392  17.862  17.334  15.415  15.885  14.338  15.873  25.001  18.294 
24.326  15.281  16.169  15.352  22.299  28.721 

1337876  25.001  16.423  17.864  17.291  15.361  15.830  14.368  15.909  25.001  18.333 
24.164  15.296  16.151  15.346  22.393  28.829 

1337896  25.001  16.458  17.844  17.241  15.295  15.771  14.383  15.946  25.001  18.366 
24.093  15.304  16.117  15.305  22.319  28.933 

1337916  25.001  16.495  17.819  17.200  15.242  15.716  14.404  15.972  25.001  18.395 
24.064  15.318  16.083  15.264  22.276  29.048 

1337936  25.001  16.528  17.796  17.160  15.188  15.677  14.354  15.997  25.001  18.421 
24.021  15.341  16.061  15.248  22.188  29.162 

1337956  25.001  16.559  17.799  17.118  15.150  15.638  14.274  15.983  25.001  18.451 
23.944  15.370  16.044  15.256  22.322  29.302 
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APPENDIX  B.  PROPORTIONAL  RADIATION  TABLE 

The  "Proportional  Radiation",  q=f(A,  ,  T)  represents  the  fraction  of  the  radiant 
exitance  which  is  emitted  by  a  blackbody  at  temperature  T  at  all  wavelengths  up  to  the 
selected  values  of  X.  It  is  obtained  from  the  integral  over  the  Planck  Law  for  wavelengths 
up  to  the  selected  values  divided  by  the  integral  over  all  wavelengths,  i.e. 


A. 

\M{XT)dX 

_o 

00 

\M(XT)dX 


XTIcm  K 

q 

XTIcm  K 

q 

AT/cm  K 

q 

0.050 

1.3652E-9 

0.140 

7.9035E-3 

0.460 

5.8075E-1 

0.052 

3.6788E-9 

0.150 

1.3023E-2 

0.480 

6.0880E-1 

0.054 

9.1749E-9 

0.160 

1.9962E-2 

0.500 

6.3494E-1 

0.056 

2.1358E-8 

0.170 

2.8858E-2 

0.520 

6.5912E-1 

0.058 

4.6745E-8 

0.180 

3.9754E-2 

0.540 

6.8146E-1 

0.060 

9.6798E-8 

0.190 

5.2613E-2 

0.560 

7.0209E-1 

0.062 

1.9069E-7 

0.200 

6.7331E-2 

0.580 

7.2116E-1 

0.064 

3.5907E-7 

0.210 

8.3750E-2 

0.600 

7.3877E-1 

0.066 

6.4902E-7 

0.220 

1.0168E-1 

0.620 

7.5507E-1 

0.068 

1.1302E-6 

0.230 

1.2091E-1 

0.660 

7.8402E-1 

0.070 

1.9025E-6 

0.240 

1.4122E-1 

0.700 

8.0885E-1 

0.072 

3.1045E-6 

0.250 

1.6239E-1 

0.740 

8.3020E-1 

0.074 

4.9236E-6 

0.260 

1.8423E-1 

0.780 

8.4861E-1 

0.076 

7.6070E-6 

0.270 

2.0653E-1 

0.820 

8.6455E-1 

0.078 

1.1473E-5 

0.280 

2.291 1E-1 

0.860 

8.7840E-1 

0.080 

1.6923E-5 

0.290 

2.5183E-1 

0.900 

8.9048E-1 

0.082 

2.4453E-5 

0.300 

2.7454E-1 

0.940 

9.0100E-1 

0.084 

3.4668E-5 

0.310 

2.9712E-1 

0.980 

9.1033E-1 

0.086 

4.8287E-5 

0.320 

3.1947E-1 

1.00 

9.1455E-1 

0.088 

6.6159E-5 

0.330 

3.4150E-1 

1.10 

9.3217E-1 

0.090 

8.9269E-5 

0.340 

3.6314E-1 

1.20 

9.4532E-1 

67 


0.092 

1.1874E-4 

0.350 

3.8432E-1 

1.30 

9.5331E-1 

0.094 

1.5586E-4 

0.360 

4.0502E-1 

1.40 

9.6304E-1 

0.096 

2.0204E-4 

0.370 

4.2518E-1 

1.50 

9.6909E-1 

0.098 

2.5885E-4 

0.380 

4.4479E-1 

1.60 

9.7390E-1 

0.100 

3.2804E-4 

0.390 

4.6382E-1 

1.70 

9.7777E-1 

0.110 

9.2957E-4 

0.400 

4.8227E-1 

1.80 

9.8091E-1 

0.120 

2.1727E-3 

0.420 

5.1738E-1 

1.90 

9.8349E-1 

0.130 

4.3866E-3 

0.440 

5.5012E-1 

2.00 

9.8563E-1 

Adapted  from  "Optoelectronics;  Theory  and  Practice"  A.  Chappell.  Ed.,  McGraw-Hill 
Book  Company,  1978 
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APPENDIX  C.  LOOK-UP  TABLE 

In-band  Radiance  (8-12um) 
and  Apparent  Temperature  Conversion  Look-up  Table 

%    Numerical  Integration  for  Planck's  law 

%    Integration  by  midpoint  rule 

%    Matlab  program  with  look-up  table  outputs 

clear; 

Tmin=291 .;    %  < —  Guess  T  to  obtain  the  corresponding  radiance 

dT=0.01; 

n=200; 

disp('  Temperature        Radiance'); 
disp(' '); 

h=6.626e-34;    c=2.998e8;     k=1.382e-23; 
a=8e-6;    b=12e-6;    dx=0.001e-6; 

for  T=Tmin:dT:Tmin+n*dT 

sfM).0; 

for  x=a:dx:b-dx 

xm=x+dx/2;  %  midpoint  of  f(x) 

fl=(2*h*cA2)/(xmA5); 

f2=l/(exp((h*c)/(xm*k*T))-l); 

fx=fl*f2; 

sf=sf+fx*dx; 

end; 

disp([T,sf]); 
end; 
end; 
(Liu,  1996) 
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Look-up  Table  for  Temperature  and  Radiance  Conversions 

Temp   Radiance  Temp   Radiance  Temp   Radiance  Temp   Radiance 


283.0000  28.8410  283.0100  28.8463  283.0200  28.8515  283.0300  28.8568 

283.0400  28.8621  283.0500  28.8673  283.0600  28.8726  283.0700  28.8779 

283.0800  28.8832  283.0900  28.8884  283.1000  28.8937  283.1100  28.8990 

283.1200  28.9042  283.1300  28.9095  283.1400  28.9148  283.1500  28.9201 

283.1600  28.9253  283.1700  28.9306  283.1800  28.9359  283.1900  28.9412 

283.2000  28.9464  283.2100  28.9517  283.2200  28.9570  283.2300  28.9623 

283.2400  28.9676  283.2500  28.9728  283.2600  28.9781  283.2700  28.9834 

283.2800  28.9887  283.2900  28.9940  283.3000  28.9993  283.3100  29.0045 

283.3200  29.0098  283.3300  29.0151  283.3400  29.0204  283.3500  29.0257 

283.3600  29.0310  283.3700  29.0363  283.3800  29.0415  283.3900  29.0468 

283.4000  29.0521  283.4100  29.0574  283.4200  29.0627  283.4300  29.0680 

283.4400  29.0733  283.4500  29.0786  283.4600  29.0839  283.4700  29.0892 

283.4800  29.0945  283.4900  29.0998  283.5000  29.1051  283.5100  29.1104 

283.5200  29.1157  283.5300  29.1210  283.5400  29.1263  283.5500  29.1316 

283.5600  29.1369  283.5700  29.1422  283.5800  29.1475  283.5900  29.1528 

283.6000  29.1581  283.6100  29.1634  283.6200  29.1687  283.6300  29.1740 

283.6400  29.1793  283.6500  29.1846  283.6600  29.1899  283.6700  29.1952 

283.6800  29.2005  283.6900  29.2058  283.7000  29.2111  283.7100  29.2164 

283.7200  29.2217  283.7300  29.2270  283.7400  29.2323  283.7500  29.2377 

283.7600  29.2430  283.7700  29.2483  283.7800  29.2536  283.7900  29.2589 

283.8000  29.2642  283.8100  29.2695  283.8200  29.2749  283.8300  29.2802 

283.8400  29.2855  283.8500  29.2908  283.8600  29.2961  283.8700  29.3014 

283.8800  29.3068  283.8900  29.3121  283.9000  29.3174  283.9100  29.3227 

283.9200  29.3280  283.9300  29.3334  283.9400  29.3387  283.9500  29.3440 

283.9600  29.3493  283.9700  29.3547  283.9800  29.3600  283.9900  29.3653 
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284.0000 
284.0400 
284.0800 
284.1200 
284.1600 
284.2000 
284.2400 
284.2800 
284.3200 
284.3600 
284.4000 
284.4400 
284.4800 
284.5200 
284.5600 
284.6000 
284.6400 
284.6800 
284.7200 
284.7600 
284.8000 
284.8400 
284.8800 
284.9200 
284.9600 
285.0000 


29.3706 
29.3920 
29.4133 
29.4346 
29.4559 
29.4773 
29.4987 
29.5200 
29.5414 
29.5628 
29.5842 
29.6056 
29.6270 
29.6485 
29.6699 
29.6913 
29.7128 
29.7343 
29.7557 
29.7772 
29.7987 
29.8202 
29.8418 
29.8633 
29.8848 
29.9064 


285.0400 
285.0800 


29.9279 
29.9495 


284.0100 
284.0500 
284.0900 
284.1300 
284.1700 
284.2100 
284.2500 
284.2900 
284.3300 
284.3700 
284.4100 
284.4500 
284.4900 
284.5300 
284.5700 
284.6100 
284.6500 
284.6900 
284.7300 
284.7700 
284.8100 
284.8500 
284.8900 
284.9300 
284.9700 
285.0100 
285.0500 
285.0900 


29.3760 
29.3973 
29.4186 
29.4399 
29.4613 
29.4826 
29.5040 
29.5254 
29.5468 
29.5681 
29.5895 
29.6110 
29.6324 
29.6538 
29.6753 
29.6967 
29.7182 
29.7396 
29.7611 
29.7826 
29.8041 
29.8256 
29.8471 
29.8687 
29.8902 
29.9118 
29.9333 
29.9549 


284.0200 
284.0600 
284.1000 
284.1400 
284.1800 
284.2200 
284.2600 
284.3000 
284.3400 
284.3800 
284.4200 
284.4600 
284.5000 
284.5400 
284.5800 
284.6200 
284.6600 
284.7000 
284.7400 
284.7800 
284.8200 
284.8600 
284.9000 
284.9400 
284.9800 
285.0200 
285.0600 
285.1000 


29.3813 
29.4026 
29.4239 
29.4453 
29.4666 
29.4880 
29.5093 
29.5307 
29.5521 
29.5735 
29.5949 
29.6163 
29.6377 
29.6592 
29.6806 
29.7021 
29.7235 
29.7450 
29.7665 
29.7880 
29.8095 
29.8310 
29.8525 
29.8740 
29.8956 
29.9171 
29.9387 
29.9603 


284.0300 
284.0700 
284.1100 
284.1500 
284.1900 
284.2300 
284.2700 
284.3100 
284.3500 
284.3900 
284.4300 
284.4700 
284.5100 
284.5500 
284.5900 
284.6300 
284.6700 
284.7100 
284.7500 
284.7900 
284.8300 
284.8700 
284.9100 
284.9500 
284.9900 
285.0300 
285.0700 
285.1100 


29.3866 
29.4079 
29.4293 
29.4506 
29.4720 
29.4933 
29.5147 
29.5361 
29.5574 
29.5788 
29.6003 
29.6217 
29.6431 
29.6645 
29.6860 
29.7074 
29.7289 
29.7504 
29.7719 
29.7934 
29.8149 
29.8364 
29.8579 
29.8794 
29.9010 
29.9225 
29.9441 
29.9657 
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285.1200 

29.9711 

285.1300 

29.9765 

285.1400 

29.9819 

285.1500 

29.9872 

285.1600 

29.9926 

285.1700 

29.9980 

285.1800 

30.0034 

285.1900 

30.0088 

285.2000 

30.0142 

285.2100 

30.0196 

285.2200 

30.0250 

285.2300 

30.0304 

285.2400 

30.0358 

285.2500 

30.0412 

285.2600 

30.0467 

285.2700 

30.0521 

285.2800 

30.0575 

285.2900 

30.0629 

285.3000 

30.0683 

285.3100 

30.0737 

285.3200 

30.0791 

285.3300 

30.0845 

285.3400 

30.0899 

285.3500 

30.0953 

285.3600 

30.1007 

285.3700 

30.1061 

285.3800 

30.1115 

285.3900 

30.1169 

285.4000 

30.1224 

285.4100 

30.1278 

285.4200 

30.1332 

285.4300 

30.1386 

285.4400 

30.1440 

285.4500 

30.1494 

285.4600 

30.1548 

285.4700 

30.1603 

285.4800 

30.1657 

285.4900 

30.1711 

285.5000 

30.1765 

285.5100 

30.1819 

285.5200 

30.1874 

285.5300 

30.1928 

285.5400 

30.1982 

285.5500 

30.2036 

285.5600 

30.2090 

285.5700 

30.2145 

285.5800 

30.2199 

285.5900 

30.2253 

285.6000 

30.2307 

285.6100 

30.2362 

285.6200 

30.2416 

285.6300 

30.2470 

285.6400 

30.2524 

285.6500 

30.2579 

285.6600 

30.2633 

285.6700 

30.2687 

285.6800 

30.2741 

285.6900 

30.2796 

285.7000 

30.2850 

285.7100 

30.2904 

285.7200 

30.2959 

285.7300 

30.3013 

285.7400 

30.3067 

285.7500 

30.3122 

285.7600 

30.3176 

285.7700 

30.3230 

285.7800 

30.3285 

285.7900 

30.3339 

285.8000 

30.3393 

285.8100 

30.3448 

285.8200 

30.3502 

285.8300 

30.3557 

285.8400 

30.3611 

285.8500 

30.3665 

285.8600 

30.3720 

285.8700 

30.3774 

285.8800 

30.3829 

285.8900 

30.3883 

285.9000 

30.3937 

285.9100 

30.3992 

285.9200 

30.4046 

285.9300 

30.4101 

285.9400 

30.4155 

285.9500 

30.4210 

285.9600 

30.4264 

285.9700 

30.4319 

285.9800 

30.4373 

285.9900 

30.4427 

286.0000 

30.4482 

286.0100 

30.4536 

286.0200 

30.4591 

286.0300 

30.4645 

286.0400 

30.4700 

286.0500 

30.4754 

286.0600 

30.4809 

286.0700 

30.4864 

286.0800 

30.4918 

286.0900 

30.4973 

286.1000 

30.5027 

286.1100 

30.5082 

286.1200 

30.5136 

286.1300 

30.5191 

286.1400 

30.5245 

286.1500 

30.5300 

286.1600 

30.5355 

286.1700 

30.5409 

286.1800 

30.5464 

286.1900 

30.5518 

286.2000 

30.5573 

286.2100 

30.5628 

286.2200 

30.5682 

286.2300 

30.5737 
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286.2400 

30.5791 

286.2500 

30.5846 

286.2600 

30.5901 

286.2700 

30.5955 

286.2800 

30.6010 

286.2900 

30.6065 

286.3000 

30.6119 

286.3100 

30.6174 

286.3200 

30.6229 

286.3300 

30.6283 

286.3400 

30.6338 

286.3500 

30.6393 

286.3600 

30.6448 

286.3700 

30.6502 

286.3800 

30.6557 

286.3900 

30.6612 

286.4000 

30.6666 

286.4100 

30.6721 

286.4200 

30.6776 

286.4300 

30.6831 

286.4400 

30.6885 

286.4500 

30.6940 

286.4600 

30.6995 

286.4700 

30.7050 

286.4800 

30.7105 

286.4900 

30.7159 

286.5000 

30.7214 

286.5100 

30.7269 

286.5200 

30.7324 

286.5300 

30.7379 

286.5400 

30.7433 

286.5500 

30.7488 

286.5600 

30.7543 

286.5700 

30.7598 

286.5800 

30.7653 

286.5900 

30.7708 

286.6000 

30.7762 

286.6100 

30.7817 

286.6200 

30.7872 

286.6300 

30.7927 

286.6400 

30.7982 

286.6500 

30.8037 

286.6600 

30.8092 

286.6700 

30.8147 

286.6800 

30.8201 

286.6900 

30.8256 

286.7000 

30.8311 

286.7100 

30.8366 

286.7200 

30.8421 

286.7300 

30.8476 

286.7400 

30.8531 

286.7500 

30.8586 

286.7600 

30.8641 

286.7700 

30.8696 

286.7800 

30.8751 

286.7900 

30.8806 

286.8000 

30.8861 

286.8100 

30.8916 

286.8200 

30.8971 

286.8300 

30.9026 

286.8400 

30.9081 

286.8500 

30.9136 

286.8600 

30.9191 

286.8700 

30.9246 

286.8800 

30.9301 

286.8900 

30.9356 

286.9000 

30.9411 

286.9100 

30.9466 

286.9200 

30.9521 

286.9300 

30.9576 

286.9400 

30.9631 

286.9500 

30.9686 

286.9600 

30.9741 

286.9700 

30.9796 

286.9800 

30.9851 

286.9900 

30.9906 

287.0000 

30.9962 

287.0000 

30.9962 

287.0100 

31.0017 

287.0200 

31.0072 

287.0300 

31.0127 

287.0400 

31.0182 

287.0500 

31.0237 

287.0600 

31.0292 

287.0700 

31.0347 

287.0800 

31.0403 

287.0900 

31.0458 

287.1000 

31.0513 

287.1100 

31.0568 

287.1200 

31.0623 

287.1300 

31.0678 

287.1400 

31.0734 

287.1500 

31.0789 

287.1600 

31.0844 

287.1700 

31.0899 

287.1800 

31.0954 

287.1900 

31.1010 

287.2000 

31.1065 

287.2100 

31.1120 

287.2200 

31.1175 

287.2300 

31.1231 

287.2400 

31.1286 

287.2500 

31.1341 

287.2600 

31.1396 

287.2700 

31.1452 

287.2800 

31.1507 

287.2900 

31.1562 

287.3000 

31.1617 

287.3100 

31.1673 

287.3200 

31.1728 

287.3300 

31.1783 

287.3400 

31.1839 
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APPENDIX  D.  AUXILIARY  PROGRAM  (  SEARAD  CODE  ) 

A.  INTRODUCTION 

"SEARAD"  is  a  FORTRAN  computer  code  that  predicts  the  radiance  (brightness) 
of  the  ocean  surface.  SEARAD"  is  valid  for  a  spectral  range  extending  from  the  visible  to 
the  far  infrared  regime.  It  is  a  self-contained,  DOS-compatible  program  that  runs  on  a 
personal  computer  and  calculates  sea  radiance.  It  is  a  modified  version  of  the  Air  Force 
program  MODTRAN2  code  that  predicts  sea  radiance  between  52.63  and  25000  cm"1. 
MODTRAN  is  a  computer  code  designed  to  determine  atmospheric  transmission  and 
radiance  at  moderate  resolution  from  0  to  50,000  cm"1  .  SEARAD  is  based  on  the  Cox- 
Munk  (Cox  and  Munk,  1954,1956)  statistical  model  for  wind-driven  capillary  wave 
facets.  Preliminary  comparisons  show  that  SEARAD  agrees  to  within  several  °C  with 
actual  sea  radiance  measurements  in  the  mid-wave  and  long-wave  in  the  infrared 
band.[Ref.l] 

B.  UNZIP  SEARAD  (  INSTALLATION) 

The  following  statement  tells  how  to  unzip  and  run  "SEARAD".  The  zipped  files 
span  three  disks  altogether.  The  last  disk  will  be  disk  3;  the  first  disk  will  be  disk  1 .  These 
commands  assume  that  the  hard  disk  in  the  computer  is  drive  c:,  and  that  the  3  1/2  inch 
disks  will  be  loaded  into  drive  a:.  If  another  drive  is  used  instead  for  the  3  1/2  disks,  the 
drive  letter  (b:  for  example)  should  be  substituted  for  a:  in  the  following  instructions. 

Sentences  in  parentheses  are  physical  acts;  sentences  without  parentheses  are 
commands  that  should  be  typed  on  the  keyboard. 
To  unzip: 

c:\  mkdir  searad 

cd  searad  (Insert  diskl  in  drive  a: ) 

copy  a:\pkunzip.exe 

pkunzip  a:\source 
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pkunzip  a:\input 

pkunzip  a:\exe     (Follow  instructions  for  disk  insertion.) 
The  files  should  all  be  unzipped  at  this  point.  The  directory  c:\searad  should  now  contain: 

zip  utility  :  pkunzip.exe 

source  code  :  Modnn  [nn=10  to  22  ]  (  except  mod21  .for  ) 

Input  file  :  Tape5xxx.std 

Log  of  mods  :  Note.txt 

Executables  :  Searad.exe,  Dirac,  and  DOSxmsf.exe 
To  run  Searad  : 

copy  tape5rad.std  tape5 

searad 

type  out 
This  should  show  the  output  file  corresponding  to  the  input  file  "  Tape5rad.std  ".[Ref.5] 

C.         SEARAD  SAMPLE  DESCRIPTION 

In  this  section  an  example  is  provided  to  show  how  to  use  SEARAD  to  predict 
ocean  radiance.  An  input  file  called  "  Tape5rad.std  "  employs  a  1976  U.S.  standard 
atmosphere  to  calculate  ocean  radiance  observed  at  a  zenith  angle  of  1 00  degree  from  a 
height  of  23  m.  The  Navy  aerosol  model  is  used.  The  calculation  is  done  with  multiple 
scattering  at  low  spectral  resolution  (LOWTRAN  7)  for  a  single  wave  number  (945  cm"1) 
in  the  long  wave  band.  The  following  DOS  commands  will  calculate  ocean  radiance  and 
print  results 

copy  tape5rad.std  tape5 

searad 

type  out 
These  commands  produce  an  output  file  called  "out".  The  four  contributions  to  ocean 
radiance  (path  to  footprint,  sea  emission,  sky  reflection,  and  sun  glint)  are  listed  at  the  end 
of  the  "out"  file.  "TBOUND"  in  the  input  file  here  is  interpreted  as  sea  temperature. 
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D.         SEARAD  MODEL 

The  SEARAD  model  computes  four  contributions  to  sea  radiance.  The 
assumption  for  the  model  is  that  the  strength  of  interaction  between  an  optical  ray  and  a 
capillary  wave  facet  is  given  by  the  facet  area  projected  normal  to  the  ray.  It  does  not 
include  multiple  reflections,  shadowing  and  gravity  waves.  It  also  ignores  polarization. 

The  first  contribution  is  path  radiance  (Afpath).  The  footprint  of  a  single  pixel  in  an 
image  of  the  sea  is  indicated  by  the  wavy  line.  The  footprint  is  observed  by  a  receiver  at 
the  end  of  a  ray  whose  zenith  angle  at  the  footprint  is  0r  .  Np2Lth  represents  the  spectral 
radiance  in  W/m2  /sr/cm  along  the  path  from  the  footprint  to  the  receiver. 

The  second  contribution  is  reflected  sky  radiance  (Nsky).  The  spectral  radiance  Ns 
from  part  of  the  sky  reaches  the  footprint  along  a  ray  whose  zenith  angle  is  9S.  The 
contributions  are  the  summation  from  all  portions  of  the  sky  after  specular  reflection  by 
the  appropriate  facets  within  the  footprint.  Nsky  is  the  sum  leaving  the  footprint  at  zenith 
angle  ds  .  During  its  path  to  the  receiver,  Nsky  is  attenuated  by  the  path  transmission  rpath. 

The  third  contribution  is  reflected  solar  radiance  (  Nsm  ),  sun  glint  .  The  radiance 
N0  from  the  solar  center  arrives  at  the  footprint  along  a  path  whose  zenith  angle  is  60. 
Within  the  footprint  most  facets  deflect  the  solar  ray  away  from  the  receiver  and  are 
rejected,  but  some  facets  are  retained  because  they  deflect  the  ray  specularly  toward  the 
receiver  along  a  path  with  zenith  angle  9r.  Nsun  is  the  spectral  radiance  leaving  the 
footprint  after  summation  over  rays  arriving  from  all  portions  of  the  solar  disk.  Nsun  is 
also  attenuated  by  the  path  transmission  rpath. 

The  fourth  contribution  is  thermal  black  body  emission  (  Nsea  ).  In  this  portion 
each  facet  emits  a  spectral  radiance  Nbh  given  by  Planck's  equation  for  a  black  body.  The 
sea  temperature  is  equal  to  the  value  of  TBOUND  in  the  input  file.  The  spectral 
emissivity  of  a  given  facet  in  the  direction  of  the  receiver  is  specified  by  the  slope  of  that 
facet  and  the  value  of  6r  Nsea  is  the  thermal  spectral  radiance  leaving  the  footprint  for  the 
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receiver  after  summation  over  all  facets  within  the  footprint.  7Vsea  is  also  attenuated  by 
path  transmission  rpath    after  leaving  the  footprint. 

The  symbol  p  represents  the  reflectivity  of  sea  water.  It  is  calculated  from 
Fresnel's  equations.  The  total  spectral  radiance  N(v)  received  at  wave  number  v  (cm1) 
is  given  by 

N(v)  =  Npalh(v)f(v)  +  [nJv)  +  N^M  +  NsJv)]rpalh(v)f(v)  ,  (C.l) 

where  f(v)  stands  for  the  spectral  responsivity  of  the  receiver. 

The  design  of  SEARAD  is  such  that  the  path  (Npath  ,  tpath)  and  source  (Ns,  N0, 
Nbb)  values  are  taken  from  the  original  MODTRAN  2  while  Fresnel  reflection  (p)  and 
slope  integrated  values  (Nsky  ,Nsun,  NseJ  are  introduced  in  new  subroutines.  Integration  of 
Equation  (C.l)  over  the  wave  number  band  specified  in  the  input  file  (Tape5)  is 
implemented  in  a  modification  of  subroutine  "TRAN"  to  produce  the  band  integrated 
values  for  sea  radiance  given  in  the  output  file  (out).[Ref.l] 
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APPENDIX  E.  DESCRIPTION  OF  EACH  CARD  IN  SEARAD  INPUT  FILE 

The  format  of  input  file  is  very  important.  Each  number  must  be  in  the  correct  row 
and  column.  Below  is  a  listing  of  the  input  file  parameters  in  SEARAD  and  the  values 
used  in  evaluating  the  sea  radiance. 
CARD  1  (  line  one  ) 
F       =>    useLOWTRAN7 
7       =>    MODEL  :  radiosonde  data  used 
3       =>    I  TYPE  :  slant  path  to  space 
1       =>    IEMSCT  :  program  execution  in  radiance  mode 
1       =>    IMULT  :  program  executed  with  multiple  scattering 

"M1,M2,M3"    are   used   to    modify   or   supplement   the   altitude   profiles    of 
temperature  and  pressure,  water  vapor,  and  ozone  "M4,M5,M6" 
6        =>    Ml  :  default  temperature  and  pressure  to  specified  MODEL 

atmosphere  (Ml =1-6) 
6        =>    M2  :  default  H20  to  specified  MODEL  atmosphere  (  M2=l~6  ) 

6        =>    M3  :  default  ozone  to  specified  MODEL  atmosphere  (  M3=l~6  ) 
6        =>    M4  :  default  CH4  to  specified  MODEL  atmosphere  (  M4=l~6  ) 
6        =>    M5  :  default  N20  to  specified  MODEL  atmosphere  (  M5=l~6  ) 

6        =>    M6  :  default  CO  to  specified  MODEL  atmosphere  (  M6=l~6  ) 

1        =>    MDEF  use  default  profile  for  C02,  02,  NO,  S02,  N02,  NH3,  HN03 

1        =>    IM  :  radiosonde  data  to  be  read  initially 

1        =>    NORPT  :  minimize  printing  of  transmittance 

xx      =>    TBOUND  :  boundary  temperature  ( °  k  ) 

0.00    =5-    SALB  :  surface  albedo  (  0.0  to  1.0  ) 
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CARD  2  (line  2) 

3        =>  IHAZE  :  navy  maritime  extinction,  sets  own  VIS. 

0        =>  ISEASN  :  Summer  for  model  7 

0        =>  IVULCN  :  default  to  stratospheric  background 

3        =>  ICSTL  :  air  mass  (1  for  open  ocean  -  10  for  strong  continental 

influence)  ICSTL  now  allow  non-integer  values.  Integers  from 
1  to  10  inclusive  will  be  taken  as  is,  while  integers  from  1 1  to  300 
inclusive  will  be  divided  by  10.  Integers  outside  the  range  of  1  to 
300  will  be  reset  to  a  default  value  of  "3.0". 

0        =>  ICLD  :  no  clouds  or  rain 

0        =>  IVSA:  not  used 

20      =>  VIS  :  meteorological  range 

1 .897=>  WSS  :  specifies  current  wind  speed  (  m/s  ) 

1 .603  =>  WHH  :  24  hour  average  wind  speed  (  m/s  ) 

0.00  =>  RAINT  :  rain  rate  (  mm/hr ) 

0.00  =>  GND ALT  altitude  of  surface  relative  to  sea  level  (  km  ) 

CARD  2c  ( line  3  ) 

33       =>  ML  :  number  of  atmospheric  layers  to  be  inserted  (Maximum  of  34) 

0         =>  IRD1  :  no  read 

0         =>  IRD2  :  no  read 

text  TITLE  :  user  defined  text 

CARD  2c  1  (next  33  line) 

xx       =>  ZMDL  :  altitude  of  layer  (  km  ) 

xx      =>  P  :  pressure  at  layer  (  mb  ) 

xx      =>  T  ambient  temperature  ( °C  ) 
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xx      =>  JCHAR(RH)  :  relative  humidity  (  %  ) 

ABH=>  set  units  for  P,T,RH 

CARD  3  ( line  37  ) 

0.006=>  HI  :  initial  altitude  (FLIR  sensor  altitude)  (  km  ) 

0.00  =>  H2  :  final  altitude  (seasurface  altitude)  (  km  ) 

xx      =>  ANGLE  :  zenith  angle  (  degree  ) 

0.000=>  RANGE  :  path  length  (  km  ) 

0.000=>  BETA  :  earth  center  angle  subtented  by  HI  and  H2  (  degree  ) 

0.000=>  RO  :  default  radius  of  earth  (  km  ) 

1  =>  LEN  for  long  path  through  tangent  height 

90.00=>  Psi  :  azimuth  of  the  upwind  direction  from  line  of  sight  positive  East  of  North 
(degree)  90  °  means  that  the  wind  is  blowing  from  right  to  left, 
perpendicular  to  the  direction  of  observation 

T  =>  True  :  sea  radiance  will  be  carried  out  (  F  means  that  sea  radiance  will  be 
prevented  under  all  conditions  ) 

CARD  4  (  line  38  ) 

830     =>  VI  :  initial  wave  number  (  cm1  ) 

1250  =>  V2  :  final  wave  number  (  cm-1 ) 

10       =>  DV  :  wave  number  increment  (  cm-1  ) 

2  =>  filter 
CARD  5 

0         =>  IRPT  :  end  program  (  3  means  read  CARD  3  again  ) 
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APPENDIX  F.  SEARAD  INPUT  DATA 

EOMET95  Data  Set  from  NPS  Boundary  Layer  Meteorology  Group 

Header  Entries  form  Appendix  E 


7     3 

1     1 

6     6 

6     6 

6     6 

111 

286.60 

0.00 

3     0 

0     3 

0     0 

10.000 

1.897 

1.603 

000 

,000 

33     0 

0      11: 

:37:53  MAY 

1995 

0.003 

1013.400 

12.100 

90.0 

ABH 

0.017 

1011.700 

13.300 

70.0 

ABH 

0.021 

1011.300 

13.000 

72.0 

ABH 

0.031 

1010.100 

12.900 

74.0 

ABH 

0.042 

1008.800 

12.800 

75.0 

ABH 

0.046 

1008.300 

12.800 

76.0 

ABH 

0.056 

1007.100 

12.600 

77.0 

ABH 

0.070 

1005.300 

12.500 

78.0 

ABH 

0.075 

1004.900 

12.500 

78.0 

ABH 

0.083 

1003.900 

12.500 

77.0 

ABH 

0.093 

1002.600 

12.500 

75.0 

ABH 

0.103 

1001.400 

12.500 

74.0 

ABH 

0.113 

1000.200 

12.500 

73.0 

ABH 

0.119 

999.500 

12.400 

73.0 

ABH 

0.128 

998.500 

12.400 

73.0 

ABH 

0.137 

997.400 

12.500 

73.0 

ABH 

0.150 

995.800 

12.500 

73.0 

ABH 

0.153 

995.500 

12.700 

72.0 

ABH 

0.160 

994.600 

12.800 

70.0 

ABH 

0.164 

994.100 

12.800 

70.0 

ABH 

0.185 

991.700 

12.800 

69.0 

ABH 

0.190 

991.100 

12.800 

69.0 

ABH 

0.201 

989.800 

12.800 

69.0 

ABH 

0.208 

988.900 

12.700 

69.0 

ABH 

0.218 

987.800 

12.600 

69.0 

ABH 

0.228 

986.600 

12.600 

69.0 

ABH 

0.234 

985.900 

12.500 

69.0 

ABH 

0.244 

984.700 

12.400 

69.0 

ABH 

0.256 

983.300 

12.300 

68.0 

ABH 

0.266 

982.100 

12.300 

69.0. 

ABH 

0.275 

981.100 

12.100 

69.0 

ABH 

0.283 

980.100 

12.100 

69.0 

ABH 

0.291 

979.200 

12.100 

69.0 

ABH 

0.006 

0.000 

90.240 

.000 

.000 

0.00     1 

252.000 

T 

830 
0 

1250 

10 

1 
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7     3 

1     1 

6     6 

6     6 

6     6 

1     1     1 

286.60 

0.00 

3     0 

0     3 

0     0 

10.000 

1.897 

1.603 

000 

000 

33     0 

0      typ 

ical  case 

0.003 

1013.400 

12.100 

90.0 

ABH 

0.017 

1011.700 

13.300 

70.0 

ABH 

0.021 

1011.300 

13.000 

72.0 

ABH 

0.031 

1010.100 

12.900 

74.0 

ABH 

0.042 

1008.800 

12.800 

75.0 

ABH 

0.046 

1008.300 

12.800 

76.0 

ABH 

0.056 

1007.100 

12.600 

77.0 

ABH 

0.070 

1005.300 

12.500 

78.0 

ABH 

0.075 

1004.900 

12.500 

78.0 

ABH 

0.083 

1003.900 

12.500 

77.0 

ABH 

0.093 

1002.600 

12.500 

75.0 

ABH 

0.103 

1001.400 

12.500 

74.0 

ABH 

0.113 

1000.200 

12.500 

73.0 

ABH 

0.119 

999.500 

12.400 

73.0 

ABH 

0.128 

998.500 

12.400 

73.0 

ABH 

0.137 

997.400 

12.500 

73.0 

ABH 

0.150 

995.800 

12.500 

73.0 

ABH 

0.153 

995.500 

12.700 

72.0 

ABH 

0.160 

994.600 

12.800 

70.0 

ABH 

0.164 

994.100 

12.800 

70.0 

ABH 

0.185 

991.700 

12.800 

69.0 

ABH 

0.190 

991.100 

12.800 

69.0 

ABH 

0.201 

989.800 

12.800 

69.0 

ABH 

0.208 

988.900 

12.700 

69.0 

ABH 

0.218 

987.800 

12.600 

69.0 

ABH 

0.228 

986.600 

12.600 

69.0 

ABH 

0.234 

985.900 

12.500 

69.0 

ABH 

0.244 

984.700 

12.400 

69.0 

ABH 

0.256 

983.300 

12.300 

68.0 

ABH 

0.266 

982.100 

12.300 

69.0 

ABH 

0.275 

981.100 

12.100 

69.0 

ABH 

0.283 

980.100 

12.100 

69.0 

ABH 

0.291 

979.200 

12.100 

69.0 

ABH 

0.006 

0.000 

90.043 

.000 

.000 

0.00     1 

252.000 

T 

830 

0 

1250 

10 

1 

84 


7     3 

1     1 

6     6 

6     6 

6     6 

111   286.60    0.00 

3     0 

0     3 

0     0 

10.000 

1.897 

1.603       .000       .000 

33     0 

0     MDR 

Estimate-1 

0.003 

1013.400 

12.100 

90.0 

ABH 

0.017 

1011.700 

13.300 

70.0 

ABH 

0.021 

1011.300 

13.000 

72.0 

ABH 

0.031 

1010.100 

12.900 

74.0 

ABH 

0.042 

1008.800 

12.800 

75.0 

ABH 

0.046 

1008.300 

12.800 

76.0 

ABH 

0.056 

1007.100 

12.600 

77.0 

ABH 

0.070 

1005.300 

12.500 

78.0 

ABH 

0.075 

1004.900 

12.500 

78.0 

ABH 

0.083 

1003.900 

12.500 

77.0 

ABH 

0.093 

1002.600 

12.500 

75.0 

ABH 

0.103 

1001.400 

12.500 

74.0 

ABH 

0.113 

1000.200 

12.500 

73.0 

ABH 

0.119 

999.500 

12.400 

73.0 

ABH 

0.128 

998.500 

12.400 

73.0 

ABH 

0.137 

997.400 

12.500 

73.0 

ABH 

0.150 

995.800 

12.500 

73.0 

ABH 

0.153 

995.500 

12.700 

72.0 

ABH 

0.160 

994.600 

12.800 

70.0 

ABH 

0.164 

994.100 

12.800 

70.0 

ABH 

0.185 

991.700 

12.800 

69.0 

ABH 

0.190 

991.100 

12.800 

69.0 

ABH 

0.201 

989.800 

12.800 

69.0 

ABH 

0.208 

988. 900 

12.700 

69.0 

ABH 

0.218 

987.800 

12.600 

69.0 

ABH 

0.228 

986.600 

12.600 

69.0 

ABH 

0.234 

985.900 

12.500 

69.0 

ABH 

0.244 

984.700 

12.400 

69.0 

ABH 

0.256 

983.300 

12.300 

68.0 

ABH 

0.266 

982.100 

12.300 

69.0 

ABH 

0.275 

981.100 

12.100 

69.0 

ABH 

0.283 

980.100 

12.100 

69.0 

ABH 

0.291 

979.200 

12.100 

69.0 

ABH 

0.250 

0.000 

95.000 

.000 

.000 

0.00     1    252.000     T 

830 

0 

1250 

10 

1 

85 


7     3 

1     1 

6     6 

6     6 

6     6 

111   286.60    0.00 

3     0 

0     3 

0     0 

10.000 

1.897 

1.603       .000       .000 

3     0 

0     MDR 

Estimate-2 

0.003 

1013.400 

12.100 

90.0 

ABH 

0.017 

1011.700 

13.300 

70.0 

ABH 

0.021 

1011.300 

13.000 

72.0 

ABH 

0.031 

1010.100 

12.900 

74.0 

ABH 

0.042 

1008.800 

12.800 

75.0 

ABH 

0.046 

1008.300 

12.800 

76.0 

ABH 

0.056 

1007.100 

12.600 

77.0 

ABH 

0.070 

1005.300 

12.500 

78.0 

ABH 

0.075 

1004. 900 

12.500 

78.0 

ABH 

0.083 

1003.900 

12.500 

77.0 

ABH 

0.093 

1002.600 

12.500 

75.0 

ABH 

0.103 

1001.400 

12.500 

74.0 

ABH 

0.113 

1000.200 

12.500 

73.0 

ABH 

0.119 

999.500 

12.400 

73.0 

ABH 

0.128 

998.500 

12.400 

73.0 

ABH 

0.137 

997.400 

12.500 

73.0 

ABH 

0.150 

995.800" 

12.500 

73.0 

ABH 

0.153 

995.500 

12.700 

72.0 

ABH 

0.160 

994.600 

12,800 

70.0 

ABH 

0.164 

994.100 

12.800 

70.0 

ABH 

0.185 

991.700 

12.800 

69.0 

ABH 

0.190 

991.100 

12.800 

69.0 

ABH 

0.201 

989.800 

12.800 

69.0 

ABH 

0.208 

988.900 

12.700 

69.0 

ABH 

0.218 

987.800 

12.600 

69.0 

ABH 

0.228 

986.600 

12.600 

69.0 

ABH 

0.234 

985.900 

12.500 

69.0 

ABH 

0.244 

984.700 

12.400 

69.0 

ABH 

0.256 

983.300 

12.300 

68.0 

ABH 

0.266 

982.100 

12.300 

69.0 

ABH 

0.275 

981.100 

12.100 

69.0 

ABH 

0.283 

980.100 

12.100 

69.0 

ABH 

0.291 

979.200 

12.100 

69.0 

ABH 

0.250 

0.000 

90.550 

.000 

.000 

0.00     1    252.000     T 

830 

0 

1250 

10 

1 

86 


APPENDIX  G.  SEARAD  OUTPUT  DATA 


DATE:  11/21/1997 

THERMAL  RADIANCE  MODE 

MULTIPLE  SCATTERING  USED 

MARINE  AEROSOL  MODEL  USED 

WIND  SPEED 
WIND  SPEED 
RELATIVE  HUMIDITY  = 
AIRMASS  CHARACTER  = 
VISIBILITY 

SLANT  PATH  TO  SPACE 


*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 

TIME:  12:00:50.69 


1.90  M/SEC 

1.60  M/SEC,  24  HR  AVERAGE 
8  9.93  PERCENT 

3 
10.00  KM 


HI 

= 

.006 

KM 

HMIN 

= 

.000 

KM 

ANGLE 

= 

90.240 

DEG 

FREQUENCY  RANGE 


IV1 

= 

.  830  CM-1 

IV2 

= 

1250  CM-1 

IDV 

= 

10  CM-1 

IFWHM 

= 

1  CM-1 

IFILTER 

= 

0 

12.05  MICROMETERS) 
8.00  MICROMETERS) 


SUMMARY  OF  THE  GEOMETRY  CALCULATION 


HI 

= 

.006 

KM 

H2 

= 

.000 

KM 

ANGLE 

= 

90 

.240 

DEG 

RANGE 

= 

1 

.444 

KM 

BETA 

= 

.013 

DEG 

PHI 

= 

89 

.764 

DEG 

HMIN 

= 

.000 

KM 

BENDING 

= 

.009 

DEG 

LEN 

= 

0 

SEA  AT  2  8  6.60  K  REPLACES  BLACK  BODY  BOUNDARY 

UPWIND   =     252.000  DEG  EAST  OF  LINE  OF  SIGHT 

ZERO  RANGE  VALUES 

SEA  EMISSION  =   14.42182  W  M-2  SR-1   (AV.  EMISS 


4636' 


SKY  REFLECTION 
SUN  GLINT 


5.22488  W  M-2  SR-1 
.00000  W  M-2  SR-1 
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TOTAL  RADIANCE 
BLACK  BODY  TEMP. 

FULL  RANGE  VALUES 

PATH  TO  FOOTPRINT 


.7181) 


SEA  EMISSION 
SKY  REFLECTION 
SUN  GLINT 

TOTAL  RADIANCE 
BLACK  BODY  TEMP, 


19. 64670  W  M-2  SR-1 
-10.1  C 


7.95021  W  M-2  SR-1 

10.50470  W  M-2  SR-1 

3.70776  W  M-2  SR-1 

.00000  W  M-2  SR-1 

22.16267  W  M-2  SR-1 

-4.3  C 


[AV.  TRANS 
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DATE:  11/21/1997 

THERMAL  RADIANCE  MODE 

MULTIPLE  SCATTERING  USED 

MARINE  AEROSOL  MODEL  USED 

WIND  SPEED 
WIND  SPEED 
RELATIVE  HUMIDITY 
AIRMASS  CHARACTER 
VISIBILITY 

SLANT  PATH  TO  SPACE 


*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 

TIME:  12:03:37.72 


1. 90  M/SEC 

1.60  M/SEC,  24  HR  AVERAGE 

8  9.93  PERCENT 

3 

10.00  KM 


HI 

= 

.006  KM 

HMIN 

= 

.000  KM 

ANGLE 

= 

90.043  DEG 

FREQUENCY  RANGE 


IV1 

830 

CM-1 

IV2 

1250 

CM-1 

IDV 

10 

CM-1 

IFWHM 

1 

CM-1 

IFILTER  = 

0 

12.05  MICROMETERS) 
8.00  MICROMETERS) 


SUMMARY  OF  THE  GEOMETRY  CALCULATION 


HI 

= 

.006 

KM 

H2 

= 

.291 

KM 

ANGLE 

= 

90 

.043 

DEG 

RANGE 

= 

97 

.261 

KM 

BETA 

= 

.875 

DEG 

PHI 

= 

90 

.487 

DEG 

HMIN 

= 

.000 

KM 

BENDING 

= 

.344 

DEG 

LEN 

= 

1 

TBOUND  SET  TO 


.10  K  FOR  MARINE  SKY 


INTEGRATED  ABSORPTION 
AVERAGE  TRANSMITTANCE 
MAXIMUM  RADIANCE 
MINIMUM  RADIANCE 
BOUNDARY  TEMPERATURE 
BOUNDARY  EMISSIVITY 
FILTERED  RADIANCE 
BLACKBODY  TEMPERATURE 


420.00  CM-1  FROM   830  TO  1250  CM-1 

.0000 

1. 

.033E-01  W  M-2  SR-1  (CM-l)-l  AT    830.0  CM-1 

4 

.248E-02  W  M-2  SR-1  (CM-l)-l  AT   1250.0  CM-1 

.10  K 

1.000 

2 

. 904E+01  W  M-2  SR-1 

9.8  C 
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*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 

TIME:  12:06:44.19 


DATE:  11/21/1997 

THERMAL  RADIANCE  MODE 

MULTIPLE  SCATTERING  USED 

MARINE  AEROSOL  MODEL  USED 

WIND  SPEED 
WIND  SPEED 
RELATIVE  HUMIDITY  = 
AIRMASS  CHARACTER  = 
VISIBILITY 

SLANT  PATH  TO  SPACE 


1. 90  M/SEC 

1.60  M/SEC,  24  HR  AVERAGE 

8  9.93  PERCENT 

3 

10.00  KM 


HI 

= 

.250 

KM 

HMIN 

= 

.000 

KM 

ANGLE 

= 

95.000 

DEG 

FREQUENCY  RANGE 


IV1 

= 

830 

CM-1 

IV2 

= 

1250 

CM-1 

IDV 

= 

10 

CM-1 

IFWHM 

= 

1 

CM-1 

IFILTER 

= 

0 

12.05  MICROMETERS) 
8.00  MICROMETERS) 


SUMMARY  OF  THE  GEOMETRY  CALCULATION 


HI 

= 

.250 

KM 

H2 

= 

.000 

KM 

ANGLE 

= 

95 

.000 

DEG 

RANGE 

= 

2 

.875 

KM 

BETA 

= 

.026 

DEG 

PHI 

= 

85 

.020 

DEG 

HMIN 

= 

.000 

KM 

BENDING 

= 

.006 

DEG 

LEN 

= 

0 

SEA  AT  286.60  K  REPLACES  BLACK  BODY  BOUNDARY 

UPWIND   =     252.000  DEG  EAST  OF  LINE  OF  SIGHT 
ZERO  RANGE  VALUES 


.5882 


SEA  EMISSION 

SKY  REFLECTION 
SUN  GLINT 

TOTAL  RADIANCE 
BLACK  BODY  TEMP, 


=   18.29033  W  M-2  SR-1   (AV.  EMISS. 

3.71435  W  M-2  SR-1 
.00000  W  M-2  SR-1 

=   22.00468  W  M-2  SR-1 
-4.7  C 


FULL  RANGE  VALUES 


90 


6154-) 


PATH  TO  FOOTPRINT 

SEA  EMISSION 
SKY  REFLECTION 
SUN  GLINT 


=   10. 91463  W  M-2  SR-1 

=   11.46502  W  M-2  SR-1 

2.24949  W  M-2  SR-1 

. 00000  W  M-2  SR-1 


(AV.  TRANS 


TOTAL  RADIANCE 
BLACK  BODY  TEMP. 


=   24 . 62915  W  M-2  SR-1 
1.1  C 
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DATE:  11/21/1997 

THERMAL  RADIANCE  MODE 

MULTIPLE  SCATTERING  USED 

MARINE  AEROSOL  MODEL  USED 

WIND  SPEED 
WIND  SPEED 
RELATIVE  HUMIDITY  = 
AIRMASS  CHARACTER  = 
VISIBILITY 

SLANT  PATH  TO  SPACE 


*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 

TIME:  12:09:01.17 


1.90  M/SEC 

1.60  M/SEC,  24  HR  AVERAGE 

8  9.93  PERCENT 

3 

10.00  KM 


HI 

= 

.250  KM 

HMIN 

= 

.000  KM 

ANGLE 

= 

90.550  DEG 

FREQUENCY  RANGE 


IV1 

= 

830 

CM-1 

IV2 

= 

1250 

CM-1 

IDV 

= 

10 

CM-1 

IFWHM 

= 

1 

CM-1 

IFILTER 

= 

0 

12.05  MICROMETERS) 
8.00  MICROMETERS) 


SUMMARY  OF  THE  GEOMETRY  CALCULATION 


HI 

= 

.250 

KM 

H2 

= 

.000 

KM 

ANGLE 

= 

90 

.550 

DEG 

RANGE 

= 

33 

.415 

KM 

BETA 

= 

.300 

DEG 

PHI 

= 

89 

.683 

DEG 

HMIN 

= 

.000 

KM 

BENDING 

= 

.067 

DEG 

LEN 

= 

0 

SEA  AT  286.60  K  REPLACES  BLACK  BODY  BOUNDARY 

UPWIND   =     252.000  DEG  EAST  OF  LINE  OF  SIGHT 
ZERO  RANGE  VALUES 


.4649! 


SEA  EMISSION 

SKY  REFLECTION 
SUN  GLINT 

TOTAL  RADIANCE 
BLACK  BODY  TEMP, 


=   14.46294  W  M-2  SR-1   (AV.  EMISS 

5.26282  W  M-2  SR-1 
.00000  W  M-2  SR-1 

=   19.72575  W  M-2  SR-1 
-9.9  C 


FULL  RANGE  VALUES 


92 


0171 


PATH  TO  FOOTPRINT 

SEA  EMISSION 
SKY  REFLECTION 
SUN  GLINT 


=   28.50598  W  M-2  SR-1   (AV.  TRANS 

.24282  W  M-2  SR-1 
.07884  W  M-2  SR-1 
.00000  W  M-2  SR-1 


TOTAL  RADIANCE 
BLACK  BODY  TEMP. 


28.82764  W  M-2  SR-1 
9.5  C 
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APPENDIX  H.  GPS  DATA 


Source:  EOMET95  Data  Set 


ws:  wind  speed 

ta:  air  temperature 

lat: 

latitude 

wd:  wind  direction 

ts:  sea  temperature 

Ion: 

longitudi 

•*■ 

shs:  ship  seep 

rh:  relative  humidity 

shd:  ship  direction 

pr:  pressure 

date 

time 

ws 

wd 

ta 

rh 

Pr 

shs 

shd 

lat 

Ion 

ts 

950515 

2030 

1.00 

153 

14.55 

70.40 

1010 

0. 

01 

264 

36. 

80 

-121. 

79 

16. 

90 

950515 

2040 

0.45 

204 

13.98 

72.84 

1010 

0. 

02 

30 

36. 

80 

-121. 

79 

16. 

50 

950515 

2050 

0.01 

309 

14.22 

72.70 

1011 

0. 

02 

48 

36. 

80 

-121. 

79 

17. 

60 

950515 

2100 

0.64 

243 

15.02 

69.92 

1011 

0. 

01 

300 

36. 

80 

-121. 

79 

18. 

20 

950515 

2110 

0.24 

274 

15.74 

70.76 

1011 

0. 

05 

37 

36. 

80 

-121. 

79 

19. 

10 

950515 

2120 

2.14 

275 

16.28 

68.09 

1012 

1. 

24 

340 

36. 

81 

-121. 

79 

14. 

20 

950515 

2130 

5.17 

240 

14.77 

76.57 

1010 

4. 

77 

247 

36. 

80 

-121. 

80 

12. 

30 

950515 

2140 

4.25 

248 

13.91 

82.30 

1011 

5. 

11 

254 

36. 

80 

-121. 

83 

13. 

10 

950515 

2150 

6.64 

258 

13.06 

87.38 

1011 

4. 

84 

257 

36. 

79 

-121. 

87 

13. 

40 

950515 

2200 

6.12 

261 

13.05 

87.31 

1011 

4. 

71 

256 

36. 

78 

-121. 

90 

12. 

40 

950515 

2210 

6.09 

274 

13.19 

85.85 

1011 

4. 

56 

254 

36. 

78 

-121. 

93 

12. 

30 

950515 

2220 

5.66 

289 

13.01 

86.30 

1011 

4. 

53 

250 

36. 

77 

-121. 

96 

12. 

70 

950515 

2230 

5.46 

299 

13.07 

83.25 

1011 

4. 

46 

250 

36. 

76 

-121. 

98 

13. 

00 

950515 

2240 

4.68 

304 

13.30 

81.02 

1011 

2. 

59 

257 

36. 

75 

-122. 

01 

13. 

50 

950515 

2250 

3.94 

300 

13.52 

82.33 

1011 

0. 

50 

275  ■ 

36. 

75 

-122. 

02 

13. 

80 

950515 

2300 

2.98 

298 

13.48 

84.72 

1011 

0. 

51 

282 

36. 

75 

-122. 

02 

13. 

90 

950515 

2310 

2.68 

302 

13.40 

84.98 

1011 

0. 

59 

284 

36. 

76 

-122. 

02 

13. 

,90 

950515 

2320 

2.62 

314 

13.47 

84.37 

1011 

0. 

58 

279 

36. 

76 

-122. 

,03 

13. 

,90 

950515 

2330 

2.74 

331 

13.53 

84.75 

1011 

0. 

.59 

280 

36. 

,76 

-122. 

.03 

13. 

,90 

950515 

2340 

2.90 

308 

13.54 

85.24 

1011 

0. 

,55 

291 

36. 

,76 

-122. 

,04 

14. 

,10 

950515 

2350 

3.00 

294 

13.48 

86.07 

1011 

0. 

,49 

293 

36, 

,76 

-122. 

,04 

13, 

,90 

950516 

0 

3.01 

276 

13.23 
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.16 

88 

.63 

1011 

2 

.69 

128 

36 

.77 

-122 

.05 

13 

.00 

950524 

1200 

4 

.85 

287 

11 

.16 

88 

.23 

1011 

2 

.32 

130 

36 

.77 

-122 

.04 

12 

.90 

950524 

1210 

5 

.27 

283 

11 

.17 

88 

.74 

1011 

2 

.23 

130 

36 

.76 

-122 

.03 

12 

.90 

950524 

1220 

4 

.22 

285 

11 

.19 

88 

.82 

1011 

1 

.90 

134 

36 

.75 

-122 

.01 

13 

.10 

950524 

1230 

5 

.74 

310 

11 

.19 

89 

.80 

1011 

0 

.65 

320 

36 

.75 

-122 

.01 

13 

.00 

950524 

1240 

6 

.22 

309 

11 

.12 

88 

.37 

1011 

0 

.59 

352 

36 

.75 

-122 

.01 

13 

.00 

950524 

1250 

4 

.87 

291 

11 

.16 

87 

.92 

1011 

4 

.56 

73 

36 

.76 

-122 

.00 

13 

.10 

950524 

1300 

4 

.93 

276 

11 

.14 

88 

.38 

1011 

4 

.63 

74 

36 

.76 

-121 

.97 

13 

.40 

950524 

1310 

5 

.20 

272 

11 

.16 

88 

.67 

1011 

4 

.68 

76 

36 

.77 

-121 

.94 

13 

.40 

950524 

1320 

4 

.93 

265 

11 

.12 

88 

.63 

1011 

4 

.72 

76 

36 

.78 

-121 

.91 

13 

.30 

950524 

1330 

5 

.25 

267 

11 

.12 

89 

.04 

1011 

4 

.68 

75 

36 

.78 

-121 

.88 

13 

.40 

950524 

1340 

4 

.83 

257 

11 

.10 

89 

.62 

1011 

4 

.63 

74 

36 

.79 

-121 

.85 

13 

.60 

950524 

1350 

4 

.72 

252 

11 

.24 

89 

.47 

1011 

4 

.74 

71 

36 

.80 

-121 

.82 

13 

.60 

950524 

1400 

4 

.82 

253 

11 

.17 

90 

.23 

1011 

2 

.76 

67 

36 

.81 

-121 

.79 

13 

.30 

950524 

1410 

1 

.39 

281 

11 

.01 

91 

.94 

1011 

0 

.46 

196 

36 

.80 

-121 

.79 

13 

.30 

950524 

1420 

1 

.05 

275 

11 

.04 

90 

.94 

1011 

0 

.00 

242 

36 

.80 

-121 

.79 

13 

.00 

113 
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APPENDIX  I.  RADIOSONDE  METEOROLOGICAL  DATA 


Sounding  program  REV  7 . 62  using  Omega 

Ship      : PT_SUR 

Location  :  36.80  N  122.00  W       3m 

Rejected  Sigma  stations:  a,b,c 

Phase  fitting  length  is  250  s  from    0  min  to  120  min 

Sounding  :  4 

RS-number:  158735751 

No  PTU  editing 

No  wind  editing 

Started  at:       16  MAY  95  14:59  GMT 


e 

AscRat 

e  Hgt/IY 

ISL  Pressure    Temp   RH 

Dewp   Dir  Speed 

WndStat 

n 

s 

m/s 

m 

hPa 

degC 

Q. 
O 

degC 

deg 

m/s 

0 

0 

0.0 

3 

1013.8 

11.4 

94 

10.5 

36 

2.5 

0 

5 

7.8 

42 

1009.0 

11.7 

75 

7.5 

/// 

//// 

— CD-FGH- 

0 

10 

5.7 

60 

1006.8 

11.5 

77 

7.7 

345 

1.  9 

— CD-FGH- 

0 

15 

4.9 

82 

1004.1 

11.6 

72 

6.8 

/// 

//// 

— CD-FGH- 

0 

20 

6.8 

139 

997.3 

12.1 

68 

6.4 

342 

2.4 

— CD-FGH- 

0 

25 

6.1 

159 

995.0 

12.2 

60 

4.7 

/// 

//// 

— CD-FGH- 

0 

30 

6.0 

184 

991.9 

12.3 

60 

4.8 

345 

2.5 

— CD-FGH- 

0 

35 

5.7 

216 

988.2 

12.2 

62 

5.2 

/// 

//// 

— CD-FGH- 

0 

40 

6.3 

248 

984.5 

12.1 

61 

4.9 

345 

2.9 

--CD-FGH- 

0 

45 

6.6 

276 

981.1 

11.9 

62 

4.9 

/// 

//// 

--CD-FGH- 

0 

50 

5.7 

309 

977.3 

11.7 

63 

4.9 

344 

3.1 

— CD-FGH- 

0 

55 

6.0 

336 

974.1 

11.6 

64 

5.1 

/// 

//// 

— CD-FGH- 

1 

0 

6.0 

364 

970.9 

11.6 

64 

5.1 

342 

3.6 

— CD-FGH- 

1 

5 

5.7 

390 

967.8 

11.7 

63 

4.9 

/// 

//// 

— CD-FGH- 

1 

10 

6.0 

427 

963.5 

11.6 

63 

4.9 

343 

3.2 

--CD-FGH- 

1 

15 

6.2 

455 

960.2 

12.0 

60 

4.5 

/// 

//// 

— CD-FGH- 

1 

20 

5.7 

480 

957.4 

12.4 

58 

4.4 

342 

3.3 

— CD-FGH- 

1 

25 

5.7 

510 

953.9 

12.4 

57 

4.2 

/// 

//// 

--CD-FGH- 

1 

30 

5.9 

541 

950.4 

12.1 

58 

4.1 

343 

3.5 

— CD-FGH- 

1 

35 

6.1 

572 

946.9 

11.8 

58 

3.9 

/// 

//// 

— CD-FGH- 

1 

40 

5.8 

600 

943.7 

11.5 

59 

3.8 

343 

3.9 

— CD-FGH- 

1 

45 

5.7 

626 

940.8 

11.3 

60 

3.9 

/// 

//// 

--CD-FGH- 

1 

50 

5.7 

651 

938.0 

11.0 

61 

3.8 

343 

4.0 

--CD-FGH- 

1 

55 

5.7 

685 

934.1 

10.8 

62 

3.9 

/// 

//// 

— CD-FGH- 

2 

0 

5.7 

712 

931.2 

10.6 

62 

3.7 

34  4 

4.2 

— CD-FGH- 

2 

5 

5.7 

740 

928.0 

10.4 

62 

3.5 

/// 

//// 

— CD-FGH- 

2 

10 

5.7 

772 

924.4 

10.4 

60 

3.0 

345 

4.3 

— CD-FGH- 

2 

15 

5.9 

804 

920.9 

10.4 

59 

2.8 

/// 

//// 

— CD-FGH- 

2 

20 

6.0 

831 

917.8 

10.1 

59 

2.5 

346 

4.4 

— CD-FGH- 

2 

25 

6.1 

863 

914.4 

9.9 

59 

2.3 

/// 

//// 

— CD-FGH- 

2 

30 

6.1 

896 

910.7 

9.6 

59 

2.0 

348 

4.4 

— CD-FGH- 

2 

35 

5.6 

913 

908.8 

9.4 

60 

2.1 

/// 

//// 

— CD-FGH- 

2 

40 

5.6 

940 

905.8 

9.2 

60 

1.9 

350 

4.2 

— CD-FGH- 

2 

45 

5.3 

970 

902.6 

9.1 

60 

1.8 

/// 

//// 

— CD-FGH- 

2 

50 

5.3 

991 

900.3 

8.9 

60 

1.6 

352 

4  .3 

— CD-FGH- 
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Sounding  program  REV  7 . 62  using  Omega 

Ship      :PT_SUR 

Location  :  36.80  N  121.80  W       3m 

Rejected  Sigma  stations:  a,b,c 

Phase  fitting  length  is  250  s  from    0  min  to  120  min 

Sounding  :  5 

RS-number:  309140543 

No  PTU  editing 

No  wind  editing 

Started  at:       16  MAY  95  18:01  GMT 


e 

AscRat 

e  Hgt/I 

4SL  Pressure 

Temp   RH 

Dewp   Dir  Spe< 

2d 

WndStat 

n 

s 

m/s 

m 

hPa 

degC 

Q. 
O 

degC 

deg 

m/s 

0 

0 

0.0 

3 

1013. 

4 

12.1 

90 

10.5 

279 

2.2 

0 

5 

2.7 

17 

1011. 

7 

13.3 

70 

8.0 

/// 

//// 

-BCD — GH- 

0 

10 

1.8 

21 

1011. 

3 

13.0 

72 

8.1 

285 

2.3 

-BCD — GH- 

0 

15 

1.7 

31 

1010. 

1 

12.9 

74 

8.4 

/// 

//// 

-BCD — GH- 

0 

20 

2.0 

42 

1008. 

8 

12.8 

75 

8.5 

281 

2.1 

-BCD--GH- 

0 

25 

1.7 

46 

1008. 

3 

12.8 

76 

8.7 

/// 

//// 

-BCD-FGH- 

0 

30 

1.8 

56 

1007. 

1 

12.  6 

77 

8.7 

281 

2.3 

-BCD-FGH- 

0 

35 

1.8 

70 

1005. 

3 

12.5 

78 

8.8 

/// 

//// 

-BCD-FGH- 

0 

40 

1.8 

75 

1004. 

9 

12.5 

78 

8.8 

261 

2.3 

-BCD-FGH- 

0 

45 

1.8 

83 

1003. 

9 

12.5 

77 

8.6 

/// 

//// 

-BCD-FGH- 

0 

50 

1.7 

93 

1002. 

6 

12.5 

75 

8.2 

257 

2.3 

-BCD-FGH- 

0 

55 

2.0 

103 

1001. 

4 

12.5 

74 

8.0 

/// 

//// 

-BCD-FGH- 

1 

0 

1.9 

113 

1000. 

2 

12.5 

73 

7.8 

249 

2.4 

-BCD-FGH- 

1 

5 

1.6 

119 

999. 

5 

12.4 

73 

7.7 

/// 

//// 

-BCD-FGH- 

1 

10 

1.8 

128 

998. 

5 

12.4 

73 

7.7 

249 

2.5 

-BCD-FGH- 

1 

15 

1.9 

137 

997. 

4 

12.5 

73 

7.8 

/// 

//// 

-BCD-FGH- 

1 

20 

1.9 

150 

995. 

,8 

12.5 

73 

7.8 

249 

2.7 

-BCD-FGH- 

1 

25 

1.6 

153 

995. 

,5 

12.7 

72 

7.8 

/// 

//// 

-BCD-FGH- 

1 

30 

1.6 

160 

994. 

,6 

12.8 

70 

7.5 

249 

2.7 

-BCD-FGH- 

1 

35 

1.4 

164 

994. 

.1 

12.8 

70 

7.5 

/// 

//// 

-BCD-FGH- 

1 

40 

1.9 

185 

991. 

,7 

12.8 

69 

7.3 

252 

2.7 

-BCD-FGH- 

1 

45 

1.7 

190 

991, 

,1 

12.8 

69 

7.3 

/// 

//// 

-BCD-FGH- 

1 

50 

1.7 

201 

989. 

.8 

12.8 

69 

7.3 

250 

2.7 

-BCD-FGH- 

1 

55 

1.8 

208 

988, 

.9 

12.7 

69 

7.2 

/// 

//// 

-BCD-FGH- 

2 

0 

1.9 

218 

987, 

.8 

12.6 

69 

7.1 

247 

2.8 

-BCD-FGH- 

2 

5 

2.3 

228 

986, 

.6 

12.6 

69 

7.1 

/// 

//// 

-BCD-FGH- 

2 

10 

1.6 

234 

985 

.9 

12.5 

69 

7.0 

247 

2.9 

-BCD-FGH- 

2 

15 

1.8 

244 

984 

.7 

12.4 

69 

6.9 

/// 

//// 

-BCD-FGH- 

2 

20 

1.8 

256 

983 

.3 

12.3 

68 

6.6 

247 

2.9 

-BCD-FGH- 

2 

25 

2.0 

266 

982 

.1 

12.3 

69 

6.8 

/// 

//// 

-BCD-FGH- 

2 

30 

1.9 

275 

981 

.1 

12.1 

69 

6.6 

253 

2.9 

-BCD-FGH- 

2 

35 

1.7 

283 

980 

.1 

12.1 

69 

6.6 

/// 

//// 

-BCD-FGH- 

2 

40 

1.9 

291 

979 

.2 

12.1 

69 

6.6 

257 

3.0 

-BCD-FGH- 

2 

45 

1.9 

302 

977 

.9 

11.9 

70 

6.7 

/// 

//// 

-BCD-FGH- 

2 

50 

1.8 

311 

976 

.8 

11.8 

70 

6.6 

251 

3.1 

-BCD-FGH- 

2 

55 

1.8 

320 

975 

.7 

11.9 

71 

6.9 

/// 

//// 

-BCD-FGH- 
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Start  Up  Date  16  MAY  95  20:09  GMT 

System  test  passed   -No  errors  found 

Sounding  program  REV  7 . 62  using  Omega 

Ship      :PT_SUR 

Location  :  36.80  N  122.00  W       3m 

Rejected  Sigma  stations:  a,b,c 

Phase  fitting  length  is  250  s  from 

Sounding  :  6 

RS-number:  309140540 

No  PTU  editing 

No  wind  editing 

Started  at:       16  MAY  95  20:17  GMT 


0  min  to  120  min 


e 

Asc 

Rate  Hgt/IY 

ISL  Pressure 

Temp   RH 

Dewp   Dir  Spe< 

3d 

WndStat 

n 

s 

m/s 

m 

hPa 

degC 

O 

degC 

deg 

m/s 

0 

0 

0.0 

3 

1011. 

3 

12.1 

92 

10.8 

292 

10.4 

0 

5 

-1.7 

-7 

1012. 

5 

15.5 

59 

7.6 

/// 

//// 

— CD-FGH- 

0 

10 

-1.0 

-7 

1012. 

5 

15.6 

61 

8.2 

300 

10.5 

--CD-FGH- 

0 

15 

-0.6 

-7 

1012. 

5 

15.5 

63 

8.5 

/// 

//// 

— CD-FGH- 

0 

20 

1.5 

34 

1007. 

6 

14.9 

64 

8.2 

301 

10.7 

— CD-FGH- 

0 

25 

4.6 

119 

997. 

4 

13.6 

70 

8.3 

/// 

//// 

--CD-FGH- 

0 

30 

4.9 

149 

993. 

9 

12.9 

74 

8.4 

298 

10.2 

— CD-FGH- 

0 

35 

6.4 

178 

990. 

5 

12.7 

75 

8.4 

/// 

//// 

— CD — GH- 

0 

40 

6.9 

200 

987. 

8 

12.4 

76 

8.3 

299 

9.9 

--CD--GH- 

0 

45 

8.0 

228 

984. 

6 

12.6 

69 

7.1 

/// 

//// 

--CD-FGH- 

0 

50 

7.2 

251 

981. 

8 

13.2 

58 

5.2 

298 

9.7 

--CD-FGH- 

0 

55 

4.9 

265 

980. 

1 

•  13.4 

54 

4.3 

/// 

//// 

-BCD-FGH- 

1 

0 

4.6 

286 

977. 

7 

13.7 

50 

3.5 

299 

9.7 

-BCD-FGH- 

1 

5 

4  .4 

312 

974. 

,8 

13.9 

46 

2.5 

/// 

//// 

-BCD-FGH- 

1 

10 

4.3 

329 

972. 

,7 

14.1 

43 

1.8 

297 

9.1 

-BCD-FGH- 

1 

15 

3.6 

337 

971. 

,8 

14.2 

42 

1.5 

/// 

//// 

-BCD-FGH- 

1 

20 

3.4 

352 

970. 

,1 

14.3 

41 

1.3 

298 

9.1 

-BCD-FGH- 

1 

25 

3.3 

366 

968. 

.5 

14.4- 

41 

1.4 

/// 

//// 

-BCD-FGH- 

1 

30 

3.  1 

380 

966. 

.8 

14.4 

41 

1.4 

299 

9.0 

-BCD-FGH- 

1 

35 

2.8 

396 

965. 

.1 

14.5 

40 

1.1 

/// 

//// 

-BCD-FGH- 

1 

40 

2.7 

409 

963. 

.5 

14.6 

40 

1.2 

300 

9.1 

-BCD-FGH- 

1 

45 

2.8 

421 

962, 

.1 

14.5 

40 

1.1 

/// 

//// 

-BCD-FGH- 

1 

50 

2.9 

440 

960. 

.0 

14.5 

40 

1.1 

298 

8.7. 

-BCD-FGH- 

1 

55 

3.0 

451 

958, 

.7 

14.4 

39 

0.7 

/// 

//// 

-BCD-FGH- 

2 

0 

2.8 

465 

957, 

.2 

14.3 

39 

0.6 

298 

8.8 

-BCD-FGH- 

2 

5 

2.8 

482 

955 

.2 

14.3 

38 

0.2 

/// 

//// 

-BCD-FGH- 

2 

10 

3.0 

498 

953 

.5 

14.2 

38 

0.2 

299 

9.0 

-BCD-FGH- 

2 

15 

3.1 

514 

951 

.6 

14.2 

38 

0.2 

/// 

//// 

-BCD-FGH- 

2 

20 

2.8 

524 

950 

.4 

14.1 

37 

-0.3 

300 

8.9 

-BCD-FGH- 

2 

25 

2.8 

539 

948 

.8 

14.1 

37 

-0.3 

/// 

//// 

-BCD-FGH- 

2 

30 

2.9 

551 

947 

.4 

14.0 

37 

-0.4 

299 

8.9 

-BCD-FGH- 

2 

35 

2.8 

565 

945 

.8 

14.0 

37 

-0.4 

/// 

//// 

-BCD-FGH- 
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Sounding  program  REV  7 . 62  using  Loran-C 

Ship      : PT_SUR 

Location  :  36.80  N  121.80  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from    0  min  to   10  min 

Phase  fitting  length  is  120  s  from   10  min  to   45  min 

Phase  fitting  length  is  240  s  from   45  min  to  120  min 

Sounding  :  10 

RS-number:  108234745 

No  PTU  editing 

No  wind  editing 

Started  at:       17  MAY  95  17:47  GMT 

Time  AscRate  Hgt/MSL  Pressure    Temp   RH    Dewp   Dir  Speed      WndStat 
min   s     m/s       m      hPa   degC    %    degC   deg   m/s 


0 

0 

0. 

,0 

3 

1016. 

8 

12. 

4 

88 

10. 

5 

278 

1.0 

000000000 

0 

5 

4. 

2 

24 

1014. 

1 

12. 

7 

63 

5. 

9 

/// 

//// 

111100000 

0 

10 

4. 

1 

44 

1011. 

7 

12. 

4 

66 

6. 

3 

268 

2.2 

111100000 

0 

15 

3. 

8 

60 

1009. 

8 

12. 

3 

69 

6. 

8 

/// 

//// 

111100000 

0 

20 

3. 

5 

74 

1008. 

0 

12. 

1 

70 

6. 

8 

263 

2.5 

111100000 

0 

25 

3. 

5 

91 

1006. 

1 

12. 

0 

72 

7. 

2 

/// 

//// 

111100000 

0 

30 

3. 

5 

107 

1004. 

1 

11. 

8 

72 

7. 

0 

263 

2.3 

111100000 

0 

35 

3. 

4 

125 

1001. 

9 

11. 

6 

73 

7. 

0 

/// 

//// 

111100000 

0 

40 

3. 

,3 

144 

999. 

7 

11. 

4 

74 

7, 

0 

259 

2.4 

111100000 

0 

45 

3. 

,3 

158 

998. 

0 

11. 

3 

75 

7. 

1 

/// 

//// 

111100000 

0 

50 

3. 

,3 

174 

996. 

1 

11. 

1 

76 

7. 

1 

256 

2.3 

111100000 

0 

55 

3. 

,7 

198 

993. 

1 

10. 

,9 

77 

7. 

,1 

/// 

//// 

111100000 

1 

0 

3. 

,5 

213 

991. 

5 

10. 

7 

78 

7. 

,1 

254 

2.3 

111100000 

1 

5 

3. 

,6 

234 

988. 

9 

10. 

5 

79 

7, 

,1 

/// 

//// 

111100000 

1 

10 

3. 

,7 

255 

986. 

4 

10. 

3 

80 

7. 

,1 

250 

2.3 

111100000 

1 

15 

4, 

.0 

280 

983. 

5 

10. 

,0 

82 

7. 

,1 

/// 

//// 

111100000 

1 

20 

4. 

.3 

302 

980. 

9 

9. 

,8 

83 

7. 

,1 

251 

2.3 

111100000 

1 

25 

3. 

.9 

318 

978. 

,9 

9. 

,7 

85 

7. 

,3 

/// 

//// 

111100000 

1 

30 

4  , 

.0 

332 

977. 

,3 

9. 

,5 

85 

7, 

,2 

250 

2.2 

111100000 

1 

35 

3. 

,8 

346 

975. 

,6 

9. 

,4 

86 

7, 

,2 

/// 

//// 

111100000 

1 

40 

4, 

.0 

375 

972. 

,3 

9. 

,1 

88 

7, 

.3 

246 

1.8 

111100000 

1 

45 

4  , 

.0 

399 

969. 

.4 

8. 

.9 

89 

7 

.2 

/// 

//// 

111100000 

1 

50 

3, 

.8 

417 

967, 

,3 

8, 

.6 

89 

6 

.9 

247 

1.3 

111100000 

1 

55 

4 

.0 

434 

965, 

.3 

8, 

.6 

90 

7, 

.1 

/// 

//// 

111100000 

2 

0 

4 

.2 

457 

962, 

,6 

8, 

.4 

90 

6 

.9 

241 

1.7 

111100000 

2 

5 

4 

.6 

483 

959, 

.5 

8, 

.2 

90 

6 

.7 

/// 

//// 

111100000 

2 

10 

4 

.5 

509 

956, 

.5 

8 

.1 

91 

6 

.8 

237 

2.1 

111100000 

2 

15 

4 

.6 

537 

953, 

.4 

8 

.1 

88 

6 

.3 

/// 

//// 

111100000 

2 

20 

4 

.4 

548 

952 

.1 

8 

.7 

79 

5 

.3 

240 

2.9 

111100000 

2 

25 

4 

.2 

562 

950 

.4 

10 

.2 

70 

5 

.0 

/// 

//// 

111100000 

2 

30 

4 

.0 

576 

948 

.8 

11 

.8 

63 

5 

.0 

240 

3.4 

111100000 

2 

35 

3 

.4 

590 

947 

.2 

12 

.6 

59 

4 

.9 

/// 

//// 

111100000 

2 

40 

3 

.2 

605 

945 

.6 

12 

.7 

55 

4 

.0 

237 

4.0 

111100000 

2 

45 

2 

.8 

620 

943 

.9 

12 

.8 

52 

3 

.3 

/// 

//// 

111100000 

2 

50 

2 

.9 

635 

942 

.1 

13 

.1 

50 

3 

.0 

236 

4.2 

111100000 

118 


Sounding  program  REV  7 . 62  using  Loran-C 

Ship      : PT_SUR 

Location  :  36.80  N  122.00  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from    0  min  to   10  min 

Phase  fitting  length  is  120  s  from   10  min  to   45  min 

Phase  fitting  length  is  240  s  from   45  min  to  120  min 

Sounding  :  11 

RS-number:  108234743 

No  PTU  editing 

No  wind  editing 

Started  at:       17  MAY  95  21:54  GMT 

Time  AscRate  Hgt/MSL  Pressure    Temp   RH    Dewp   Dir  Speed     WndStat 
min   s      m/s        m      hPa    degC    %    degC   deg    m/s 


0 

0 

0. 

0 

3 

1016. 

3 

12. 

2 

96 

11. 

6 

261 

5.4 

000000000 

0 

5 

3. 

8 

20 

1014. 

1 

13. 

3 

72 

8. 

4 

/// 

//// 

111100000 

0 

10 

4. 

1 

44 

1011. 

3 

12. 

6 

76 

8. 

5 

309 

5.3 

111100000 

0 

15 

4. 

2 

67 

1008. 

4 

12. 

4 

78 

8. 

7 

/// 

//// 

111100000 

0 

20 

4. 

1 

85 

1006. 

,3 

12. 

2 

79 

8. 

.1 

302 

5.7 

111100000 

0 

25 

3. 

8 

100 

1004. 

5 

12. 

1 

81 

9. 

,0 

/// 

//// 

111100000 

0 

30 

3. 

7 

115 

1002. 

,6 

11. 

,9 

83 

9. 

2 

300 

5.8 

111100000 

0 

35 

3. 

7 

129 

1000. 

9 

11. 

8 

83 

9. 

,1 

/// 

//// 

111100000 

0 

40 

3. 

5 

148 

998. 

7 

11. 

6 

84 

9. 

,0 

286 

7.0 

111100000 

0 

45 

3. 

2 

160 

997. 

,3 

11. 

5 

84 

8. 

,9 

/// 

//// 

110100000 

0 

50 

3. 

3 

185 

994. 

,4 

11. 

3 

85 

8. 

,9 

272 

9.7 

110100000 

0 

55 

3. 

3 

197 

992. 

,9 

11. 

,2 

86 

9. 

,0 

/// 

//// 

111100000 

1 

0 

3. 

3 

215 

990. 

,7 

11. 

,0 

87 

9. 

,0 

281 

8.0 

111100000 

1 

5 

3. 

3 

229 

989. 

,0 

10. 

,9 

86 

8. 

,7 

/// 

//// 

111100000 

1 

10 

3. 

4 

250 

986. 

,6 

10. 

,7 

87 

8. 

.7 

291 

7.1 

111100000 

1 

15 

3. 

,5 

265 

984, 

,8 

10. 

,7 

87 

8. 

.7 

/// 

//// 

111100000 

1 

20 

3. 

,3 

284 

982. 

,5 

10. 

,5 

87 

8. 

.5 

291 

7.3 

111100000 

1 

25 

3. 

,4 

296 

981. 

.1 

10, 

.4 

87 

8, 

.4 

/// 

//// 

111100000 

1 

30 

3. 

.2 

312 

979. 

.2 

10, 

,2 

87 

8, 

.2 

291 

7.2 

111100000 

1 

35 

3. 

,5 

336 

976. 

,4 

10. 

,1 

88 

8, 

.2 

/// 

//// 

111100000 

1 

40 

3. 

,4 

351 

974, 

,7 

9, 

,9 

88 

8, 

.1 

292 

7.1 

111100000 

1 

45 

3. 

,4 

368 

972, 

,6 

9. 

.7 

89 

8, 

.0 

/// 

//// 

111100000 

1 

50 

3, 

,2 

379 

971, 

,3 

9, 

.6 

89 

7, 

.9 

290 

7.2 

111100000 

1 

55 

3. 

.4 

401 

968, 

,7 

9, 

,6 

90 

8 

.1 

/// 

//// 

111100000 

2 

0 

3, 

,5 

416 

967, 

.1 

9, 

.4 

90 

7, 

.9 

288 

7.0 

111100000 

2 

5 

3. 

,4 

439 

964, 

.4 

9, 

.3 

91 

8 

.0 

/// 

//// 

111100000 

2 

10 

3, 

.4 

452 

962, 

.8 

9, 

.2 

91 

7, 

.  9 

287 

6.9 

111100000 

2 

15 

3, 

.3 

468 

960 

.9 

9 

.1 

92 

7 

.9 

/// 

//// 

111100000 

2 

20 

3, 

.7 

4  90 

958 

.4 

8 

.9 

91 

7 

.6 

285 

6.5 

111100000 

2 

25 

3, 

.4 

505 

956 

.6 

8 

.9 

91 

7 

.6 

/// 

//// 

111100000 

2 

30 

3, 

.5 

522 

954 

.6 

9 

.4 

89 

7 

.7 

285 

5.6 

111100000 

2 

35 

3. 

.3 

538 

952 

.9 

9 

.5 

88 

7 

.7 

/// 

//// 

111100000 

2 

40 

3 

.4 

555 

950 

.9 

9 

.8 

83 

7 

.1 

283 

4.9 

111100000 

2 

45 

3 

.4 

571 

949 

.1 

10 

.2 

78 

6 

.6 

/// 

//// 

111100000 

119 


3  m 


Sounding  program  REV  7 . 62  using  Loran-C 

Ship      :PT_SUR 

Location  :  36.70  N  121.90  W 

Loran-C  chain  1:   9940 

Phase  fitting  length  is 

Phase  fitting  length  is 

Phase  fitting  length  is 

Sounding  :  14 

RS-number:  108234746 

No  PTU  editing 

No  wind  editing 

Started  at:       18  MAY  95  15:17  GMT 


60 

s  from 

0 

min 

to 

10 

min 

120 

s  from 

10 

min 

to 

45 

min 

240 

s  from 

45 

min 

to 

120 

min 

Time  AscRate  Hgt/MSL  Pressure    Temp   RH 
min   s     m/s       m      hPa    degC 


Dewp   Dir  Speed 
i    degC   deg    m/s 


WndStat 


0 

0 

0. 

0 

3 

1016. 

6 

10. 

4 

93 

9. 

3 

237 

0.6 

000000000 

0 

5 

4. 

0 

25 

1014. 

0 

11. 

6 

72 

6. 

8 

/// 

//// 

111100000 

0 

10 

3. 

5 

38 

1012. 

5 

11. 

0 

75 

6. 

8 

240 

3.1 

111100000 

0 

15 

3. 

8 

57 

1010. 

1 

10. 

9 

78 

7. 

3 

/// 

//// 

111100000 

0 

20 

3. 

5 

72 

1008. 

3 

10. 

8 

79 

7. 

4 

241 

2.6 

111100000 

0 

25 

3. 

3 

85 

1006. 

7 

10. 

7 

81 

7. 

6 

/// 

//// 

111100000 

0 

30 

3. 

3 

103 

1004. 

,6 

10. 

,5 

82 

7. 

6 

237 

1.0 

111100000 

0 

35 

3. 

3 

123 

1002. 

2 

10. 

,3 

83 

7. 

6 

/// 

//// 

111100000 

0 

40 

3. 

5 

143 

999. 

,1 

10, 

,1 

84 

7. 

6 

254 

0.5 

111100000 

0 

45 

3. 

5 

163 

997. 

,3 

10, 

,0 

85 

7. 

6 

/// 

//// 

111100000 

0 

50 

3. 

7 

182 

995, 

,1 

9, 

,7 

86 

7. 

5 

243 

2.1 

111100000 

0 

55 

3. 

8 

198 

993. 

,1  ■ 

9, 

,6 

88 

7. 

8 

/// 

//// 

110100000 

1 

0 

3. 

7 

214 

991, 

,2 

9, 

,4 

88 

7. 

6 

246 

2.5 

110100000 

1 

5 

3. 

7 

234 

988. 

,8 

9, 

,3 

89 

7. 

6 

/// 

//// 

111100000 

1 

10 

3. 

8 

256 

986. 

,1 

9, 

,1 

90 

7, 

,6 

254 

1.8 

111100000 

1 

15 

3. 

.6 

274 

984. 

,0 

8, 

,9 

91 

7. 

,6 

/// 

//// 

111100000 

1 

20 

3. 

,6 

291 

982. 

.1 

8, 

,8 

92 

7, 

,6 

274 

0.6 

111100000 

1 

25 

3. 

.6 

304 

980, 

,5 

8, 

,  7- 

92 

7. 

,5 

/// 

//// 

111100000 

1 

30 

3. 

,6 

321 

978, 

.5 

8, 

.6 

92 

7, 

,4 

313 

0.8 

111100000 

1 

35 

3, 

,6 

341 

976, 

,1 

8, 

.5 

93 

7, 

,5 

/// 

//// 

111100000 

1 

40 

3, 

,7 

367 

973, 

.0 

8 

.4 

93 

7, 

,4 

326 

0.8 

111100000 

1 

45 

3, 

.8 

385 

970 

.9 

8 

.3 

93 

7, 

.3 

/// 

//// 

111100000 

1 

50 

3, 

.8 

405 

968 

.5 

8 

.1 

93 

7, 

.1 

334 

0.9. 

111100000 

1 

55 

3, 

.9 

423 

966 

.4 

8 

.1 

94 

7, 

.2 

/// 

//// 

111100000 

2 

0 

3, 

.9 

437 

964 

.7 

8 

.0 

94 

7 

,1 

347 

1.0 

111100000 

2 

5 

3. 

.7 

454 

962 

.7 

7 

.9 

94 

7, 

.0 

/// 

//// 

111100000 

2 

10 

3 

.5 

472 

960 

.7 

7 

.8 

94 

6 

,9 

338 

1.0 

111100000 

2 

15 

3 

.6 

493 

958 

.3 

7 

.7 

94 

6 

.8 

/// 

//// 

111100000 

2 

20 

3 

.5 

510 

956 

.3 

7 

.6 

94 

6 

.7 

336 

0.9 

111100000 

2 

25 

3 

.5 

528 

954 

.2 

7 

.5 

94 

6 

.6 

/// 

//// 

111100000 

2 

30 

3 

.7 

548 

951 

.8 

7 

.4 

94 

6 

.5 

338 

1.0 

111100000 

2 

35 

3 

.9 

567 

949 

.6 

7 

.3 

94 

6 

.4 

/// 

//// 

111100000 

2 

40 

3 

.7 

582 

947 

.9 

7 

.2 

94 

6 

.3 

348 

1.1 

111100000 

2 

45 

3 

.7 

604 

945 

.4 

7 

.0 

94 

6 

.1 

/// 

//// 

111100000 

120 


Sounding  program  REV  7 . 62  using  Omega 

Ship      :PT_SUR 

Location  :  36.80  N  121.90  W       3m 

Rejected  Sigma  stations:  a,b,c 

Phase  fitting  length  is  250  s  from 

Sounding  :  16 

RS-number:  158735752 

No  PTU  editing 

No  wind  editing 

Started  at:       18  MAY  95  21:53  GMT 


0  min  to  120  min 


e 

AscRat 

e  Hqt/K 

ISL  Pressure 

Temp   RH 

Dewp   Dir  Spei 

3d 

WndStat 

n 

s 

m/s 

m 

hPa 

degC 

O 

degC 

deg 

m/s 

0 

0 

0.0 

3 

1015. 

5 

11.2 

90 

9.6 

255 

6.3 

0 

5 

3.7 

21 

1013. 

2 

12.4 

63 

5.6 

/// 

//// 

— CD — GH 

0 

10 

4.5 

48 

1010. 

0 

11.8 

67 

5.9 

287 

4.3 

--CD--GH 

0 

15 

5.4 

77 

1006. 

5 

11.6 

70 

6.4 

/// 

//// 

--CD--GH 

0 

20 

5.2 

107 

1002. 

9 

11.2 

72 

6.4 

288 

3.9 

--CD--GH 

0 

25 

5.3 

137 

999. 

2 

11.0 

74 

6.6 

/// 

//// 

--CD--GH 

0 

30 

5.5 

168 

995. 

6 

10.7 

76 

6.7 

283 

3.5 

--CD--GH 

0 

35 

5.8 

190 

992. 

9 

10.4 

77 

6.6 

/// 

//// 

--CD--GH 

0 

40 

5.6 

216 

989. 

8 

10.2 

78 

6.6 

283 

3.3 

--CD--GH 

0 

45 

5.5 

252 

985. 

5 

9.8 

80 

6.6 

/// 

//// 

--CD--GH 

0 

50 

5.7 

279 

982. 

3 

9.6 

81 

6.5 

277 

3.0 

— CD--GH 

0 

55 

5.5 

301 

979. 

7 

9.3 

82 

6.4 

/// 

//// 

--CD-FGH 

1 

0 

5.3 

328 

976. 

,6 

9.2 

83 

6.5 

277 

2.8 

--CD-FGH 

1 

5 

5.1 

346 

974. 

,4 

9.3 

83 

6.6 

/// 

//// 

--CD-FGH 

1 

10 

5.0 

366 

972. 

,0 

9.4 

83 

6.7 

264 

2.3 

--CD-FGH 

1 

15 

4.4 

383 

970. 

,0 

9.7 

81 

6.6 

/// 

//// 

--CD-FGH 

1 

20 

4.0 

398 

968. 

,2 

10.1 

78 

6.5 

258 

2.3 

— CD-FGH 

1 

25 

3.9 

418 

966. 

,0 

12.3 

63 

5.5 

/// 

//// 

--CD-FGH 

1 

30 

3.7 

439 

963. 

,5 

13.0 

59 

5.2 

258 

2.3 

--CD-FGH 

1 

35 

3.9 

460 

961. 

,1 

13.2 

59 

5.4 

/// 

//// 

— CD-FGH 

1 

40 

3.7 

478 

959. 

,1 

13.4 

51 

3.5 

260 

2.3 

--CD-FGH 

1 

45 

3.7 

494 

957. 

,2 

15.1 

33 

-1.0 

/// 

//// 

--CD-FGH 

1 

50 

3.8 

513 

955. 

,1 

15.9 

38 

1.7 

261 

2.3 

—CD-FGH 

1 

55 

3.5 

525 

953. 

.7 

16.0 

39 

2.1 

/// 

//// 

-BCD-FGH 

2 

0 

3.2 

536 

952, 

.5 

16.0 

39 

2.1 

258 

2.3 

-BCD-FGH 

2 

5 

3.1 

559 

949, 

.9 

15.9 

40 

2.4 

/// 

//// 

-BCD-FGH 

2 

10 

3.1 

571 

948, 

.6 

15.8 

41 

2.6 

260 

2.2 

-BCD-FGH 

2 

15 

3.3 

592 

946, 

.2 

15.8 

43 

3.3 

/// 

//// 

-BCD-FGH 

2 

20 

3.1 

606 

944 

.7 

15.8 

44 

3.6 

254 

2.1 

-BCD-FGH 

2 

25 

3.3 

625 

942 

.5 

15.9 

44 

3.7 

/// 

//// 

-BCD-FGH 

2 

30 

3.6 

643 

940 

.5 

15.9 

44 

3.7 

261 

2.0 

-BCD-FGH 

2 

35 

3.5 

658 

938 

.8 

16.1 

41 

2.9 

/// 

//// 

-BCD-FGH 

2 

40 

3.4 

674 

937 

.1 

16.3 

37 

1.7 

266 

1.8 

-BCD-FGH 

2 

45 

3.3 

693 

935 

.0 

16.3 

38 

2.0 

/// 

//// 

-BCD-FGH 

2 

50 

3.5 

711 

933 

.0 

16.2 

39 

2.3 

267 

1.8 

-BCD-FGH 

2 

55 

3.5 

730 

930 

.9 

16.1 

41 

2.9 

/// 

//// 

-BCD-FGH 

121 


Sounding  program  REV  7 . 62  using  Omega 

Ship      :PT_SUR 

Location  :  36.78  N  121.90  W  3m 

Rejected  Sigma  stations:  a,b,c 

Phase  fitting  length  is  250  s  from    0  min  to  120  min 

Sounding  :  18 

RS-number:  158735842 

No  PTU  editing 

No  wind  editing 

Started  at:       19  MAY  95  15:35  GMT 

Time  AscRate  Hgt/MSL  Pressure  Temp   RH    Dewp   Dir  Speed  WndStat 

min   s      m/s  m  hPa  degC    %  degC  deg  m/s 

0   0      0.0  3  1014.6  10.0  100  10.0  93    0.9  

0   5      3.3  22  1012.4  11.0  80  7.7  ///  ////  — CD-FGH 

0  10      3.9  42  1010.0  10.7  81  7.6  9  1.4  — CD-FGH 

0  15      3.4     •   55  1008.4  10.7  83  8.0  ///  ////  --CD-FGH 

0  20      3.3  70  1006.6  10.5  84  8.0  6  1.6  --CD-FGH 

0  25      3.3  87  1004.5  10.4  86  8.2  ///  ////  --CD-FGH 

0  30      3.3  102  1002.6  10.3  86  8.1  1  2.2  --CD-FGH 

0  35      3.4  122  1000.3  10.2  87  8.2  ///  ////  —CD-FGH 

0  40      3.3  141  998.0  10.1  87  8.1  5  2.1  —CD-FGH 

0  45      3.4  156  996.2  10.0  88  8.1  ///  ////  —CD-FGH 

0  50      3.4  173  994.1  9.9  88  8.1  15  1.9  --CD-FGH 

0  55      3.4  189  992.2  9.8  88  8.0  ///  ////  --CD-FGH 

1  0  3.4  204  990.5  9.7  88  7.9  23  1.9  --CD-FGH 
1  5  3.3  220  988.6  9.7  89  8.0  ///  ////  -BCD-FGH 
1  10  3.2  238  986.4  9.6  89  7.9  23  2.0  -BCD-FGH 
1  15  3.3  255  984.3  9.5  89  7.8  ///  ////  -BCD-FGH 
1  20  3.2  270  982.5  9.4  89  7.7  24  2.2  -BCD-FGH 
1  25  3.3  289  980.3  9.3  89  7.6  ///  ////  -BCD-FGH 
1  30  3.4  305  978.5  9.2  89  7.5  29  2.2  -BCD-FGH 
1  35  3.4  322  976.4  9.1  89  7.4  ///  ////  -BCD-FGH 
1  40  3.4  341  974.2  9.0  89  7.3  33  2.2  -BCD-FGH 
1  45  3.3  354  972.6  9.0  89  7.3  ///  ////  -BCD-FGH 
1  50      3.3  369  970.9  9.2  89  7.5  31  2.2  -BCD-FGH 

1  55      3.2  382  969.3  10.1  90  8.6  ///  ////  -BCD-FGH 

2  0  3.0  395  967.8  10.6  90  9.1  30  2.2  -BCD-FGH 
2  5  2.9  409  966.2  11.4  90  9.9  ///  ////  -BCD-FGH 
2  10  2.6  418  965.2  11.6  88  9.7  28  2.3  -BCD-FGH 
2  15  2.5  429  963.9  12.1  88  10.2  ///  ////  -BCD-FGH 
2  20  2.3  438  962.8  12.3  86  10.1  31  2.1  -BCD-FGH 
2  25  2.5  454  960.9  12.3  82  9.4  ///  ////  -BCD-FGH 
2  30  2.4  467  959.5  12.4  78  8.7  31  2.2  -BCD-FGH 
2  35  2.4  481  957.9  12.5  78  8.8  ///  ////  -BCD-FGH 
2  40  2.5  493  956.5  12.9  79  9.4  38  2.1  -BCD-FGH 
2  45  2.6  506  955.0  13.2  80  9.9  ///'  ////  -BCD-FGH 
2  50  2.7  518  953.7  13.3  81  10.2  44  2.1  -BCD-FGH 
2  55      2.6  533  952.0  13.7  77  9.8  ///  ////  -BCD-FGH 
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Sounding   program   REV   7 . 62    using    Loran-C 

Ship  :PT_SUR 

Location  :  36.80  N  121.80  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  19 

RS-number:  108234748 

No  PTU  editing 

No  wind  editing 

Started  at:       19  MAY  95  18:27  GMT 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure    Temp   RH 
min   s      m/s        m      hPa    degC 


Dewp   Dir  Speed 
i    degC   deg    m/s 


WndStat 


0 

0 

0. 

,0 

3 

1013. 

4 

10. 

5 

96 

9. 

9 

263 

1.4 

000000000 

0 

5 

2. 

7 

13 

1012. 

2 

11. 

3 

78 

7. 

7 

/// 

//// 

110100000 

0 

10 

2. 

6 

29 

1010. 

3 

11. 

0 

82 

8. 

1 

299 

1.7 

110100000 

0 

15 

3. 

2 

51 

1007. 

7 

10. 

9 

85 

8. 

5 

/// 

//// 

111100000 

0 

20 

3. 

5 

72 

1005. 

1 

10. 

7 

88 

8. 

8 

281 

2.0 

111100000 

0 

25 

3. 

.7 

95 

1002. 

,3 

10. 

6 

89 

8. 

,9 

/// 

//// 

111100000 

0 

30 

3. 

7 

115 

1000. 

,0 

10. 

5 

90 

9. 

0 

271 

2.2 

111100000 

0 

35 

3. 

,9 

136 

997. 

4 

10. 

3 

90 

8. 

8 

/// 

//// 

111100000 

0 

4  0 

3. 

,9 

147 

996. 

,1 

10. 

,3 

90 

8. 

8 

318 

1.3 

111100000 

0 

45 

3. 

,7 

161 

994, 

,4 

10. 

,3 

90 

8, 

,8 

/// 

//// 

111100000 

0 

50 

3. 

.6 

179 

992. 

,2 

10. 

1 

90 

8, 

,6 

250 

3.2 

111100000 

0 

55 

3. 

,4 

200 

989, 

,8 

10. 

,1 

90 

8, 

,6 

/// 

//// 

111100000 

1 

0 

3. 

,4 

216 

987, 

,8 

10. 

0 

90 

8, 

,5 

270 

1.4 

111100000 

1 

5 

3, 

,3 

234 

985, 

,7 

9. 

,9 

90 

8. 

,4 

/// 

//// 

111100000 

1 

10 

3, 

,2 

244 

984, 

.5 

9. 

,8 

90 

8, 

,3 

12 

1.2 

111100000 

1 

15 

3, 

.3 

257 

982, 

.9 

9. 

,8 

90 

8, 

,3 

/// 

//// 

111100000 

1 

20 

3, 

,4 

281 

980, 

,1 

9, 

,1 

89 

8, 

.0 

14 

1.1 

111100000 

1 

25 

3. 

.4 

300 

977, 

,9 

9, 

,7 

89 

8, 

.0 

/// 

//// 

111100000 

1 

30 

3, 

,0 

305 

977, 

,3 

9, 

,6 

89 

7, 

,9 

358 

0.6 

111100000 

1 

35 

3, 

,4 

337 

973, 

,5 

9, 

,5 

89 

7, 

,8 

/// 

//// 

111100000 

1 

40 

3. 

.6 

353 

971, 

,6 

9, 

,4 

89 

7, 

,7 

0 

0.0 

111100000 

1 

45 

3, 

.5 

369 

969, 

,8 

9, 

,3 

89 

7, 

.6 

/// 

//// 

111100000 

1 

50 

3, 

.5 

386 

967, 

.8 

9, 

,1 

88 

7, 

.3 

4 

0.5 

111100000 

1 

55 

3, 

.3 

398 

966, 

.4 

10, 

,2 

91 

8 

.8 

/// 

//// 

111100000 

2 

0 

3, 

.4 

406 

965, 

.4 

11, 

,6 

89 

9 

.9 

31 

1.0 

111100000 

2 

5 

2 

.6 

419 

963 

.9 

12, 

,6 

86 

10, 

.4 

/// 

//// 

111100000 

2 

10 

3 

.0 

443 

961 

.2 

13, 

.3 

77 

9 

.4 

44 

1.1 

111100000 

2 

15 

2 

.6 

449 

960 

.5 

13 

.3 

75 

9 

.0 

/// 

//// 

111100000 

2 

20 

2 

.1 

449 

960 

.5 

13, 

.2 

72 

8 

.3 

52 

1.5 

111100000 

2 

25 

3 

.1 

490 

955 

.8 

13, 

.6 

69 

8 

.1 

/// 

//// 

111100000 

2 

30 

3 

.1 

498 

954 

.9 

13 

.7 

70 

8 

.4 

66 

1.2 

111100000 

2 

35 

3 

.4 

515 

952 

.9 

14 

.1 

71 

9 

.0 

/// 

//// 

111100000 

2 

40 

2 

.8 

527 

951 

.6 

14 

.6 

66 

8 

.4 

107 

1.2 

111100000 

2 

45 

3 

.2 

544 

949 

.6 

15 

.1 

64 

8 

.4 

/// 

//// 

111100000 

123 


60 

s  from 

120 

s  from 

240 

s  from 

Sounding  program  REV  7 . 62  using  Loran-C 

Ship      : PT_SUR 

Location  :  36.80  N  122.00  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is 

Phase  fitting  length  is 

Phase  fitting  length  is 

Sounding  :  20 

RS-number:  108234840 

No  PTU  editing 

No  wind  editing 

Sounding  program  REV  7 . 62  using  Loran-C 

Ship      : PT_SUR 

Location  :  36.80  N  122.00  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  20 

RS-number:  108234753 

No  PTU  editing 

No  wind  editing 

Started  at:       19  MAY  95  21:31  GMT 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


0  min  to 
10  min  to 


10  min 
45  min 


45  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure    Temp   RH 
WndStat 
min   s      m/s       m      hPa    degC    % 


Dewp   Dir  Speed 
degC   deg    m/s 


0  0 
0  5 
0  10 
0  15 
0  20 
0  25 
0  30 
0  35 
0  40 
0  45 
0  50 
0  55 


0.0 

1.3 


0 
5 
10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
0 
5 
10 
15 
20 


3.3 

3.9 

3.8 

3. 

3. 

3. 

3. 

2. 

2. 

2. 

2. 

2, 

2, 

2, 

2 

2 

2 

2 

2 

2 

2 

2 

2 


30 
53 
72 
84 
103 
125 
143 
157 
170 
182 
192 
204 
216 
229 
241 
251 
259 
274 
285 
298 
308 
318 
328 
340 
349 
363 
374 


1012.3 

1011.5 

1009.0 

1006.2 

1003 

1002 

1000 

997 

995 

993 

992 


990.7 

989.5 
988.1 
986.6 
985.1 
983.7 
982.5 
981.6 
979.8 
978.5 
977.0 
975.8 
974.7 
973.5 
972.0 
971.1 
969.4 
968.2 


11.0 

12.4 

11.8 

11 

11 

11 


6 

4 
2 

11.0 
10.8 


10 
10 
10 
10 
10 
10 
10 
10 
10 
11 
11 
11 
11 
11 


14 
15 

15 


12.0 
13.2 


16.0 
16.4 


95  10.2  275  6.1  000000000 

71  7.3  ///  ////  111100000 

74  7.4  292  5.8  111100000 

75  7.4  ///  ////  111100000 

76  7.4  289  6.3  111100000 

78  7.6  ///  ////  111100000 

79  7.5  294  5.7  111100000 
81  7.7  ///  ////  111100000 
81  7.5  293  5.7  111100000 

81  7.4  ///  ////  111100000 

82  7.5  296  5.3  111100000 

83  7.6  ///  ////  111100000 
82  7.4  298  4.8  111100000 
82  7.3  ///  ////  111100000 
81  7.2  298  4.4  111100000 
79  7.0  ///  ////  111100000 
79  7.5  292  4.2  111100000 
79  7.7  ///  ////  111100000 
79  7.9  292  3.8  111100000 
79  8.0  ///  ////  111100000 
79  8.0  285  3.7  111100000 
79  8.1  ///  ////  111100000 

77  8.1  279  3.5  111100000 
59  5.4  ///  ////  111100000 
43  2.2  276  3.4  111100000 
43  3.0  ///  ////  111100000 
57  7.5  273  3.2  111100000 
57  7.6  ///  ////  111100000 
52  6.6  275  2.8  111100000 
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2 

25 

2. 

.2 

385 

966, 

,9 

16. 

,9 

53 

7. 

3 

/// 

//// 

111100000 

2 

30 

2, 

.3 

396 

965. 

,7 

17, 

,7 

45 

5. 

7 

277 

2.4 

111100000 

2 

35 

2, 

.2 

407 

964. 

,5 

18. 

1 

47 

6. 

7 

/// 

//// 

111100000 

2 

40 

2, 

.3 

417 

963. 

,3 

18. 

,1 

47 

6. 

7 

282 

2.0 

111100000 

2 

45 

2, 

.2 

428 

962. 

,  1 

18. 

0 

48 

6. 

9 

/// 

//// 

111100000 

2 

50 

2. 

,2 

440 

960. 

,7 

18. 

0 

49 

7. 

2 

284 

1.7 

111100000 

2 

55 

2, 

,2 

450 

959. 

,5 

17. 

,9 

49 

7. 

,1 

/// 

//// 

111100000 

3 

0 

2. 

.2 

461 

958. 

,4 

17. 

9 

49 

7. 

1 

285 

1.6 

111100000 

3 

5 

2, 

.1 

470 

957. 

,3 

18. 

,1 

49 

7. 

3 

/// 

//// 

111100000 

3 

10 

2, 

.0 

478 

956. 

,5 

18. 

,2 

48 

7. 

1 

296 

1.4 

111100000 

3 

15 

1. 

.9 

488 

955. 

,3 

18. 

,4 

47 

7. 

,0 

/// 

//// 

111100000 

3 

20 

1, 

.9 

498 

954. 

,2 

18. 

,4 

46 

6. 

,6 

306 

1.3 

111100000 

3 

25 

1, 

,9 

509 

953. 

,0 

18. 

,4 

46 

6. 

,6 

/// 

//// 

111100000 

3 

30 

1. 

.9 

517 

952. 

,  1 

18. 

,4 

45 

6. 

,3 

301 

1.2 

111100000 

3 

35 

1, 

.8 

525 

951. 

,1 

18. 

,5 

45 

6. 

,4 

/// 

//// 

111100000 

3 

40 

1, 

.9 

534 

950. 

,2 

18. 

,4 

44 

6. 

,0 

283 

1.2 

111100000 

3 

45 

1, 

,8 

542 

949. 

,3 

18. 

,5 

42 

5. 

,4 

/// 

//// 

111100000 

3 

50 

1. 

,8 

551 

948. 

,4 

18. 

,6 

41 

5, 

,2 

285 

1.1 

111100000 

3 

55 

1. 

.8 

561 

947. 

,2 

18. 

,7 

41 

5. 

,3 

/// 

//// 

111100000 

4 

0 

1. 

,7 

569 

946. 

,3 

18. 

,9 

42 

5. 

,8 

288 

1.1 

111100000 

4 

5 

1, 

.7 

575 

945. 

,  6 

19. 

,1 

41 

5. 

,6 

/// 

//// 

111100000 

4 

10 

1, 

.7 

584 

944. 

,7 

19. 

,1 

38 

4. 

,5 

283 

1.3 

111100000 

4 

15 

1. 

,7 

592 

943, 

,7 

19. 

,2 

37 

4, 

.2 

/// 

//// 

111100000 

4 

20 

1, 

.7 

601 

942, 

.8 

19. 

,2 

37 

4  , 

.2 

281 

1.5 

111100000 

4 

25 

1, 

,6 

609 

941, 

.9 

19. 

,1 

37 

4, 

.1 

/// 

//// 

111100000 

4 

30 

1, 

.5 

615 

941, 

.2 

19. 

,1 

36 

3, 

.8 

271 

1.8 

111100000 

4 

35 

1 

.5 

620 

940, 

.6 

19, 

,1 

36 

3, 

.8 

/// 

//// 

111100000 

4 

40 

1, 

.6 

632 

939, 

.4 

19, 

,0 

36 

3, 

.7 

265 

2.0 

111100000 

4 

45 

1, 

.5 

638 

938, 

.7 

18, 

,9 

35 

3, 

.2 

/// 

//// 

111100000 

4 

50 

1. 

.5 

647 

937, 

.8 

18, 

,9 

35 

3, 

.2 

259 

2.6 

111100000 

4 

55 

1 

,5 

653 

937, 

.1 

18, 

.9 

35 

3, 

.2 

/// 

//// 

111100000 

5 

0 

1 

.5 

659 

936, 

.4 

18, 

,8 

35 

3 

.1 

256 

3.2 

111100000 

5 

5 

1 

.6 

668 

935, 

.4 

18, 

.9 

36 

3, 

.6 

/// 

//// 

111100000 

5 

10 

1 

.3 

672 

935, 

.0 

19, 

.0 

35 

3. 

.3 

252 

3.9 

111100000 

5 

15 

1 

.5 

681 

934 

.0 

19, 

.1 

31 

1 

.7 

/// 

//// 

111100000 

5 

20 

1 

.3 

686 

933 

.4 

19, 

.2 

28 

0 

.3 

253 

3.9 

111100000 

5 

25 

1 

.3 

694 

932 

.6 

19 

.3 

24 

-1 

.7 

/// 

//// 

111100000 

5 

30 

1 

.3 

699 

932 

.1 

19 

.3 

23 

-2 

.3 

250 

4.0 

111100000 

5 

35 

1 

.2 

705 

931 

.4 

19, 

.4 

19 

-4 

.7 

/// 

//// 

111100000 

5 

40 

1 

.3 

711 

930 

.7 

19 

.5 

18 

-5 

.4 

251 

4.1 

111100000 

5 

45 

1 

.3 

718 

930 

.0 

19 

.5 

17 

-6 

.1 

/// 

//// 

111100000 

5 

50 

1 

.3 

724 

929 

.3 

19 

.5 

17 

-6 

.1 

247 

4.5 

111100000 

5 

55 

1 

.2 

730 

928 

.7 

19 

.5 

17 

-6 

.1 

/// 

//// 

111100000 

125 


Sounding   program   REV   7 . 62    using    Loran-C 

Ship  :PT_SUR 

Location  :  36.80  N  122.00  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  21 

RS-number:  108234651 

No  PTU  editing 

No  wind  editing 

Started  at:       19  MAY  95  21:55  GMT 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


Time 

AscRate  Hgt 

/MSL 

Pressure 

Temp 

RH 

Dewp 

Dir 

Speed 

n 

.  a.  u 
S 

m/s 

m 

hPa 

degC 

"5 

degC 

deg 

m/s 

0 

0 

0.0 

3 

1013.9 

11.0 

95 

10.2 

276 

5.8 

000000000 

0 

5 

3.0 

20 

1011.9 

12.0 

71 

7.0 

/// 

//// 

111100000 

0 

10 

3.5 

38 

1009.8 

11.7 

74 

7.3 

304 

6.5 

111100000 

0 

15 

3.9 

60 

1007.1 

11.5 

76 

7.5 

/// 

//// 

111100000 

0 

20 

3.5 

74 

1005.3 

11.3 

77 

7.5 

308 

5.9 

111100000 

0 

25 

3.5 

93 

1003.1 

11.1 

79 

7.6 

/// 

//// 

111100000 

0 

30 

3.5 

109 

1001.2 

11.0 

80 

7.7 

305 

5.7 

111100000 

0 

35 

3.5 

125 

999.2 

10.8 

81 

7.7 

/// 

//// 

111100000 

0 

40 

3.4 

141 

997.3 

10.6 

81 

7.5 

312 

5.4 

111100000 

0 

45 

3.2 

157 

995.3 

10.5 

82 

7.6 

/// 

//// 

111100000 

0 

50 

3.3 

172 

993.6 

10.4 

83 

7.7 

309 

5.1 

111100000 

0 

55 

3.3 

191 

991.3 

10.2 

85 

7.8 

/// 

//// 

111100000 

1 

0 

3.6 

218 

988.1 

10.0 

85 

7.6 

315 

4.6 

111100000 

1 

5 

3.9 

239 

985.5 

10.0 

84 

7.5 

/// 

//// 

111100000 

1 

10 

4.1 

263 

982.8 

10.3 

82 

7.4 

308 

4.1 

111100000 

1 

15 

4.0 

278 

981.0 

10.6 

80 

7.3 

/// 

//// 

111100000 

1 

20 

4.1 

295 

978.9 

11.1 

79 

7.6 

294 

3.7 

111100000 

1 

25 

4.0 

311 

977.0 

11.3 

81 

8.2 

/// 

//// 

111100000 

1 

30 

3.8 

332 

974.7 

11.4 

82 

8.5 

287 

3.5 

111100000 

1 

35 

3.6 

348 

972.7 

11.9 

80 

8.6 

/// 

//// 

111100000 

1 

40 

3.4 

364 

970.9 

13.8 

40 

0.5 

286 

3.2 

111100000 

1 

45 

3.4 

379 

969.2 

15.0 

41 

1.9 

/// 

//// 

111100000 

1 

50 

3.3 

395 

967.3 

15.6 

51 

5.6 

284 

2.7 

111100000 

1 

55 

3.3 

412 

965.4 

16.3 

53 

6.8 

/// 

//// 

111100000 

2 

0 

3.3 

431 

963.3 

17.3 

47 

6.0 

279 

2.5 

111100000 

2 

5 

3.2 

445 

961.7 

17.9 

49 

7.1 

/// 

//// 

111100000 

2 

10 

3.1 

458 

960.2 

17.9 

49 

7.1 

274 

2.5 

111100000 

2 

15 

3.0 

474 

958.4 

17.9 

49 

7.1 

/// 

//// 

111100000 

2 

20 

3.2 

491 

956.5 

17.8 

50 

7.3 

273 

2.3 

111100000 

2 

25 

3.2 

508 

954.5 

17.9 

49 

7.1 

/// 

//// 

111100000 

2 

30 

3.2 

526 

952.5 

18.1 

47 

6.7 

272 

2.2 

111100000 

2 

35 

2.9 

535 

951.6 

18.4 

45 

6.3 

/// 

//// 

111100000 

2 

40 

2.9 

545 

950.4 

18.7 

42 

5.6 

273 

2.1 

111100000 

2 

45 

3.2 

569 

947.8 

18.8 

41 

5.3 

/// 

//// 

111100000 
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Sounding  program  REV  7 . 62  using  Loran-C 

Ship      : PT_SUR 

Location  :  36.70  N  122.40  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  23 

RS-number:  108234843 

No  PTU  editing 

No  wind  editing 

Started  at:       20  MAY  95  15:19  GMT 


0  min  to 
10  min  to 


10  min 
4  5  min 


45  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure    Temp   RH 
WndStat 
min   s      m/s        m      hPa    degC    % 


Dewp   Dir  Speed 
degC  deg   m/s 


0  0 
0  5 
0  10 
0  15 
0  20 
0  25 
0  30 
0  35 
0  40 
0  45 
0  50 
0  55 


0.0 
6.5 


0 
5 

10 
15 
20 
25 
30 


1  35 


40 

45 
50 
55 
0 
5 
10 
15 
20 
25 
30 
35 
40 
45 


3 
31 
58 
93 
123 
156 
186 
214 
247 
269 
305 
331 
356 
391 
420 
452 
480 
513 
544 
572 
613 
629 
655 
692 
738 
760 
783 
801 
821 
830 
848 
865 
883 
902 


1014 

1011 

1008 

1004 

1000 
996 
992 
989 
985.7 
983.0 
978.7 
975.6 
972.7 
968.6 
965 
961 
958 
954 
950.7 
947.4 
942.8 
941.0 
938.0 
933.7 
928.4 
925 
923 
921 
919 
918 
916 
914 
912 
910 


11.1 

12.1 
11.6 

11.3 
10.9 
10.7 


10, 

10. 
9. 
9, 
9, 


3 

1 
7 
5 
2 
9 
7 
4 
0 
8 
6 
4 
2 
7.0 
6.8 


6, 

6, 

6. 

6, 

7. 
10.2 
13.2 
15.3 
15.9 
16.1 
16.0 
16.0 
16.0 


87 
65 
67 
71 
72 
74 
75 
77 
78 
79 
80 
82 
83 
85 
86 
88 
89 
89 
89 
90 
90 
90 
90 
90 
27 
19 
14 
10 
9 
11 
11 
11 
10 


9.0   326 


5. 

5. 

6. 

6. 

6. 

6. 

6. 

6. 

6. 

6. 

6. 

6. 

6. 

5. 

6, 

5, 

5, 
'5, 

5, 

5, 

5, 

4, 

4, 
■11. 
■14 
•16, 


7 
3 
1 
3 
1 
3 
1 
1 
0 
0 
0 

1 

8 
0 
9 
8 
6 
5 
3 
2 
9 
6 
5 
9 
1 

■17.7 
■17.4 


•14, 

■14, 

■14 

■15 

•18 


/// 
326 
/// 
329 
/// 
325 
/// 
320 
/// 
311 
/// 
316 
/// 
323 
/// 
330 
/// 
324 
/// 
344 
/// 
309 
/// 
265 
/// 
241 
/// 
237 
/// 
230 
/// 
226 
/// 


4.2  000000000 
////  111100000 

3.7  111100000 
////  111100000 

3.5  111100000 
////  111100000 

3.5  111100000 
////  111100000 

3.5  111100000 
////  111100000 

3.5  111100000 
////  111100000 

3.3  111100000 
////  111100000 

3.2  111100000 
////  111100000 

3.0  111100000 
////  111100000 

2.6  111100000 
////  111100000 

2.4  111100000 
////  111100000 

1.6  111100000 
////  111100000 

1.9  111100000 
////  111100000 

1.8  111100000 
////  111100000 

1.9  111100000 
////  111100000 

2.1  111100000 
////  111100000 

2.3  111100000 
////  111100000 
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Sounding  program  REV  7 . 62  using  Loran-C 

Ship      :PT_SUR 

Location  :  36.70  N  122.40  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  24 

RS-number:  108234747 

No  PTU  editing 

No  wind  editing 

Sounding  program  REV  7 . 62  using  Loran-C 

Ship      : PT_SUR 

Location  :  36.70  N  122.40  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  24 

RS-number:  108234747 

No  PTU  editing 

No  wind  editing 

Started  manually  by  operator 

Started  at:       20  MAY  95  17:26  GMT 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


0  min  to 
10  min  to 


10  min 
45  min 


45  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure 
WndStat 
min   s      m/s        m      hPa 


Temp   RH    Dewp   Dir  Speed 


degC 


degC   deg    m/s 


0 

0 

0. 

0 

3 

1015. 

2 

11. 

2 

84 

8. 

6 

355 

2.3 

000000000 

0 

5 

175. 

5 

1057 

895. 

4 

15. 

6 

11 

-14. 

7 

/// 

//// 

000100000 

0 

10 

104. 

2 

1045 

896. 

6 

15. 

8 

11 

-14. 

6 

/// 

//// 

000100000 

0 

15 

64. 

7 

1041 

897. 

0 

15. 

9 

11 

-14  . 

5 

/// 

//// 

011100000 

0 

20 

50. 

9 

1021 

899. 

2 

16. 

0 

10 

-15. 

6 

204 

5.6 

011100000 

0 

25 

39. 

0 

1017 

899. 

5 

16. 

0 

10 

-15. 

6 

/// 

//// 

111100000 

0 

30 

33. 

4 

1006 

900. 

7 

16. 

0 

10 

-15. 

,6 

211 

3.2 

111100000 

0 

35 

-2. 

1 

997 

901. 

7 

16. 

0 

11 

-14  . 

,4 

/// 

//// 

111100000 

0 

40 

-1. 

,8 

990 

902. 

,5 

15. 

9 

11 

-14. 

,5 

225 

4.8 

111100000 

0 

45 

-2, 

,1 

979 

903. 

,6 

15. 

9 

11 

-14. 

,5 

/// 

//// 

111100000 

0 

50 

-1. 

,7 

969 

904, 

,7 

16. 

,0 

12 

-13. 

,3 

232 

4.  1 

111100000 

0 

55 

-2, 

,0 

960 

905. 

,7 

16. 

,1 

11 

-14, 

.3 

/// 

//// 

111100000 

1 

0 

-1. 

,9 

950 

906. 

,7 

16. 

,1 

11 

-14, 

.3 

222 

2.2 

111100000 

1 

5 

-1, 

.8 

941 

907, 

.7 

16. 

,2 

11 

-14, 

.3 

/// 

//// 

111100000 

1 

10 

-2, 

.0 

930 

908, 

.9 

16. 

.3 

11 

-14, 

.2 

231 

2.8 

111100000 

1 

15 

-1, 

.9 

921 

909, 

.8 

16, 

.3 

11 

-14 

.2 

/// 

//// 

111100000 

1 

20 

-1 

.9 

913 

910, 

.7 

16, 

.3 

12 

-13, 

.1 

231 

2.0 

111100000 

1 

25 

-1 

.9 

902 

911 

.9 

16, 

.4 

13 

-12 

.0 

/// 

//// 

111100000 

1 

30 

-1 

.9 

892 

912 

.9 

16 

.4 

12 

-13 

.0 

232 

2.5 

111100000 

1 

35 

-2 

.0 

882 

914 

.0 

16 

.5 

12 

-12 

.9 

/// 

//// 

111100000 

1 

40 

-1 

.8 

876 

914 

.7 

16 

.5 

12 

-12 

.9 

234 

3.0 

111100000 

1 

45 

-1 

.8 

865 

915 

.8 

16 

.6 

12 

-12 

.9 

/// 

//// 

111100000 

1 

50 

-1 

.9 

855 

916 

.9 

16 

.5 

12 

-12 

.9 

235 

3.2 

111100000 

1 

55 

-2 

.0 

844 

918 

.2 

16 

.5 

12 

-12 

.9 

/// 

//// 

111100000 
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Sounding  program  REV  7.62  using  Loran-C 

Ship      : PT_SUR 

Location  :  36.70  N  122.40  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  25 

RS-number:  108234749 

No  PTU  editing 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


No  wind 

editing 

Grour 

id  check 

Ref 

RS    C 

orr 

Pressure 

1016 

.2  1016.2 

0.0 

Temperature 

22 

.7 

22.7 

0.0 

Humidity 

42 

42 

0 

Started 

manually 

by 

operator 

Started 

at: 

20 

MAY 

95  17:52 

GMT 

Time 

AscRate 

Hgt/MSL 

Pressure 

Temp 

RH 

Dewp 

Dir 

Speed 

min 

-UL 

S 

m/s 

m 

hPa 

degC 

Q. 
0 

degC 

deg 

m/s 

0 

0 

0.0 

3 

1016.2 

11.2 

88 

9.3 

25 

1.9 

000000000 

0 

5 

0.0 

2 

1016.3 

12.9 

73 

8.2 

/// 

//// 

111100000 

0 

10 

-0.2 

1 

1016.5 

12.9 

72 

8.0 

219 

2.2 

111100000 

0 

15 

-0.1 

1 

1016.5 

12.9 

72 

8.0 

/// 

//// 

111100000 

0 

20 

-0.1 

1 

1016.5 

12.8 

73 

8.1 

228 

2.8 

111100000 

0 

25 

-0.2 

1 

1016.5 

•  12.8 

73 

8.1 

/// 

//// 

111100000 

0 

30 

-0.1 

1 

1016.5 

12.7 

73 

8.0 

234 

2.7 

111100000 

0 

35 

-0.1 

1 

1016.5 

12.7 

74 

8.2 

/// 

//// 

111100000 

0 

40 

0.0 

1 

1016.5 

12.6 

74 

8.1 

242 

2.4 

111100000 

0 

45 

0.0 

1 

1016.5 

12.6 

74 

8.1 

/// 

//// 

111100000 

0 

50 

0.0 

1 

1016.5 

12.6 

74 

8.1 

287 

0.8 

111100000 

0 

55 

0.1 

1 

1016.5 

12.6- 

74 

8.1 

/// 

//// 

111100000 

1 

0 

0.0 

1 

1016.5 

12.5 

74 

8.0 

324 

1.0 

111100000 

1 

5 

-0.1 

-2 

1016.8 

12.5 

75 

8.2 

/// 

//// 

111100000 

1 

10 

-0.1 

-1 

1016.7 

12.5 

75 

8.2 

325 

1.6 

111100000 

1 

15 

-0.1 

-2 

1016.8 

12.4 

74 

7.9 

/// 

//// 

111100000 

1 

20 

-0.1 

-1 

1016.7 

12.4 

77 

8.5 

323 

1.9. 

111100000 

1 

25 

-0.1 

-1 

1016.7 

12.3 

77 

8.4 

/// 

//// 

111100000 

1 

30 

-0.1 

-1 

1016.7 

12.3 

77 

8.4 

311 

2.2 

111100000 

1 

35 

0.0 

-3 

1016.9 

12.5 

78 

8.8 

/// 

//// 

111100000 

1 

40 

0.1 

3 

1016.2 

12.2 

77 

8.3 

316 

2.5 

111100000 

1 

45 

0.4 

8 

1015.6 

12.4 

77 

8.5 

/// 

//// 

111100000 

1 

50 

0.3 

9 

1015.5 

12.3 

76 

8.2 

318 

2.7 

111100000 

1 

55 

0.3 

8 

1015.6 

12.4 

76 

8.3 

/// 

//// 

111100000 

2 

0 

0.3 

9 

1015.5 

12.6 

76 

8.5 

318 

2.7 

111100000 

2 

5 

0.4 

9 

1015.5 

12.4 

75 

8.1 

/// 

//// 

111100000 

2 

10 

0.1 

7 

1015.7 

12.5 

75 

8.2 

314 

2.9 

111100000 

2 

15 

0.1 

10 

1015.3 

12.5 

76 

8.4 

/// 

//// 

111100000 

2 

20 

-0.1 

7 

1015.7 

12.4 

75 

8.1 

313 

3.2 

111100000 

2 

25 

-0.1 

7 

1015.7 

12.6 

76 

8.5 

/// 

//// 

111100000 

2 

30 

-0.1 

7 

1015.7 

12.6 

76 

8.5 

311 

3.2 

111100000 

2 

35 

-0.1 

5 

1016.0 

12.7 

75 

8.4 

/// 

//// 

111100000 

2 

40 

-0.2 

1 

1016.5 

12.6 

75 

8.3 

312 

3.1 

111100000 

2 

45 

-0.4 

0 

1016.6 

12.6 

75 

8.3 

/// 

//// 

111100000 

2 

50 

-0. 

3 

-1 

1016.7 

12.9 

75 

8.6 

317 

3.1 

111100000 
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Sounding  program  REV  7 . 62  using  Loran-C 

Ship      :PT_SUR 

Location  :  36.70  N  122.40  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  25 

RS-number:  108234749 

No  PTU  editing 

No  wind  editing 

Ground  check   :     Ref      RS 

Pressure 

Temperature 

Humidity     :      71      71 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


1016.9  1016.9 
13.3    13.3 


Corr 

0.0 

-0.0 

0 


Started  manually  by  operator 
Started  at:       20  MAY  95  19:54  GMT 


Time  AscRate  Hgt/MSL  Pressure 
WndStat 
min   s      m/s        m      hPa 


Temp   RH    Dewp   Dir  Speed 


degC 


degC   deg   m/s 


0 

0 

0. 

0 

3 

1016. 

9 

13.3 

71 

8 

2 

124 

0.5 

000000000 

0 

5 

0. 

0 

3 

1016. 

9 

13.2 

72 

8 

3 

/// 

//// 

111100000 

0 

10 

0. 

4 

7 

1016. 

5 

13.1 

72 

8 

2 

218 

1.2 

111100000 

0 

15 

0. 

0 

3 

1016. 

9 

13.2 

72 

8 

3 

/// 

//// 

111100000 

0 

20 

0. 

0 

3 

1016. 

9 

13.0 

71 

7 

9 

217 

1.2 

111100000 

0 

25 

0. 

0 

4 

1016. 

8 

13.0 

72 

8 

1 

/// 

//// 

111100000 

0 

30 

0. 

0 

3 

1016. 

9 

12.7 

70 

7 

4 

218 

1.1 

111100000 

0 

35 

0. 

0 

2 

1017. 

1 

13.0 

71 

7 

9 

/// 

//// 

111100000 

0 

40 

-0. 

1 

3 

1016. 

9 

13.1 

71 

8 

0 

216 

0.8 

111100000 

0 

45 

-0. 

1 

3 

1016. 

9 

13.0 

71 

7 

9 

/// 

//// 

111100000 

0 

50 

0. 

0 

3 

1016. 

9 

13.4 

72 

8 

5 

208 

0.5 

111100000 

0 

55 

0. 

1 

4 

1016. 

8 

13.2 

71 

8 

1 

/// 

//// 

111100000 

1 

0 

0 

0 

3 

1016. 

9 

12.9 

71 

7 

8 

208 

0.6 

111100000 

1 

5 

0. 

0 

3 

1016. 

9 

12.8 

72 

7 

9 

/// 

//// 

111100000 

1 

10 

0 

0 

3 

1016 

9 

12.9 

75 

8 

6 

187 

0.5 

111100000 

1 

15 

0 

1 

5 

1016. 

7 

13.0 

73 

8 

3 

/// 

//// 

111100000 

1 

20 

0 

1 

7 

1016 

5 

12.8 

71 

7 

7 

119 

0.6 

111100000 

1 

25 

0 

0 

5 

1016 

7 

13.1 

72 

8 

2 

/// 

//// 

111100000 

1 

30 

0 

1 

7 

1016 

5 

12.9 

70 

7 

6 

139 

0.6 

111100000 

1 

35 

0 

2 

8 

1016 

3 

12.8 

69 

7 

3 

/// 

//// 

111100000 

1 

40 

0 

1 

7 

1016 

5 

12.8 

69 

7 

3 

106 

0.7 

111100000 

1 

45 

0 

1 

9 

1016 

2 

12.7 

70 

7 

4 

/// 

//// 

111100000 

1 

50 

0 

1 

9 

1016 

2 

12.6 

70 

7 

.3 

133 

0.5 

111100000 

1 

55 

0 

1 

7 

1016 

5 

12.7 

70 

7 

.4 

/// 

//// 

111100000 

2 

0 

0 

1 

9 

1016 

2 

12.6 

71 

7 

.5 

99 

0.8 

111100000 

2 

5 

0 

.1 

9 

1016 

2 

12.7 

70 

7 

.4 

/// 

//// 

111100000 

2 

10 

0 

.3 

15 

1015 

.5 

12.6 

70 

7 

.3 

108 

0.8 

111100000 

2 

15 

0 

.3 

17 

1015 

.2 

12.6 

72 

7 

.7 

/// 

//// 

111100000 

2 

20 

0 

.2 

15 

1015 

.5 

12.5 

72 

7 

.6 

95 

1.0 

111100000 

2 

25 

0 

.4 

18 

1015 

.  1 

12.5 

71 

7 

.4 

/// 

//// 

111100000 

2 

30 

0 

.3 

17 

1015 

.2 

12.5 

72 

7 

.6 

98 

0.7 

111100000 

2 

35 

0 

.3 

19 

1015 

.0 

12.6 

72 

7 

.7 

/// 

//// 

111100000 

2 

40 

0 

.2 

21 

1014 

.7 

12.3 

72 

7 

.5 

136 

0.5 

111100000 

2 

45 

0 

.3 

25 

1014 

.2 

12.3 

72 

7 

.5 

/// 

//// 

111100000 

2 

50 

0 

.6 

32 

1013 

.5 

12.3 

72 

7 

.5 

112 

0.5 

111100000 
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Start  Up  Date  22  MAY  95  15:42  GMT 

System  test  passed   -No  errors  found 

Sounding  program  REV  7.62  using  Loran-C 

Ship      :PT_SUR 

Location  :  36.80  N  121.80  W 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  31 

RS-number:  108234842 

No  PTU  editing 

No  wind  editing 

Started  at :   .     22  MAY  95  15:52  GMT 


3 

m 

from 

0  min  to 

10  min 

from 

10  min  to 

45  min 

45  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure 
WndStat 
min   s      m/s        m      hPa 


Temp   RH    Dewp   Dir  Speed 


degC 


degC   deg    m/s 


0 

0 

0.0 

3 

1015. 

1 

11. 

2 

94 

10.3 

279 

0 

5 

3.7 

23 

1012. 

7 

12. 

1 

75 

7.9 

/// 

0 

10 

3.7 

40 

1010. 

8 

11. 

9 

77 

8.0 

270 

0 

15 

3.4 

56 

1008. 

8 

11. 

7 

79 

8.2 

/// 

0 

20 

3.3 

70 

1007. 

1 

11. 

6 

80 

8.3 

273 

0 

25 

3.2 

81 

1005. 

7 

11. 

,5 

81 

8.4 

/// 

0 

30 

3.1 

97 

1003. 

9 

11. 

,3 

82 

8.4 

277 

0 

35 

3.1 

115 

1001. 

,7 

11. 

2 

83 

8.5 

/// 

0 

40 

3.1 

133 

999. 

,5 

11. 

,0 

84 

8.4 

285 

0 

45 

3.1 

149 

997. 

,5 

10. 

,9 

85 

8.5 

/// 

0 

50 

3.3 

168 

995. 

,3 

10. 

,7 

86 

8.5 

276 

0 

55 

3.5 

188 

992, 

,9 

10, 

,5 

85 

8.1 

/// 

1 

0 

3.7 

208 

990. 

,5 

10, 

,3 

82 

7.4 

275 

1 

5 

4.0 

229 

987. 

,9 

10, 

.2 

83 

7.5 

/// 

1 

10 

3.7 

245 

986, 

.1 

10, 

,1 

82 

7.2 

275 

1 

15 

3.6 

259 

984, 

,5 

9. 

,9 

83 

7.2 

/// 

1 

20 

3.5 

273 

982, 

,8 

9, 

,7 

83 

7.0 

275 

1 

25 

3.5 

293 

980, 

,4 

9. 

,6 

84 

7.1 

/// 

1 

30 

3.4 

311 

978, 

.2 

9, 

,5 

81 

6.5 

289 

1 

35 

3.2 

326 

976, 

.6 

9, 

,5 

81 

6.5 

/// 

1 

40 

3.2 

340 

974, 

.9 

9, 

.4 

81 

6.4 

292 

1 

45 

3.1 

352 

973 

.5 

9, 

.3 

81 

6.3 

/// 

1 

50 

3.2 

370 

971 

.3 

9 

.1 

81 

6.1 

290 

1 

55 

3.0 

381 

970 

.0 

9 

.0 

82 

6.1 

/// 

2 

0 

2.8 

394 

968 

.5 

8 

.8 

82 

5.9 

281 

2 

5 

2.7 

409 

966 

.7 

8 

.7 

83 

6.0 

/// 

2 

10 

2.7 

422 

965 

.2 

8 

.5 

83 

5.8 

281 

2 

15 

2.8 

434 

963 

.8 

8 

.5 

84 

6.0 

/// 

2 

20 

2.6 

447 

962 

.4 

8 

.3 

85 

6.0 

290 

4.1  000000000 

////  111100000 

4.4  111100000 
////  111100000 

4.3  111100000 
////  111100000 

4.0  111100000 

////  111100000 

3.5  111100000 
////  111100000 

3.9  111100000 

////  111100000 

4.0  111100000 

////  111100000 

3.8  111100000 

////  111100000 

3.6  111100000 
////  111100000 

2.7  111100000 
////  111100000 

2.4  111100000 
////  111100000 

2.5  111100000 
////  111100000 

2.6  111100000 
////  111100000 

2.6  111100000 

////  111100000 

2.3  111100000 
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Sounding  program  REV  7.62  using  Loran-C 


Ship      :PT_SUR 

Location  :  36.80  N  121.80  W 


3  m 


Loran-C  chain  1 


9940 


Phase  fitting  length  is  60  s  from  0  min  to  10  min 
Phase  fitting  length  is  120  s  from  10  min  to  45  min 
Phase  fitting  length  is  240  s  from   45  min  to  120  min 


Sounding  :  32 
RS-number:  108234754 


No  PTU  editing 
No  wind  editing 

Started  at:       22  MAY  95  17:57  GMT 

Time  AscRate  Hgt/MSL  Pressure    Temp 
WndStat 
min   s     m/s       m      hPa   degC 


RH    Dewp   Dir  Speed 
%    degC   deg   m/s 


0 

0 

0. 

0 

3 

1016.7 

11. 

4 

89 

9.7 

306 

0 

5 

5. 

8 

29 

1013.5 

12. 

6 

70 

7.3 

/// 

0 

10 

5. 

1 

54 

1010.5 

12. 

1 

73 

7.5 

313 

0 

15 

5. 

2 

80 

1007.3 

11. 

9 

74 

7.5 

/// 

0 

20 

4. 

6 

95 

1005.6 

11. 

7 

75 

7.5 

326 

0 

25 

4. 

8 

123 

1002.2 

11. 

4 

77 

7.6 

/// 

0 

30 

4. 

9 

150 

999.0 

11. 

1 

79 

7.6 

317 

0 

35 

4. 

7 

173 

996.2 

10. 

9 

79 

7.5 

/// 

0 

40 

4  . 

6 

192 

993.9 

10. 

7 

78 

7.1 

332 

0 

45 

4. 

4 

216 

991.1 

10. 

.6 

78 

7.0 

/// 

0 

50 

5. 

1 

247 

987.3 

10. 

,4 

77 

6.6 

318 

0 

55 

4. 

,7 

264 

985.3 

10. 

,3 

77 

6.5 

/// 

1 

0 

4  , 

,4 

282 

983.3 

10. 

,2 

76 

6.2 

309 

1 

5 

4  . 

,1 

301 

981.0 

10. 

,2 

76 

6.2 

/// 

1 

10 

4. 

,5 

326 

978.0 

10. 

,0 

77 

6.2 

303 

1 

15 

4. 

,5 

351 

975.0 

9. 

.7 

78 

6.1 

/// 

1 

20 

4. 

,3 

375 

972.3 

9. 

.5 

79 

6.1 

298 

1 

25 

4. 

.5 

397 

969.7 

9, 

.4 

80 

6.2 

/// 

1 

30 

4. 

.4 

415 

967.5 

9, 

.2 

81 

6.2 

291 

1 

35 

4, 

.6 

436 

965.1 

9, 

.0 

82 

6.1 

/// 

1 

40 

4, 

.7 

467 

961.4 

8 

.7 

84 

6.2 

277 

1 

45 

4 

.9 

496 

958.1 

8 

.5 

86 

6.3 

/// 

1 

50 

4 

.6 

512 

956.3 

8 

.3 

87 

6.3 

285 

1 

55 

4 

.4 

532 

953.9 

8 

.1 

88 

6.3 

/// 

2 

0 

4 

.5 

550 

951.8 

8 

.0 

89 

6.3 

281 

2 

5 

4 

.5 

572 

949.3 

7 

.8 

90 

6.3 

/// 

2 

10 

4 

.2 

592 

947.0 

7 

.7 

90 

6.2 

274 

2 

15 

3 

.7 

607 

945.2 

7 

.6 

91 

6.3 

/// 

2 

20 

3 

.8 

626 

943.1 

7 

.5 

91 

6.2 

261 

2 

25 

3 

.9 

650 

940.3 

7 

.4 

91 

6.1 

/// 

2 

30 

3 

.9 

668 

938.2 

7 

.3 

91 

6.0 

257 

2 

35 

3 

.9 

687 

936.1 

7 

.2 

91 

5.9 

/// 

2 

40 

3 

.8 

706 

933.9 

7 

.1 

91 

5.8 

249 

2 

45 

3 

.9 

726 

931.6 

7 

.0 

92 

5.8 

/// 

1.6  000000000 
////  110100000 

4.5  110100000 
////  110100000 

4.2  110100000 
////  110100000 

4.2  110100000 
////  110100000 

4.0  110100000 
////  110100000 

3.8  110100000 
////  110100000 

3.6  110100000 
////  110100000 

3.5  110100000 
////  110100000 

3.1  110100000 
////  110100000 

3.1  110100000 
////  110100000 

3.3  110100000 
////  111100000 

2.1  111100000 
////  111100000 

1.8  111100000 
////  111100000 

1.9  111100000 
////  111100000 

2.2  111100000 
////  111100000 

2.3  111100000 
////  111100000 

2.4  111100000 
////  111100000 
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May  23 

1600 

0 

5 

4. 

,5 

23 

1011. 

5 

12. 

9 

69 

7. 

4 

/// 

//// 

— CD-FGH 

0 

10 

4. 

7 

50 

1008. 

3 

12. 

2 

72 

7. 

4 

263 

2.6 

— CD-FGH 

0 

15 

4. 

6 

71 

1005. 

7 

12. 

0 

74 

7. 

6 

/// 

//// 

— CD-FGH 

0 

20 

4. 

6 

95 

1002. 

9 

11. 

7 

75 

7. 

5 

228 

2.3 

— CD-FGH 

0 

25 

4. 

5 

115 

1000. 

4 

11. 

6 

76 

7. 

6 

/// 

//// 

— CD-FGH 

0. 

30 

4. 

4 

136 

998. 

0 

11. 

3 

77 

7. 

5 

226 

2.3 

— CD-FGH 

0 

35 

4. 

3 

153 

995. 

9 

11. 

2 

78 

7. 

6 

/// 

//// 

— CD-FGH 

0 

40 

4. 

0 

170 

993. 

9 

11. 

0 

78 

7. 

4 

234 

2.2 

— CD-FGH 

0 

45 

4. 

0 

191 

991. 

3 

10. 

8 

79 

7. 

4 

/// 

//// 

— CD-FGH 

0 

50 

3. 

8 

209 

989. 

3 

10. 

6 

80 

7. 

3 

231 

2.2 

— CD-FGH 

0 

55 

3. 

9 

230 

986. 

7 

10. 

4 

81 

7. 

3 

/// 

//// 

-BCD-FGH 

1 

0 

3. 

7 

247 

984. 

7 

10. 

2 

81 

7. 

1 

222 

2.2 

-BCD-FGH 

1 

5 

4. 

3 

276 

981. 

2 

10. 

0 

83 

7. 

3 

/// 

//// 

-BCD-FGH 

1 

10 

4. 

3 

298 

978. 

7 

9. 

7 

84 

7. 

2 

230 

2.0 

-BCD-FGH 

1 

15 

4. 

1 

315 

976. 

7 

9. 

6 

86 

7. 

4 

/// 

//// 

-BCD-FGH 

1 

20 

4. 

0 

330 

974. 

9 

9. 

5 

86 

7. 

3 

214 

2.1 

-BCD-FGH 

1 

25 

3. 

8 

•  350 

972. 

5 

9. 

4 

87 

7. 

4 

/// 

//// 

-BCD-FGH 

1 

30 

3. 

8 

362 

971. 

1 

9. 

2 

87 

7. 

2 

226 

1.8 

-BCD-FGH 

1 

35 

3. 

9 

394 

967. 

4 

9. 

1 

87 

7. 

1 

/// 

//// 

-BCD-FGH 

1 

40 

3. 

8 

413 

965. 

2 

9. 

0 

87 

7. 

0 

228 

1.6 

-BCD-FGH 

1 

45 

4. 

0 

436 

962. 

5 

8. 

9 

88 

7. 

1 

/// 

//// 

-BCD-FGH 

1 

50 

4. 

,2 

457 

960. 

0 

8. 

8 

88 

7. 

0 

233 

1.5 

-BCD-FGH 

1 

55 

4. 

A 

471 

958. 

4 

8. 

,7 

88 

6. 

9 

/// 

//// 

-BCD-FGH 

2 

0 

4. 

,5 

497 

955. 

3 

8. 

,6 

88 

6. 

8 

240 

1.3 

-BCD-FGH 

2 

5 

4. 

,0 

517 

953. 

1 

8. 

,5 

88 

6. 

7 

/// 

//// 

-BCD-FGH 

2 

10 

4. 

,5 

548 

949. 

5 

8. 

,4 

88 

6. 

6 

253 

1.1 

-BCD-FGH 

2 

15 

4, 

.4 

567 

947. 

3 

8. 

,3 

88 

6. 

,5 

/// 

//// 

-BCD-FGH 

2 

20 

4, 

,6 

594 

944. 

2 

8. 

,1 

88 

6. 

,3 

259 

1.0 

-BCD-FGH 

2 

25 

4. 

,3 

612 

942. 

1 

8. 

,1 

88 

6. 

,3 

/// 

//// 

-BCD-FGH 

2 

30 

4. 

,7 

638 

939. 

,1 

7. 

,9 

87 

5. 

,9 

273 

0.8 

-BCD-FGH 

2 

35 

4, 

,3 

647 

938. 

,1 

7. 

,9 

87 

5. 

,9 

/// 

//// 

-BCD-FGH 

2 

40 

4. 

,3 

678 

934. 

,6 

8. 

,3 

89 

6. 

,6 

303 

0.8 

-BCD-FGH 

2 

45 

4. 

.7 

708 

931. 

,3 

9. 

,8 

90 

8, 

,3 

/// 

//// 

-BCD-FGH 

2 

50 

4, 

,2 

721 

929. 

,8 

10. 

,5 

90 

9. 

.0 

331 

1.1 

-BCD-FGH 

2 

55 

4. 

,2 

737 

928. 

,0 

10, 

.8 

90 

9, 

.3 

/// 

//// 

-BCD-FGH 

3 

0 

3, 

.9 

755 

925. 

,9 

10. 

.8 

87 

8, 

.8 

345 

1.5 

-BCD-FGH 

3 

5 

4, 

.2 

772 

924. 

,0 

10.. 

,8 

85 

8, 

.4 

/// 

//// 

-BCD-FGH 

3 

10 

3, 

.6 

786 

922, 

,6 

11, 

.1 

85 

8, 

.7 

350 

1.9 

-BCD-FGH 

3 

15 

3 

.6 

811 

919, 

,7 

11, 

.6 

84 

9 

.0 

/// 

//// 

-BCD-FGH 

3 

20 

3 

.8 

835 

917, 

.2 

12, 

.1 

82 

9 

.2 

356 

2.6 

-BCD-FGH 

3 

25 

3 

.8 

847 

915, 

.8 

12 

.2 

80 

8 

.9 

/// 

//// 

-BCD-FGH 

3 

30 

3 

.6 

863 

914, 

.0 

12 

.3 

79 

8 

.8 

6 

2.9 

-BCD-FGH 

3 

35 

3 

.7 

883 

911, 

.8 

12 

.3 

80 

9 

.0 

/// 

//// 

-BCD-FGH 

3 

40 

3 

.7 

896 

910 

.5 

12, 

.5 

81 

9 

.4 

6 

3.5 

-BCD-FGH 

3 

45 

3 

.4 

914 

908 

.5 

12 

.7 

81 

9 

.6 

/// 

//// 

-BCD-FGH 

3 

50 

3 

.3 

933 

906 

.5 

12 

.8 

80 

9 

.5 

6 

4.0 

-BCD-FGH 

3 

55 

3 

.5 

953 

904 

.3 

13 

.0 

80 

9 

.7 

/// 

//// 

-BCD-FGH 

4 

0 

3 

.9 

980 

901 

.4 

13 

.0 

81 

9 

.9 

5 

4.4 

-BCD-FGH 

4 

5 

3 

.4 

991 

900 

.2 

12 

.9 

80 

9 

.6 

/// 

//// 

-BCD-FGH 

4 

10 

3 

.7 

1008 

898 

.3 

12 

.8 

79 

9 

.3 

360 

5.5 

-BCD-FGH 

4 

15 

4 

.0 

1040 

894 

.9 

13 

.1 

79 

9 

.6 

/// 

//// 

-BCD-FGH 

4 

20 

4 

.7 

1074 

891 

.3 

13 

.0 

80 

9 

.7 

359 

6.1 

-BCD-FGH 

4 

25 

4 

.5 

1091 

889 

.5 

13 

.0 

79 

9 

.5 

/// 

//// 

-BCD-FGH 

4 

30 

4 

.0 

1100 

888 

.5 

13 

.1 

74 

8 

.6 

358 

6.8 

-BCD-FGH 

4 

35 

4 

.5 

1121 

886 

.3 

13 

.1 

72 

8 

.2 

/// 

//// 

-BCD-FGH 

4 

40 

4 

.4 

1141 

884 

.2 

13 

.3 

72 

8 

.4 

357 

7.1 

-BCD-FGH 

133 


Sounding  program  REV  7 . 62  using  Loran-C 

Ship      :PT_SUR 

Location  :  36.80  N  121.80  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is 


60 
Phase  fitting  length  is  120 
Phase  fitting  length  is  240 
Sounding  :  35 
RS-number:  108234749 
No  PTU  editing 
No  wind  editing 

Start  Up  Date  23  MAY  95  16:57  GMT 
System  test  passed   -No  errors  found 
Sounding  program  REV  7 . 62  using  Loran-C 
Ship      : PT_SUR 

Location  :  36.80  N  121.80  W       3m 
Loran-C  chain  1:   9940 
Phase  fitting  length  is   60  s  from 
Phase  fitting  length  is  120  s  from 
Phase  fitting  length  is  240  s  from 
Sounding  :  35 
RS-number:  108234749 
No  PTU  editing 
Started  manually  by  operator 
Started  at:       23  MAY  95  17:05  GMT 


s  from  0  min  to  10  min 
s  from  10  min  to  45  min 
s  from   45  min  to  120  min 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure    Temp   RH 
WndStat 
min   s      m/s        m      hPa    degC    % 


Dewp   Dir  Speed 
degC   deg    m/s 


0 

0 

0. 

0 

3 

1015. 

0 

12. 

,7 

76 

8. 

6 

184 

0. 

7 

000000000 

0 

5 

-0. 

7 

-1 

1015. 

5 

12. 

,8 

77 

8. 

9 

134 

2. 

6 

111100000 

0 

10 

-0. 

6 

-3 

1015. 

7 

12. 

,7 

78 

9. 

0 

109 

3. 

6 

111100000 

0 

15 

-0. 

3 

-2 

1015. 

6 

12. 

,8 

79 

9. 

3 

97 

3. 

,9 

111100000 

0 

20 

-0. 

2 

-1 

1015. 

5 

12. 

,7 

78 

9. 

0 

86 

3. 

,9 

111100000 

0 

25 

-0. 

1 

2 

1015. 

1 

12. 

,8 

78 

9. 

1 

82 

3. 

,6 

111100000 

0 

30 

-0. 

1 

1 

1015. 

2 

12. 

,7 

78 

9. 

0 

82 

3. 

2 

111100000 

0 

35 

0. 

1 

2 

1015. 

1 

12. 

,7 

78 

9. 

0 

85 

2. 

,9 

111100000 

0 

40 

0. 

.1 

1 

1015. 

,2 

12. 

.6 

78 

8. 

9 

90 

2. 

,4 

111100000 

0 

45 

0, 

,1 

2 

1015. 

,1 

12. 

,7 

79 

9. 

,2 

95 

2. 

,1 

111100000 

0 

50 

0. 

,1 

1 

1015. 

,2 

12. 

.7 

78 

9. 

,0 

100 

1. 

.7 

111100000 

0 

55 

0. 

,0 

2 

1015. 

,1 

12, 

.7 

79 

9. 

,2 

100 

1. 

.8 

111100000 

1 

0 

0, 

.1 

3 

1015. 

,0 

13, 

.1 

79 

9. 

,6 

101 

1, 

.6 

111100000 

1 

5 

0. 

.1 

6 

1014, 

,6 

13, 

.2 

78 

9. 

,5 

103 

1, 

.5 

111100000 

1 

10 

0, 

.2 

7 

1014. 

.5 

13, 

.3 

77 

9. 

,4 

103 

1, 

,4 

111100000 

1 

15 

0. 

,1 

8 

1014, 

,3 

13, 

.5 

77 

9. 

,6 

103 

1, 

.3 

111100000 

1 

20 

0, 

.1 

5 

1014. 

,7 

14, 

.0 

76 

9. 

,9 

101 

1, 

.2 

111100000 

1 

25 

0, 

.2 

7 

1014, 

.5 

13 

.6 

76 

9, 

.5 

99 

1, 

.1 

111100000 

1 

30 

0 

.2 

9 

1014, 

.2 

13, 

.0 

75 

8, 

.7 

96 

1, 

.0 

111100000 

1 

35 

0 

.3 

12 

1013, 

.8 

12 

.6 

77 

8, 

.7 

95 

1, 

.0 

111100000 

1 

40 

0 

.1 

11 

1014 

.0 

12 

.5 

77 

8 

.6 

94 

1 

.0 

111100000 

1 

45 

0 

.3 

14 

1013 

.6 

12 

.4 

79 

8 

.9 

93 

1 

.0 

111100000 

1 

50 

0 

.3 

15 

1013 

.5 

12 

.4 

81 

9 

.3 

93 

0 

.9 

111100000 

1 

55 

0 

.2 

13 

1013 

.7 

12 

.3 

81 

9 

.2 

93 

0 

.9 

111100000 

2 

0 

0 

.1 

13 

1013 

.7 

12 

.2 

80 

8 

.9 

93 

0 

.9 

111100000 

2 

5 

-0 

.1 

12 

1013 

.8 

12 

.3 

81 

9 

.2 

93 

0 

.9 

111100000 

2 

10 

0 

.0 

11 

1014 

.0 

12 

.3 

81 

9 

.2 

94 

0 

.9 

111100000 
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Sounding  program  REV  7 . 62  using  Loran-C 

Ship      :PT_SUR 

Location  :  36.80  N  121.86  W       3m 

Loran-C  chain  1:   9940 

Phase  fitting  length  is   60  s  from 

Phase  fitting  length  is  120  s  from 

Phase  fitting  length  is  240  s  from 

Sounding  :  36 

RS-number:  10a234749 

No  PTU  editing 

Started  manually  by  operator 

Started  at:       23  MAY  95  19:16  GMT 


0  min  to  10  min 
10  min  to  45  min 
45  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure    Temp   RH 
WndStat 
min   s     m/s       m      hPa    degC    % 


Dewp   Dir  Speed 
degC   deg   m/s 


0 

0 

0. 

0 

3 

1015. 

8 

12. 

6 

86 

10. 

3 

269 

3. 

0 

000000000 

0 

5 

0. 

0 

2 

1015. 

8 

12. 

6 

85 

10. 

2 

243 

2. 

0 

111100000 

0 

10 

0. 

0 

3 

1015. 

7 

12. 

6 

84 

10. 

0 

219 

1. 

4 

111100000 

0 

15 

0. 

0 

4 

1015. 

6 

12. 

6 

84 

10. 

0 

215 

1. 

1 

111100000 

0 

20 

0. 

0 

3 

1015. 

7 

12. 

6 

84 

10. 

0 

199 

0. 

8 

111100000 

0 

25 

-0. 

1 

1 

1016. 

0 

12. 

5 

84 

9. 

9 

165 

0. 

6 

111100000 

0 

30 

-0. 

1 

1 

1016. 

0 

12. 

8 

84 

10. 

2 

116 

0. 

5 

111100000 

0 

35 

-0. 

1 

-1 

1016. 

,2 

12. 

4 

85 

10. 

0 

109 

0. 

5 

111100000 

0 

40 

-0. 

1 

-1 

1016. 

2 

12. 

4 

85 

10. 

0 

112 

0. 

5 

111100000 

0 

45 

0. 

.1 

4 

1015. 

6 

12. 

8 

85 

10. 

4 

103 

0. 

5 

111100000 

0 

50 

0. 

,1 

7 

1015. 

,2 

12. 

6 

85 

10. 

2 

78 

0. 

6 

111100000 

0 

55 

0. 

,2 

8 

1015. 

1 

12. 

8 

84 

10. 

2 

68 

0. 

7 

111100000 

1 

0 

0. 

,1 

5 

1015. 

5 

12. 

6 

84 

10. 

0 

56 

0. 

6 

111100000 

1 

5 

0. 

,3 

6 

1015. 

,3 

12. 

,5 

85 

10. 

1 

42 

0. 

,5 

111100000 

1 

10 

0. 

,2 

5 

1015. 

,5 

12. 

,4 

84 

9. 

,8 

0 

0, 

,0 

111100000 

1 

15 

-0. 

.1 

5 

1015. 

,5 

12. 

,5 

85 

10. 

,1 

0 

0, 

,0 

111100000 

1 

20 

0, 

,1 

9 

1015. 

,0 

12. 

,6 

85 

10. 

,2 

0 

0, 

,0 

111100000 

1 

25 

0, 

,1 

12 

1014. 

,6 

12, 

,4 

84 

9. 

,8 

0 

0. 

,0 

111100000 

1 

30 

0, 

.4 

17 

1014. 

,0 

12, 

.1 

85 

9. 

,7 

0 

0, 

,0 

111100000 

1 

35 

0, 

,5 

22 

1013, 

.4 

12, 

,0 

85 

9. 

,6 

0 

0, 

,0 

111100000 

1 

40 

0, 

.7 

25 

1013. 

.0 

12, 

,0 

85 

9, 

,6 

0 

0, 

.0 

111100000 

1 

45 

1, 

.0 

35 

1011. 

.9 

11, 

,9 

86 

9. 

,7 

0 

0, 

.0 

111100000 

1 

50 

1, 

.0 

40 

1011, 

.3 

11, 

,8 

86 

9, 

,6 

328 

0, 

.5 

111100000 

1 

55 

1, 

.2 

46 

1010, 

.5 

11, 

,8 

86 

9, 

,6 

332 

0, 

,5 

111100000 

2 

0 

1, 

.2 

54 

1009, 

.5 

11, 

.7 

86 

9, 

.5 

332 

0, 

.5 

111100000 

2 

5 

1, 

.3 

60 

1008, 

.8 

11, 

.7 

87 

9, 

.7 

326 

0, 

.5 

111100000 

2 

10 

1 

.3 

64 

1008, 

.3 

11, 

.6 

87 

9 

.6 

320 

0 

,5 

111100000 

2 

15 

1 

.3 

73 

1007, 

.2 

11 

.5 

87 

9 

.5 

314 

0 

.5 

111100000 

2 

20 

1 

.2 

76 

1006 

.8 

11 

.5 

87 

9 

.5 

0 

0 

.0 

111100000 

2 

25 

1 

.1 

81 

1006 

.2 

11 

.4 

88 

9 

.5 

0 

0 

.0 

111100000 

2 

30 

0 

.9 

82 

1006 

.1 

11 

.4 

88 

9 

.5 

0 

0 

.0 

111100000 

2 

35 

1 

.0 

89 

1005 

.3 

11 

.3 

88 

9 

.4 

0 

0 

.0 

111100000 

2 

40 

1 

.0 

95 

1004 

.6 

11 

.3 

88 

9 

.4 

0 

0 

.0 

111100000 

2 

45 

1 

.0 

103 

1003 

.6 

11 

.2 

88 

9 

.3 

0 

0 

.0 

111100000 
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Sounding  program  REV  7 . 62  using  Omega 

Ship      :PT_SUR 

Location  :  36.77  N  121.90  W       3m 

Rejected  Sigma  stations:  a,b,c 

Phase  fitting  length  is  250  s  from 

Sounding  :  37 

RS-number:  309140455 

No  PTU  editing 

No  wind  editing 

Started  at:       23  MAY  95  21:06  GMT 


0  min  to  120  min 


Time  AscRate  Hgt/MSL  Pressure 
WndStat 
min   s     m/s       m      hPa 


Temp   RH    Dewp   Dir  Speed 


degC 


degC   deg   m/s 


0 

0 

0 

5 

0 

10 

0 

15 

0 

20 

0 

25 

0 

30 

0 

35 

0 

40 

0 

45 

0 

50 

0 

55 

1 

0 

1 

5 

1 

10 

1 

15 

1 

20 

1 

25 

1 

30 

1 

35 

1 

40 

1 

45 

1 

50 

1 

55 

2 

0 

2 

5 

2 

10 

2 

15 

2 

20 

2 

25 

2 

30 

2 

35 

2 

40 

2 

45 

2 

50 

2 

55 

0.0  3 

-461.8  ////// 

-277.1  ////// 

-173.2  ////// 

-138.6  ////// 

-106.6  ////// 

-92.4  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 

0.0  ////// 


104. 

3 

11. 

9 

93 

10. 

8 

236 

1.9 

1011. 

4 

12. 

9 

76 

8. 

8 

/// 

//// 

— CD--GH 

1009. 

0 

12. 

3 

80 

9. 

0 

190 

2.3 

--CD--GH 

1005. 

6 

12. 

0 

83 

9. 

3 

/// 

//// 

— CD--GH 

1002. 

9 

11. 

8 

84 

9. 

2 

201 

1.9 

— CD--GH 

998. 

9 

11. 

6 

85 

9. 

2 

/// 

//// 

— CD — GH 

997. 

0 

11. 

4 

86 

9. 

2 

244 

1.5 

— CD--GH 

994. 

8 

11. 

2 

87 

9. 

2 

/// 

//// 

— CD-FGH 

992. 

7 

11. 

0 

88 

9. 

1 

249 

1.6 

— CD-FGH 

991. 

2 

10. 

9 

89 

9. 

2 

/// 

//// 

— CD-FGH 

989. 

5 

10. 

7 

89 

9. 

0 

285 

2.2 

— CD-FGH 

987. 

0 

10. 

6 

90 

9. 

1 

/// 

//// 

— CD-FGH 

984. 

,5 

10. 

4 

91 

9. 

0 

260 

1.6 

— CD-FGH 

983. 

,3 

10. 

3 

92 

9. 

,1 

/// 

//// 

— CD-FGH 

981. 

,1 

10. 

1 

93 

9. 

1 

257 

1.5 

— CD-FGH 

979. 

,1 

10. 

0 

94 

9. 

,1 

/// 

//// 

— CD-FGH 

977. 

,0 

9. 

9 

95 

9. 

,2 

255 

1.6 

— CD-FGH 

974. 

,3 

9. 

7 

96 

9. 

,1 

/// 

//// 

— CD-FGH 

970. 

,6 

9. 

,6 

96 

9. 

,0 

296 

2.6 

— CD-FGH 

967, 

.8 

9. 

,5 

96 

8. 

,9 

/// 

//// 

— CD-FGH 

964, 

.7 

9. 

,3 

97 

8, 

.9 

304 

3.1 

--CD-FGH 

962, 

,7 

o 

.2 

96 

8, 

.6 

/// 

//// 

— CD-FGH 

959, 

.1 

9. 

.0 

96 

8, 

.4 

311 

3.9 

— CD-FGH 

956, 

.4 

9, 

.0 

97 

8, 

.6 

/// 

//// 

— CD-FGH 

953 

.7 

8, 

.8 

97 

8, 

.4 

311 

3.9 

— CD-FGH 

950 

.3 

8, 

.7 

97 

8 

.3 

/// 

//// 

--CD-FGH 

946 

.7 

9, 

.1 

98 

8 

.8 

314 

4.1 

--CD-FGH 

944 

.2 

10, 

.0 

99 

9 

.9 

/// 

//// 

— CD-FGH 

942 

.1 

10, 

.3 

98 

10 

.0 

316 

4.6 

— CD-FGH 

940 

.3 

10 

.5 

94 

9 

.6 

/// 

//// 

— CD-FGH 

938 

.5 

10 

.6 

92 

9 

.4 

319 

5.0 

— CD-FGH 

936 

.2 

11 

.2 

95 

10 

.5 

/// 

//// 

— CD-FGH 

933 

.4 

12 

.5 

95 

11 

.8 

319 

4.8 

— CD-FGH 

932 

.0 

13 

.0 

93 

11 

.9 

/// 

//// 

— CD-FGH 

930 

.0 

13 

.2 

91 

11 

.8 

322 

5.6 

— CD-FGH 

927 

.4 

13 

.3 

88 

11 

.4 

/// 

//// 

— CD-FGH 
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APPENDIX  J.  MATLAB  CODE  FOR  COMPUTATION 


A.         GPS,  ELEVATION  AND  AZIMUTH  ANGLE  CALCULATION 

%  GPS  &  elevation  &  azimuth  angle  caculation 
%  File  name  gps.m 

ship_longitude=  121.81; 
ship_latitude=36.8 1 ; 
sensor_longitude=  121.8; 
sensor_latitude=3  6. 8 ; 
heading=120; 

dx=(ship_longitude-sensor_longitude)*47.995; 
dy =(ship_latitude-sensor_latitude)  *  5  9 .9  7 ; 

a=(atan(abs(dx/dy)))*57.3 

ifdy>0 
ifheading>(180-a) 

bearing=(heading- 1 80+a); 
else 

bearing=l  80-a-heading 
end 
else 
if  heading<a 

bearing=a-heading; 
else 

bearing=heading-a; 
end 
end 

if  bearing<90 

phi=(90-bearing)/57.3 
else 

phi=(180-bearing)/57.3 
end 

range_in_nmi=sqrt(dx  *  dx+dy  *  dy ) 
range_in_km=(range_in_nmi)*  1 .852 
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theta=atan(0 .  003  2/range_in_nmi) 


1=41.5; 

%ship  length 

w=9.75; 

%ship  width 

h=8.8; 

%ship  height 

Area=l  *  w*  sin(theta)+h*  w*  cos(theta)*  sin(phi)+h*  1  *  cos(theta)  *  cos(phi) 
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B.         MRT  &MDT  vs  Spatial  Frequency 

%  MRT  &  MDT  vs  Spatial  Frequency 
%  File  name  drtfre 

SNRT=2.5 ;  %  Shumaker  Example  8- 1 3 

NET=0.125; 

dx=0.25; 

dy=0.25; 

L=7; 

Te=0.1; 

Fr=30; 

Nos=l; 

Nss=l; 

rb=0.335; 

rs=0.456; 

v=0.000001:0.25:3.000001; 

MTF=[1  0.89  0.79  0.64  0.53  0.43  0.35  0.27  0.2  0.15  0.1  0.07  0.04]; 

omigaT=l  ./(4*v.*v); 
ROx=(l+(2*2*v.*v*rb.*rb)).A(-0.5); 

MRT=2*SNRT*NET*(ROx.A0.5).*((v.*v*dx*dy/L).A0.5).*((Te*Fr*Nos*Nss)A(- 
0.5))./MTF; 

MDT_l=SNRT*NET*(omigaT+rs*rs)*((dx*dy)A0.5); 
MDT_2=omigaT.*((((pi/4)*(rs*rs+rb*rb+omigaT)*Te*Fr*Nos*Nss)).A0.5); 
MDT=MDT_1  ./MDT_2; 

plot(v,MDT,v,MRT); 

title('MRT  &  MDT  vs  Spatial  Frequency'); 

xlabel('Spatial  Frequency'); 

ylabel('Temperature'); 

gtext('MRT'); 

gtext('MDT'); 

grid; 
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C.         MRT,  MDT  &  AT  vs  Range  for  Detection 

%  MRT  &  MDT  vs  &  Range 
%  N=l  for  Detection 
%  File  name  r  n  1   f 


SNRT=2.5; 

NET=0.125; 

dx=0.25; 

dy=0.25; 

L=7; 

Te=0.1; 

Fr=30; 

Nos=l; 

Nss=l; 

rb=0.335; 

rs=0.456: 


%  Shumaker  Example  8-13 


N=l; 

At=2 14.45 13; 

Dc=(At)A(0.5); 

X=Dc/N; 


%  physical  dimension 


dt=3.387144; 
u=0.0002132; 

v=0.000001:0.25:3.000001; 

MTF=[1  0.89  0.79  0.64  0.53  0.43  0.35  0.27  0.2  0.15  0.1  0.07  0.04]; 

Range=v*2000*X; 

omigaT=l  ./(4*v.*v); 

ROx=(l+(2*2*v.*v*rb.*rb)).A(-0.5); 

tau=exp(-u.*Range); 
dtapp=dt*tau; 

MRT=2*SNRT*NET*(ROx.A0.5).*((v.*v*dx*dy/L).A0.5).*((Te*Fr*Nos*Nss)A(- 
0.5))./MTF; 

MDT_l=SNRT*NET*(omigaT+rs*rs)*((dx*dy)A0.5); 
MDT_2=omigaT.*((((pi/4)*(rs*rs+rb*rb+omigaT)*Te*Fr*Nos*Nss)).A0.5); 
MDT=MDT  l./MDT  2; 
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plot(Range,MDT,Range,MRT,Range,dtapp); 

title("MRT,  MDT  &  dT  vs  Range  for  Detection'); 

xlabel('Range  (m)'); 

ylabel(Temperature'); 

gtext('MRT'); 

gtext('MDT*); 

grid; 
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D.         MRT  &  AT  vs  Range  for  Classificaton 

%  MRT  &  dt  vs  Range 
%  N=8  for  Classification 
%  File  name  r  n  8  f 


SNRT=2.5; 

NET-0.125; 

dx=0.25; 

dy=0.25; 

L=7; 

Te=0.1; 

Fr=30; 

Nos=l; 

Nss=l; 

rb=0.335; 

rs=0.456; 


%  Shumaker  Example  8-13 


N=8; 

At=214.4513; 

Dc=(At)A(0.5); 

X=Dc/N: 


%  physical  dimension 


dt=3.387144; 
u=0.0002132; 

v=0.000001:0.25:3.000001; 

MTF=[1  0.89  0.79  0.64  0.53  0.43  0.35  0.27  0.2  0.15  0.1  0.07  0.04]; 

Range=v*2000*X; 

omigaT=l./(4*v.*v); 

ROx=(l+(2*2*v.*v*rb.*rb)).A(-0.5); 

tau=exp(-u.  *  Range) ; 
dtapp=dt*tau; 

MRT=2*SNRT*NET*(ROx.A0.5).*((v.*v*dx*dy/L).A0.5).*((Te*Fr*Nos*Nss)A(- 
0.5))./MTF; 

MDT_l=SNRT*NET*(omigaT+rs*rs)*((dx*dy)A0.5); 
MDT_2=omigaT.*((((pi/4)*(rs*rs+rb*rb+omigaT)*Te*Fr*Nos*Nss)).A0.5); 
MDT=MDT  l./MDT  2: 
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plot(Range,MRT,Range,dtapp); 

title('MRT  &  dT  vs  Range  for  Classification'); 

xlabel('Range  (m)'); 

ylabel(Temperature'); 

gtext(*MRT'); 

grid; 


143 


E. 


MRTD  AND  AT  VS  RANGE  FOR  IDENTIFICATION 


%  MRT  &  dT  vs  Range 
%  N=12.8  for  Identification 
%  File  name  r  n  128  f 


SNRT=2.5; 

NET=0.125; 

dx=0.25; 

dy=0.25; 

L=7; 

Te=0.1; 

Fr=30; 

Nos=l; 

Nss=l; 

rb=0.335; 

rs=0.456; 


%  Shumaker  Example  8-13 


N=12.8; 

At=2 14.45 13; 
Dc=(At)A(0.5); 

X=Dc/N; 


%  physical  dimension 


dt=3.387144; 
u=0.0002132 

v=0.000001:0.25:3.000001; 

MTF=[1  0.89  0.79  0.64  0.53  0.43  0.35  0.27  0.2  0.15  0.1  0.07  0.04]; 

Range=v*2000*X; 

omigaT=l  ./(4*v.*v); 

ROx=(l+(2*2*v.*v*rb.*rb)).A(-0.5); 

tau=exp(-u.  *  Range); 
dtapp=dt*tau; 

MRT=2*SNRT*NET*(ROx.A0.5).*((v.*v*dx*dy/L).A0.5).*((Te*Fr*Nos*Nss)A(- 
0.5))./MTF; 

MDT_l=SNRT*NET*(omigaT+rs*rs)*((dx*dy)A0.5); 
MDT_2=omigaT.*((((pi/4)*(rs*rs+rb*rb+omigaT)*Te*Fr*Nos*Nss)).A0.5); 
MDT=MDT  1./MDT  2; 
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plot(Range,MRT,Range,dtapp); 

title('MRT  &  dT  vs  Range  for  Identification'); 

xlabel('Range  (m)'); 

ylabel('Temperature'); 

gtextCMRT*); 

grid; 
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